PROCEEDINGS 


THE 


AMERICAN SOCIETY CIVIL ENGINEERS 


VOL. OCTOBER, 1947 No. 


TECHNICAL PAPERS 


AND 


DISCUSSIONS 


Published monthly, except July and August, Prince and Lemon Streets, Lancaster, Pa., the American 
Society Civil Engineers. Editorial and General Offices West Thirty-ninth Street, New York 18, 
Reprints from this publication may made condition thai the full title paper, 
name author, page reference, and date publication the Society are given. 


Entered Second-Class Matter, September 23, 1937, the Post Office Lancaster, Pa., under the Act 
March 1879. Acceptance for mailing special rate postage provided for Section 1103, 
Act October 1917, authorized July 1918. 


(if entered before January $8.00 per annum Price $1.00 per copy 


Copyright, 1947, the AMERICAN ENGINEERS 
Printed the United States America 


4 
q 
| 
| 
q 


Current Papers aND Discussions 


Published 
Machis, Alfred. Observations Grouting Sands and Gravels..Nov., 1946 
Symposium: Cleaning and Grouting Limestone Foundations, Tennessee 
Matthes, Gerard Mississippi River Jan., 1947 
Muldrow, Forecasting Productivity Irrigable Lands.............. Feb., 1947 
Freudenthal, Alfred Reflections Standard Specifications for Structural 
Blaisdell, Fred Development and Hydraulic Design, Saint Anthony Falls 
DuVal, Miles The Marine Operating Problems, Panama Canal, and the 
Claybourn, Sea Level Plan for Panama Feb., 1947 
Hrennikoff, Alexander. Theory Inelastic Bending with Reference Moment 
Shih-Chang. Beam Deflections Second and Third Moments......... Mar., 1947 
Ryder, Frederick Rational Explanation Column Behavior........... Mar., 1947 
Rathbun, Charles, and Cunningham, Continuous Frame Analysis 
Sturm, Stability Thin Cylindrical Shells Torsion............... Apr., 1947 
Report: Water Supply Engineering: Report the Committee the Sanitary 
Engineering Division for the Two Years Ending December 31, 1946....... Apr., 1947 
Committee the San Francisco (Calif.) Section, ASCE, Timber Test Pro- 
gram. Tests Timber Structures from Golden Gate International Expo- 
Westergaard, New Formulas for Stresses Concrete Pavements Air- 
N., and Sorkin, Josef. Mississippi River Bridgé Dubuque, 
Casagrande, Arthur. Classification and Identification Soils............... June, 1947 
Barron, Reginald Consolidation Fine-Grained Soils Drain une, 1947 
Kazmann, River Infiltration Source Ground Water Supply..June, 1947 
Spangler, M.G. Underground Conduits—An Appraisal Modern Research..June, 1947 
Jens, Stifel Drainage Airport Surfaces—Some Design Considera- 
Peterson, Dean F., Jr. Lines for Continuous Structures Geo- 
Boissonnault, Frank Estimating Data for Reservoir Gates.............. Sept., 1947 
Coddington, F., and Marshall, Least Squares Adjustment Tri- 
angulation Net Between Geodetic Sept., 1947 
Kellogg, Investigation Drainage Rates Affecting Stability Earth 


registered for special extension time. 
Publication closing discussion pending. 


Discussion 
closes 


Closed* 


Closed* 


Closed* 


Closed* 


Nov. 1947 


Nov. 1947 


Nov. 1947 


Uncertain 


Uncertain 


Dec. 
Dec. 


Dec. 


Dec. 


Dec. 


Dec. 


Dec. 


Dec. 


Jan. 


Jan. 
Jan. 
Jan. 
Jan. 


Feb. 


Feb. 
Feb. 


Feb. 


Feb. 


closing dates herein published are final except when names prospective discussers are 


1947 


1947 


1947 


1947 


1947 


1947 


1947 


1947 


1948 


1948 
1948 
1948 
1948 


1948 


1948 
1948 


1948 


1948 


Sul 


The 


Stal 


For 


Clas 


Mis: 


The 


PROCEEDINGS 


CONTENTS FOR OCTOBER, 1947 


PAPERS 


Thin-Shelled Domes Loaded Eccentrically. 
Walter Voss, Dean Peabody, Jr., Howard Staley, and Albert Dietz. 1173 


Subsidence the Terminal Island-Long Beach Area, California. 
Frederic Harris and Eugene Harlow 


Variation Coefficients Simultaneous Linear Equations. 
Bernard Weiner 


Application Geology Tunneling Problems. 
Ernest Wahlstrom 


The General Principles Highway Transition Curve Design. 
John Joseph Leeming 


Analysis Stepped-Column Mill Bents. 
Daniel Ling 


DISCUSSIONS 


Theory Inelastic Bending with Reference Limit Design. 
Michael Horne 


Stability Thin Cylindrical Shells Torsion. 
Batdorf and Manuel Stein 


Forecasting Productivity Irrigable Lands. 
Lewis 


Friction Coefficients Large Tunnel. 
Tapley 


Classification and Identification Soils. 
James Stratton, and Donald Belcher 


Mississippi River Bridge Dubuque, Iowa. 


1197 
1219 
1245 
1302 
1307 
48 
are 
1319 


CONTENTS FOR OCTOBER, 1947 (Continued) 


Rational Explanation Column Behavior. 
Tests Timber Structures from Golden Gate International Exposition. 
James Carr, Jr., Howard Hansen, and 1323 


list “Current Papers and Discussions” may found the page preceding 
the table contents 


The Society not responsible for any made opinion expressed 
its publications 


HINTS AUTHORS 


Those who are planning papers for submission PROCEEDINGS will expedite 
Committee action measurably first studying the standard instructions 
style, content, and format. For reprints, address the Manager, Technica 
Publications, 39th St., New York 18, 


4 


| 
K 
Ww 
cus 
| 
q 
| 


AMERICAN SOCIETY CIVIL ENGINEERS 
Founded November 1852 


PAPERS 


THIN-SHELLED DOMES LOADED 
ECCENTRICALLY 


WALTER Assoc. ASCE, DEAN PEABODY, 
ASCE, AND HOWARD AND ALBERT 
Assoc. MEMBERS, ASCE 


SYNOPSIS 


shallow thin-shelled model dome plaster Paris was subjected 
concentrated load successively placed between the springing and the crown. 
Strains were measured near-by and distant stations the upper surface. 
Compression and transverse tests the plaster gave data concerning moduli 
elasticity and Poisson’s ratio whereby stresses were computed from the re- 
corded strains. Sand loads were distributed uniformly over eighth seg- 
ment and over quarter segment study the transition from concentrated 
load uniformly distributed load over the entire surface. 

mathematical analysis was made determine the stresses caused 
concentrated load applied the crown shallow thin-shelled spherical dome. 
Stresses any magnitude are essentially local ones, this analysis applies 
approximately eccentric load applied some distance from the springing 
any hole the dome. far the test one model dome can authorize, 
this series tests seems justify the use the mathematical analysis com- 
pute maximum stresses shallow dome. 


INTRODUCTION 


The absence information the effects eccentric concentrated loads 
thin-shelled domes with small rise prompted these tests. Domes this type 
were already existence and were being built without accurate data regard 
the stresses carried load concentrations. 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted March 1948. 


Head, Dept. Bldg. Eng. and Constr., Mass. Inst. Tech., and Consultant, Cambridge, Mass. 
Associate Prof., Dept. Bldg. Eng. and Constr., Mass. Inst. Tech., Cambridge, Mass. 

Asst. Prof., Dept. Bldg. Eng. and Constr., Mass. Inst. Tech., Cambridge, Mass. 

Asst. Prof., Dept. Bldg. Eng. and Constr., Mass. Inst. Tech., Cambridge, Mass. 
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check laboratory and theoretical studies against loading data from ac- 
tual structures, was decided build plaster model large dome, actually 
erected, scale one eight. The base the spherical dome was 
diameter; its rise was d); and the thickness the shell The 
radius the sphere was 8.5 ft. 

This exploratory test was made 1941, before the mathematical analysis 
had been initiated; and was not considered wise that time assume that 
eccentric concentrated load would produce measurable strains only the 
surface close the area load application. Consequently, sixty-four gage 
stations were located the dome surface with two wires connected each 
station. The assumption was made that the axes principal strains for highly 
strained stations would at, near, the meridional and circumferential axes 
the dome. mathematical analysis Eric confirms this as- 
sumption (see Appendix and indicates that any divergence test data 
would probably slight. would necessary add sixty-four more wired 
connections obtain third strain each station complete rosette. 

The inner form the dome was not dropped far enough permit strain 
readings the under surface. attempt was made read deflections 
0.001 in. with transit, but deflections this magnitude were recorded. 
subsequent tests are made, desirable that the number stations re- 
duced that three strains each station may taken complete rosette, 
that stations located both the upper and lower surfaces near the area 
load application, and that accurate measurement deflections attempted 
near the load application. 

Strains were measured the meridional and circumferential directions 
each station. After the completion the concentrated load tests, the dis- 
tributed sand load tests were made, the sand load remaining several days 
the dome while the time effect the strains was measured. Finally, the dome 
was tested failure concentrated load applied about three tenths the 
meridian length from the springing. 

letter symbols used this paper are defined they occur, 
and are assembled for convenience reference Appendix II. 


PART EXPERIMENTAL DETERMINATION STRAINS 
GENERAL 


objective these tests was study the effect con- 
centrated loads spherical dome. Test data models full-sized domes 
are lacking for such loads, especially for loads that are not applied the cen- 
tral axis. nearly can determined, theoretical analyses were not 
existence prior the work done Mr. Test data for eccentric 
loadings will valuable for establishing the extent the surface strained and 
the comparative magnitude these strains, and for comparison with the 
mathematical analysis. 

Eccentric distributed loading desirable supply the transition from 
eccentric concentrated load symmetrical uniformly distributed load. 


and Small Displacements Shallow Spherical Eric Reissner, Journal Mathe- 
matics and Physics, February, 1946, 80. 
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Construction the rigid wooden ring was mounted frame 
even plane. recess left the top inner part this ring acted sup- 
port and also took the reaction the model when loaded. The bottom 
this recess acted guide for the templet during the operation fabricating 
the base mold which the dome 
shell was cast. The top the 
recess was used guide for 
forming the dome shell. The 
base mold and dome shell were 
formed rotating metal 
templet about central shaft. 
testing, the base mold was 
lowered jacks. 

Strain Gages.—Electrical re- 
sistance wire gages the type 
Ruge, Assoc. ASCE, and 
deForest the Massachu- 
setts Institute Technology Line 
Cambridge were adopted best 
meeting the requirement for pre- 
cision over short gage length. 
Gages were placed the sixty- 
four stations shown Fig. 
each station, strains were 
measured the circumferential 
and meridional directions. The 
central selector panel which al- 
lowed any one gage read will. the initial tests, the gages were 
straight wires, approximately in. long, and were cemented the surface, 
shown Fig. hairpin gages were used at, near, the loaded 
stations. These hairpin gages measured the strain for 1.5-in. length. 
autumn 1940 the shorter gages, now use, were not sale. 


PROCEDURE 


Concentrated most these tests concentrated load 
was used. Lead blocks, in. square, rested four soft rubber pads, each 
in. thick and in. square. The pads were separated clearance, 
that was in. out-to-out for the pads, whereas the actual loaded area was 
6.25 in. The loads were applied the meridian other points be- 
tween the 45° and 315° meridians. 

The test failure was made applying concentrated load Station 5,6 
the meridian (see Fig. 1). this case, lever system was employed 
which consisted arm swinging about pivot through the center the 
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dome. This arm had load platform the far end, and rested leg bear- 
ing shoe placed the rubber pads (Fig. 2). 

Sand Loads.—A form plywood and sheet metal rested the dome shell 
and confined the sand definite area (Fig. 3). Dry sand was spread uni- 
formly over the segment until the desired weight was obtained. 


Compression tests were made plaster cylinders, in. diameter and 
in. high. Gages similar those used the dome were installed give longi- 


tudinal and circumferential strains. From the data, the modulus elasticity 
and Poisson’s ratio were computed. 


/ / \ 
Strains Compression 
/ 
| ——-—-— Tension 


Strain Equals Contour 
Value Times 
Inches per Inch 


Transverse tests were also made small pieces the broken dome. 
Assuming uniformly varying stresses and strains, fiber stresses and the modulus 
rupture were computed, well the modulus elasticity. For each speci- 
men plot the load-strain data can represented two straight lines 
slightly different slopes. However, the comparison moduli elasticity with 


8 ANN SZY 7 
\ 
\ 4 \ 
x \ 1 
\ 
18 WS) 


1178 THIN-SHELLED DOMES Papers 


the compression moduli (presented subsequently Table has been made 
the stress varies tenths the ultimate. 


EXPERIMENTAL AND RESULTS 


The complete report which this paper based (including Mr. Reissner’s 
complete mathematical treatment, which reported briefly Appendix 
can consulted Engineering Societies well the library 
Massachusetts Institute Technology. 

Concentrated strains measured for given position the con- 
centrated load were assembled tables such Tables 1(a) and 1(b). These 
strain data were plotted show the variation along meridian and along cir- 


(The Number Times 10-5 the Strain Inches Per Inch. 


Station Distance (in Inches) Load (in 


Cols. 

age 


See also Stations Nos. 62, inclusive, the continuation Table Load, Values 
the 26-in. distance, the 315° line and the 45° line, were identical. The average strain Col. the aver- 


cumferences. The contour plots Fig. show equal compression strain 
full lines, and equal tension dashed lines. Table gives the data for the 
test failure concentrated load Station 5,6, and Fig. shows the type 
failure. 


PART II. DISCUSSION RESULTS 


strains were measured the upper surface the dome for 
various concentrated loads applied gage station, applied measured 
distances circumference through the gage station. Thus, was possible 
record, with minimum stations, the strains definite distances from the 
loading point. For certain loadings several tests were made the same cir- 
cumference. these cases the readings have been averaged. Most the 
loads were constant lb, but certain cases loads lb, 55.0 lb, 55.1 


West 39th St., New York, 


Col. 


age 


sti 
lo: 


thi 
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and 56.1 were used. The corresponding strain data were revised com- 
mon base load lb, assuming strains vary directly with the low loads 
used. 

was not feasible locate gage points the under surface the dome 
because the form was dropped only few inches and only the gage stations 
near the outer haunch could have been reached. 

The averaged strain data are also shown the contour plots, using contour 
lines connect points the same strain value. These contours (see Fig. 
were located scaling the plotted strain values. necessity many these 
contour lines were sketched connect rather widely dispersed points 
equal strain. However, close the loaded areas, they are accurate. 


Elongation Negative and Shortening Positive) 


Pounds) from Station Sta- 


tion 
No. 


(16) 


40.5 


Ols. ols. ols. ols. ois. ois. ois, ois. 


Col. the strain readings Col. revised common base corresponding 40-lb load. The strains 
age of the two. 


The common characteristic these plots the large compressive strain 
immediately under the loaded area both meridional and circumferential 
directions. The strains decrease rapidly outside the loaded area; farther away 
small tensile strains usually occur, often followed equally small compressive 
strains. The deformations suggest those observed one presses finger 
against the surface inflated toy balloon. Except for the area close the 
loaded point, the magnitude the strains very small. Near the loaded 
area, the meridional contours are approximately circular whereas the cireum- 
ferential contours are The reader should note that 
and refer the central dome axis and not axis through 
the center the load. 

The strains are comparatively small for load applied Station 1,2 near 
the outer ring, haunch (see Fig. 5). stiffer haunch influences the strain 
distribution. the position application the 40-lb load moved the 
meridian toward the crown the dome, the strains increase maximum 


: 
(in 
-10 
ues In 
q 
7 
r 
’ 


1180 THIN-SHELLED DOMES Papers 


Station 21,22, which was the station nearest the hole the crown, the in- 
crease being about 40% more than the maxima Station numerical 
values the maximum strains the meridional and circumferential directions 
the same station are approximately equal. 

stress-strain plots the compression and transverse tests 
indicate that the plaster behaves elastic material for stresses less than 
persqin. strain data stresses, average values the modu- 
lus elasticity (1,770,000 per in.) and Poisson’s ratio (0.20) are 
used. Knowing the meridional strain and the circumferential strain 


under the load, the corresponding stresses and can computed the 
elastic relations: 


and 

TABLE 1.—(Continued) Using the averaged maxi- 
mum strains under the load, the 
maximum stresses for 40-lb 
load applied the successive 
Station gage stations are given Table 
These maximum stresses in- 

crease the load applied 
any given station, the maxi- 


garded asequal. The computed 
+02 stresses Station 21,22 are 
+0.4 about 16% the compressive 
strength the plaster cylinders 

—0.2 


that the load applied relatively 


pounds. Values Col. are the strain distant from boundary 


readings Col. revised common base corresponding hole the shell. Ap- 
40-lb load. 


a 
ve 
0 
FO25 
owes 


NOLLVIG LV Of 40 AGVOT V NOILVIUVA NIVUIG SISAIVNY—'G 


° 
| 
Bury : 
™ 


: 
vs 
# 
| 
7] 
l- 
a 
q 
a 7 
q 
i 
ue 7 
d 


1182 THIN-SHELLED DOMES Papers 


normal stresses the meridional and circumferential directions the upper and 
lower surfaces the shell due direct stress and bending. The equations are 
then applied structure with the dimensions the test dome. The theo- 
retical analysis compared with test results Appendix (see heading, 
where the loaded circular area assumed 6.25 
in. believed that these results are comparable because the actual 


Compression Tension 


180° 


Compression 
Compression 


Compression 


Tension 


Tension 
Tension 


45° 


Compression Tension (Multiply 


VARIATION ALONG THE MERIDIANS FOR SAND 


loaded area during the tests was through four thin rubber pads, each 1.25 in. 
square. Therefore, the actual loaded test area was 6.25 in. Using the 
numerical results thus determined, the normal stresses, direct and bending, 
the upper surface equal: 
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and 
Compression Tension 


180° 


Compression 


Compression 


Compression 


70° 


Tension 
Tension 
Tension 


MERIDIONAL STRAINS 
Tension 


value close the average the test stresses Station 21,22. The stresses 


the upper surface are larger than the resultant stresses the under surface. 
appears that the theoretical analysis approximates closely the 
stresses given the test data, even though Station 21,22 about 6.5 in. 
from the edge hole, approximately in. diameter, the crown the 
test dome. The results from only one test dome are not conclusive, but they 
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establish the possibility that the theoretical analysis may used with safety 
design shallow domes subjected concentrated load. 


Papers 


Fig. 6(a) examined, will noted that the curves stress variation 


are similar the plots strain variation from the test data. 
along the meridian the theoreti- 
cal stresses should equal zero 
about 4.95 in. out from the center 
the loaded area, for meridio- 
nal stresses, and about 7.15 in. 
out, for circumferential stresses. 
The strain data the tests indi- 


TABLE FAILURE, WITH 
CoNCENTRATED Loap 3-IN. 
APPLIED LEVER 
System 5,6 


(The Number Times the Strain Inches per 
Inch. Elongation Negative and 
Shortening Positive) 


Measuring 


cate that zero stresses occur, 


Line No. Load Station Station 

(1) (2) (3) (4) 
96.5 28.35 31.0 
6.. 181.7 73.9 79.4 
280 +110.0 +118.6 
178.3 76.1 84.2 
236 +108.9 +117.0 
264.5 +127.1 +137.9 
293 +160.3 +165.5 


general, about in. out for both 
stresses. 
the theoretical circular 
area and the actual four squares 
probably explains the noncon- 
formance. 
compressive stresses the up- 
per surface have reduced zero 
distance beyond the loaded 
boundary about twice the 
diameter, width, the loaded 
area. 


The 


either case, the 


failure was made applying concentrated load Station 5,6 (Fig.1). The 
load-strain data are recorded Table These data, with moduli elasticity 


and values Poisson’s ratio estimated from the 
plaster cylinder and transverse tests, were used 
compute stresses. The results are given Table 

This test was made several weeks after the 
concentrated load test Station 5,6; 
and the strains due the load 37.6 lb, even 
for the reloading (see Table are about half 
those recorded earlier Station 5,6. would 
appear that the strains for the lower loads are 
smaller than would expected from the previous 
tests. During the period testing, the plaster 
had dried thoroughly. However, with the great- 
est load applied (293 lb), the computed stresses 
approach the compressive strength the moist 
cylinders, giving about 94% the cylinder 


TABLE 3.—Maximum 
STRESSES 

PounDs PER SQUARE 
CREATED 


Station 
(Fig. 
(1) (2) 
164 167 
437 436 
462 476 
545¢ 5452 
596 631 


meridional and circumferential 
directions. 


strength for load that was 91% the load failure. 
The stresses computed the mathematical analysis Appendix ap- 
proximated the stresses obtained from the measured strains only for the load 


applied Station 21,22. For the 40-lb load applied Station 5,6, the average 
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stress from the strain data was 469 per in., 77.8% the theoretical 
stress 603 per Using this percentage reduction the theoretical 
stress applied the second loading, comparison afforded Table 
the stresses from strain data (Col. and the corrected theoretical stresses. 


(Ib) ratio Average cylinder age 
ional, erential, 
(1) (2) (3) (4) (6) (7) (8) (9) 


(a) 


1.77 0.20 2.4 11.7 
1.77 0.20 380 352 366 9.5 20.7 
1.77 0.20 676 636 656 30.0 
1.77 0.20 995 925 960 25.0 38.7 
0.20 1,300 1,260 1,280 33.3 47.5 
1.76 0.21 1,746 1,670 1,708 44.3 56.4 
1.76 0.23 2,164 2,096 2,130 55.3 65.5 
1.69 0.24 2,600 2,484 2,542 66.1 87.0 
(b) Loap RELEASED AND RELOADED 
1.77 0.20 274 254 264 441 6.9 
1.76 0.22 1,858 1,740 1,799 2,090 46.8 5.5 
1.71 0.23 2,564 2,450 2,507 2,770 65.2 73.3 
1.67 0.25 3,020 2,880 2,950 3,190 76.6 
1.65 0.26 3,662 3,598 3,630 3,440 94.3 
The values Col. were computed the formula: 
0.9415 
0.25? (4a) 


The theoretical stress 441 per in. for load 37.6 the value that 
would expected from the strain data the initial tests Station 
the dome approached failure this same ratio (0.778) test data theoretical 
analysis was approached for the strain data Table The dome failed 
close the load predicted stress computations for the upper surface. 
However, the computed stresses for the lower surface are tensile and 


For load 293 lb, this tensile stress equals 0.778 0.5539 
293 2,020 per in., compared with the compressive value 
3,440 per in. the upper surface (Col. Table The transverse 
tests Table indicate modulus rupture for the plaster about 750 
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per in., the failure being tensile. would seem that the effect shrinkage 
and plastic flow such that the computed tensile stresses the under surface 
are not realized. Until greater experience gained, the high computed tensile 
stresses are added incentive use conservative allowable stresses and 


Gage No. (Compression) 
P 


Specimen Widths 
number 


1.98 
2.00 
Gage No. (Tension) 2.02 


2.205 4.410 6.615 8.820 11.025 15.435 19.845 22.050 


+4.74 +9.61 +14.29 +18.88 +33.88 +33.40 +43.5 +48.6 


—6.15 —12.42 —18.70 —24.94 —30.94 —45.20 
+6.06 +12.25 +18.16 +24.07 +30.40 +42.85 


(c) Srress, INcH, ASSUMING THAT 
Srress AND VARY UNIFORMLY 


75.8 151.5 227.5 300 379.5 530.5 682 
80.8 161.9 242.9 320 404.0 566.0 728 
106.2 212.5 318.8 425 531.0 744.0 


THE NuMBER MULTIPLIED 


1.55 1.535 1.515 1.485 
1.55 1.54 1.535 
1.73 1.69 


disregard such stress reduction factors for position the load the value 
0.778 used Eq. 4a. Tests observe variations modulus elasticity 
are presented Table 

Continued time-load test was made with applied 
Station 7,8 for days. During this interval the compressive strains directly 
under the loaded area gave increased values. The strains the other stations 
tended become tensile, this effect being most marked from day days. 
Between the initial and the 24-hour readings, the maximum increase com- 
pressive strain occurred Station 7,8, followed slight steady decrease for 
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the next From days day few the other stations also showed 
increase compressive strain which was overcome the decrease for the 


next days. 


loads about one eighth the ultimate load were 
applied stations meridian varying position from close the outer 


haunch close the central hole. every 
test the greatest stresses occur the upper sur- 
face the shell directly under the load, and they 
increase the load moved toward the crown. 
The stresses for load close the crown check 


TABLE 
AVERAGE TESTS FOR 
ELASTICITY 
PoUNDS PER SQUARE 
(NUMBER 


closely those obtained mathematical Times 10°) 
(see Appendix I). The compressive stresses 
particles the upper surface decrease rapidly 


the distance increases away from the loaded area. 


Tan- 

The stresses located twice the width the loaded Secant 
area from the loaded boundary have ceased 

1.60 1.71 1.59 


concentrated load remains place for 
several days, the strains increase about 15% un- 
der the load during the first day. Thereafter, 


strength. 


there slow decrease; days the strains are 


about greater than the initial readings. Elsewhere the dome the 
strains showed small elongations due plastic flow. 


DISTRIBUTED 


Eighth Segment Loaded.—The sand load 267 eighth segment cor- 
gives strain data the time load application for continued loading 1.7 
days. Fig. shows the plotted strain data the time load application. 
The sand load was centered the meridian, and all other recorded stations 
were outside the loaded area. The initial circumferential strains (Fig. 
were relatively large along the meridian with maximum values Stations 
and about the magnitude those recorded for the concentrated load 
Station 21,22. Elsewhere the dome the circumferential strains were 
small tensile magnitudes, except for the innermost Station the 90° 
meridian and three stations the 315° meridian. The meridional strains 
(Fig. 6(b)) were also compressive the meridian, but the maximum values 
occurred nearer the crown the dome Stations and 22. Elsewhere the 
strains were very low tensile values. Using these strain data, with modulus 
elasticity 1.77 per in., and Poisson ratio 0.20, the com- 
puted stresses along the meridian are given Table 8(a). 

None these stresses, caused load 267 lb, exceed those obtained 
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(The Number Times the Strain Inches Per Inch. 
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the concentrated load Station 21,22. the meridional stress, 
each station plotted according its position the meridian, the 
individual values not deviate greatly from their average. This not true 
for the circumferential stresses, which greatly exceed their average from 
lucation close Station 1,2 (and the 

haunch) position about halfway TABLE 


the meridian; beyond this point the cir- RIDIAN PER SQUARE 


cumferential stresses are much lower INcH) For SAND 


Continued loading for 0.85 day gives SEGMENT 
along the central meridian. The 
amount depression had decreased but Station SEGMENT 
the general shape was unchanged, since 
the stations greater strain had pro- 
portionately greater reductions. Else- (1) (2) (3) (4) 


sion under the load was reduced, whereas 
elsewhere the dome shell bulged outward. 

the interval from 0.85 1.7 days there was further readjustment. The 
circumferential compressive strains the meridian again decreased roughly 
proportion the strain value, whereas the meridional strains increased 
Stations and with considerable decrease Stations and near the 
haunch. Elsewhere the dome, outside the loaded area, the change strain 
was largely reduction tensile strain, that about half the final strains 
were compressive. appears that the depression under the load was again 
reduced; elsewhere the previous bulge was also decreased. 

The total effect continued loading from 1.7 days was reduce the 
compressive strains under the load the meridian. The circumferential 
strains were reduced from 38% 94%; and the meridional, from 20% 68%. 
The greater strains had the greater percentage reduction that the strains 
1.7 days tended approach common value more closely. The stations out- 
side the loaded area, for the most part, had net shortenings that compressive 
strains were increased and tensile strains were reduced. the twenty-five 
circumferential stations outside the loaded area, twenty one recorded small 
initial tensile strains less than 1.55 upon application the load. 
The other four were close the hole the crown the dome. the end 
1.7 days eleven had tensile, twelve compressive, and two zero strains, com- 
pared with the unloaded dome. Similarly, for the twenty-six meridional sta- 
tions outside the loaded area, the initial strains were: Nineteen tensile, six 
compressive, and one zero. Four the compressive strains were the 180° 
meridian, the other two were near the haunch. the end 1.7 days these 
twenty-six stations divided into twelve with tensile, thirteen with compressive, 
and one with zero, strains. 
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The effect plastic flow this plaster dome appears readjustment 
which the strains the upper surface under the loaded area are markedly 
reduced, with tendency approach more uniform values. Stations outside 
the loaded area show increase tension from day 0.85 day and then 
greater decrease, that the general net result decrease the original ten- 
sile values. However, stations the circumference including Station 3,4 have 
final stresses indicating that the restraint the near-by haunch causes bending. 

The foregoing statement covers the observed readings the upper surface. 
complete knowledge the effect plastic flow cannot obtained unless 
data regard the strains the under surface and the change deflection 
are also known. The the load bending well direct stress 
and the plastic readjustment complex. possible that the dome 
longer perfect spheroid elsewhere, and has been displaced horizontally 
well vertically. 

Quarter Segment load 534 also corresponded in- 
tensity per ft, 0.278 The sand load centered the 
meridian, and its boundaries were the 45° and 315° meridians. Compared 
with eighth-segment loading, the initial strains were much reduced. Using 
modulus elasticity 1.77 per in. and Poisson ratio 0.2, the 
computed stresses along the meridian are given Table 

the boundary the 315° meridian the circumferential strains the 
top surface were tensile and the meridional strains were compressive. Outside 
the loaded area the strains were all small magnitude, whereas their reversals 
suggested wave action the shell. 

This sand load remained place for days. Strain readings measure 
plastic flow were taken approximately and days. The readings 
varied considerably from day day with many reversals strain from com- 
pression tension and vice versa. Many readings were small the 
probable error reading. the center line the loaded area (0° meridian), 
the circumferential strains indicated that proximity the haunch affected 
Station and the hole the crown may have affected Stations and 21. 
The haunch apparently affected the meridional strains far out Station 
Neglecting the readings Station 1,2, the end 2.8 days the average 
circumferential strain was 3.8 and the individual values were fairly 
uniform; 4.8 days the average strain was 2.2 10-5, but the values 
Stations and departed considerably from this average. Similarly, the 
meridional strains averaged 3.5 2.8 days and 3.2 
4.8 days. 

the boundary the sand load the 315° meridian, the circumferential 
strain curve preserved its general shape time passed but steadily shifted 
the compression side, that initial strains that were wholly tensile 2.8 
days became wholly compressive. The 4.8-day readings were slightly lower. 
The meridional strains were wholly compressive for the initial readings (except 
for those Station 50) with the greatest strain recorded for all stations 
Station near the crown (11.95 There was marked shift 0.8 
days that all strains were tensile, except those Station near the haunch. 
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Thereafter the general pattern was preserved but the shift was toward small 
compressive strains. 

The general effect plastic flow the upper surface the shell along the 
center line meridian the loaded area was produce fairly uniform compres- 
sive strains which oscillate above and below between 1.8 days and 
4.8 days, except for the station close the haunch. These strains were usually 
greater than the initial strains. the boundary the sand load (315° 
meridian) the strains reversed, the circumferential strain becoming com- 
pressive, whereas the meridional strain became tensile 0.8 day, and then 
the strains slowly returned low compressive values. 

Outside the loaded area the initial strains had low values, some compres- 
sive, some tensile. From 0.8 day 2.8 days they all shifted toward greater 
tensile values, remaining fairly constant during this time interval. 4.8 
days the strains all tended approach zero strain: Seventeen with final com- 
pressive values (maximum 3.0 10-5), twenty four with final tensile 
values (maximum 2.95 and two with zero strain. 

Uniform Load Entire Dome tendency for the stresses 
decrease larger area loaded reasonable. the same sand load 
per assumed spread over the entire surface, the circumferential 
and meridional stresses can computed, assuming the shell 
the part the shell close the haunch, the circumferential stress equals 
per in.; and the meridional stress, per in. allowance was 
made for bending stresses caused the proximity the larger and more 
rigid haunch; but load four times that the quarter segment produced 
direct stresses the general magnitude those computed Table for 
the quarter segment. 


GENERAL SUMMARY 


The distributed load tests were made clarify the transition from 
eccentric concentrated load uniformly distributed load. They emphasize 
the fact that progressive concentration causes increasingly larger strains under 
the loaded area. 

area 6.25 in. The distributed loads were eighth seg- 
ment and quarter segment. The maximum stresses for the 
eighth-segment loading were comparable those obtained for the concentrated 
load the position producing the greatest recorded strains. The maximum 
stresses for quarter-segment loading were about sixth much. The part 
the dome well outside the loaded area had negligible stresses. 

The readjustments due plastic flow tended reduce the strains the 
upper surface the loaded area. the strains tended become 
tensile the first day quiescent loading. Thereafter the strains outside 
the loaded area tended approach zero. During these adjustments the re- 
straining effect the haunch was noticed near-by stations. 


Solution Masonry Domes,” David Coyle, Transactions, ASCE, Vol. 88, 1925, 
102. 
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APPENDIX MATHEMATICAL ANALYSIS 


Eric REISSNER 


facilitate the analysis the test data the model dome, theoretical 
equations were for the direct stresses and bending stresses com- 
ponents the total normal stress shallow spherical dome treated 
thin shell. The equations apply total load distributed uniformly over 
circular area with its center the apex the dome. Let “ber,” 
and used designate tabular functions. Also let 
and denote dimensionless variable and dimensionless constant, 
respectively (see Fig. 7). For small values 


and 


Equations for maximum stresses, are: Direct— 


and bending— 


TABLE FoR which subscript denotes 

the following data must known: The 

concentrated load the radius the 

spherical dome shell; the thickness 

the shell; the radius the circular 

loaded area, equal the actual loaded 

area; and Poisson’s ratio For con- 

crete 0.20 reasonable value. The arbitrary constants and 
Eqs. are determined solving the equations Table 
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The procedure consists first determining the values the arbitrary con- 
stant the length parameter being defined as: 


BENDING AND DIRECT MAXIMUM STRESSES 
STRESS; =0.5146 AND DEFLECTIONS 


Values of 


t2 


Maximum Direct Stresses £2 


1.6 
Values 


Values 


With this value enter Fig. 7(a) and read the maximum bending stress 


and the maximum direct stress oz. Fig. 7(a) not applicable for 0.2 
1.3; but Eqs. can solved directly interpolation Table 10. 
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Illustrative Example—A numerical example was solved and discussed 
detail the original analytical For this study the dimensions the 
model were used—namely: 102 in.; 0.25 in.; and 0.20. 


EQUATIONS THE FLAT SPHERICAL PROBLEM 


, , 

1.0000 0.01000 1.7331 —0.00250 0.5000 —24.615 1.1146 
0.9996 0.04000 1.0626 —0.7038 —0.01000 0.4999 —5.8802 0.7738 
0.9980 0.8998 0.6931 —0.6374 —0.02250 0.4997 —2.4276 0.5803 
0.9936 0.1599 0.4529 —0.5664 —0.3999 0.4989 —1.2342 0.4488 
0.9844 0.2496 0.2867 —0.4950 —0.06245 0.4974 —0.6946 0.3524 
0.9211 0.5576 0.05294 —0.3314 —0.1400 —0.1961 0.1971 
0.7517 0.9723 —0.04164 —0.2024 —0.2465 0.4585 —0.05330 0.1099 
~..| —0.2214 1.9376 —0.06703 | —0.05115 | —0.5233 0.2935 | +0.00716 0.03068 
—6.2301 0.1160 —0.01148 +0.01118 —0.7691 —0.8708 0.00344 —0.00016 


act three different 1.5625 in., 6.25 in., and in. 
The case which the loaded circular area equals 6.25 in. approaches closely 
the test case 6.25 in. composed four squares with outer boundaries 
in. each direction. For this case: 0.5146; 0.9050; and 
0.1728. The maximum combined stress the upper surface the shell 
equals: 


and 


These maximum stresses are equal for the meridional and circumferential di- 
rections. Fig. shows, not only the maximum stresses, but the variation 
stress over, and beyond, the loaded aréa the loaded area increased size. 

Table gives the values the functions used the equations Table 
obtain the constants and These stress equations are derived for 
concentrated load circular area the crown the dome. The fact that 
the loaded area not circular can neglected, provided the area small 
compared the total dome surface. The equations can also applied ap- 
proximately, the load not the crown, long the loaded area not 
near the springing, hole, other discontinuity. These results are approxi- 


mate because eccentrically applied gravity load not normal the surface 
loaded. 


APPENDIX II. NOTATION 


The following letter symbols, adopted for use the paper, and its dis- 
cussion, conform essentially American Standard Letter Symbols for Me- 
chanics (ASA—Z10.3—1942), prepared Committee the American 
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Standards Association, with Society participation, and approved the Asso- 
ciation 1942: 


arbitrary constant; formulas for the determination C2, and 
are listed Table 
subscript denoting 
diameter 
modulus elasticity 
thickness shell; 
length parameter defined Eq. 
subscript denoting 
concentrated point load, the resultant distributed load assumed 
acting the centroid pressure; 
radius, with subscript denote radius area concentrated loading, 
with subscript denote radius spherical shell; 
perpendicular distance from the axis symmetry spherical shell; 
strain, with subscripts and denote and 
tial,” respectively 
Poisson’s ratio; 
dimensionless quantity, with subscripts and denote 
constant” and “dimensionless variable,” respectively and 
unit stress, with subscripts and denote and 
stress, respectively, and subscripts and denote and 
respectively. 


The various tabular functions and are defined 
Eqs. These operators primed denote first differentials. 
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SUBSIDENCE THE TERMINAL ISLAND- 
LONG BEACH AREA, CALIFORNIA 


FREDERIC ASCE, AND EUGENE 


early 1945 the subsidence taking place and around Terminal Island, 
Long Beach, Calif., had reached maximum more than and the rate was 
approaching per yr. Therefore, protect existing structures and eval- 
uate the risks future engineering enterprises the area, became imperative 
determine the answers three questions: 

(1) What the cause the subsidence? 

(2) Can its future progress predicted? 

(3) What effect will have civil engineering works? 


This paper presents answers these questions. now appears certain 
that far the predominant influence the effect pumping oil from the Wil- 
mington oil field, which coincides location and extent with the subsiding 
area. Extensive data were gathered determine the factors producing oil 
field settlements, and laboratory tests were performed samples shales from 
different depths determine the relative compressibilities. Calculated theo- 
retical subsidence was compared with observed subsidence, correction factor 
was computed, and theory was then extended into prediction the future 
course the settlement. From these studies concluded that the maximum 
rate settlement has probably been reached and that henceforth the rate will 
gradually diminish until eventually, perhaps ten fifteen years, the max- 
imum subsidence will reach about ft. 

Elevations waterfront structures, gravity flow pipe lines, etc., may have 
adjusted accordingly. Since differentials will gradual, only very long 
rigid structures will subjected warping cracking. 


comments are immediate publication; insure publication the last dis- 
cussion should submitted March 


Rear Admiral, CEC, Engr., New York, 
Foundation Engr., Frederic Harris, Inc., New York, 
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NATURE OBSERVED SUBSIDENCE 


Occurrences differential movements the Los Angeles-Long Beach har- 
bor areas California have long been recognized However, 
comparison records tide gages maintained several points the Los 
Angeles and Long Beach Harbor areas for many years with observations 


Pacific Coast 


Benchmarks 
Tidal (Datum Used) 


Jolla and San Francisco show that the gage berth 60, municipal pier No. 
Los Angeles harbor, referred has not changed elevation since 
1931 more than about 0.1 ft. 

Subsidence various bench marks the affected area during different 
periods was determined, general, subtraction corresponding elevations. 
Those for 1933 were taken representing zero subsidence, because study the 
Long Beach outer harbor tide gage record, well frequent levelings the 
Southern California Edison Company, reveals the 1933 leveling the most 

extensive survey prior the beginning the rate subsidence pre- 
vailing 1945 and 1946. 


Numerals parentheses .thus: (1), refer corresponding items the Bibliography (see Appendix). 
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Fig. shows contours the subsidence that occurred from 1933 1945. 
Profiles taken several dates along representative line through the area are 
given Fig. shows the increase subsidence certain bench marks 
plotted against time. 


GEOLOGY THE REGION 


The geologic history the region has been amply described previous re- 
ports and papers (5)(6). Generally, all the sediments the upper 7,000 are 
poorly consolidated and 


Benchmarks 


The Miocene and lower Pliocene 
strata are warped into gentle 
anticline with its axis 
northwest-southeast 
forming dome 
whose highest parts are close 
the east end Terminal Island. 
These strata are all highly 
faulted, and form the structure 
the Wilmington oil field, 
present one the largest pro- 
ducing fields California. The 
faults have been traced upward 
unconformity the base 
the Upper Pliocene, Pico, 
formation, about 2,000 below 
the surface, but not above this 
point Varying thicknesses 
beds the underlying Repetto formation indicate that the motion along these 
faults was still occurring until late Pliocene times, but did not continue into 
later periods. The upper 1,800 sediments generally contain more less 
fresh water the voids, and have been the source fresh water supply 
residents the basin for many years. 


Subsidence, Feet 


March 1946 


Factors THE SUBSIDENCE 


least five major sources settlement operating this area have con- 
tributed the observed subsidence varying degree. They are: (a) Changes 
the land surface from filling and dredging; (b) pumping water for indus- 
trial and domestic uses; (c) pumping water during dewatering the site 
dry dock No. United States naval drydocks; (d) pumping oil from five 
major oil zones the Wilmington oil field; and (e) tectonic forces. Any study 
this problem resolves itself finally into opinion based factual obser- 
vation and accepted principles statics and mechanics, concerning the relative 
importance these possible sources. Therefore, each them discussed. 

Filling and Dredging the Los Angeles and 
Long Beach harbor areas has been progress since about 1900. that time, 
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Terminal Island, then called Rattlesnake Island, was low sand spit. has 
been gradually enlarged dredged fill several times its initial size, and was 
eventually built about El. +13. However, far known, fill 
materials have been brought man into this area from any distance. Con- 


Years 
1931 1933 1934 1935 1936 1937 1938 1939 1940 1941 1942 1943 19441945 1946 1947 


BENCHMARKS 
(a) Tidal Datum 

(Los Angeles Harbor Datum) 
H-97 
(c) Ford Avenue Bridge 
(d) Pier Tidal Gage (Long Beach) 
(e) 215 


Jan. 1945 
July 1945 


Jan. 1947 


Subsidence, Feet 


sequently, part the over-all subsidence 
the area can attributed increase 
the total weight overburden, except for the 
accretion that has taken place from natural 
causes. 

Theoretical calculations determine the 
probable amount subsidence resulting from 
the changes surface loadings the vicinity 
Terminal Island involved three steps: First, 

surface deposits over different parts the area, 
and the time sequence these changes; second, the total consolidation un- 
derlying compressible strata that should ultimately occur; and, third, the in- 
fluence time effects during the consolidation, which reflect the proportion 
the settlement that would cause observed movement surface bench marks. 
From study many consolidation tests clays and loams several investi- 
gators (7)(8), was concluded that, the total subsidence measured, perhaps 
the strip along the seaward side Terminal Island, about 0.2 could 
attributed surface loading. 

General Effect Pumping the ground-water table lowered, all 
the sediments between the initial and final ground-water surfaces will experience 
loss buoyancy, and will rest more heavily the underlying soil skeleton, 
which will eventually compress. The compression may not occur suddenly, 
since reduction volume requires the escape excess water the voids. Ac- 
cording the theory consolidation advanced Karl Terzaghi, ASCE, 
the time required varies the square the thickness the stratum under- 
going compression. Since the more compressible strata are generally the fine- 
grained silts and clays, which may exist layers between more pervious sand 
beds, time needed for such volume change before static equilibrium re- 
established. 

two more pervious beds are separated continuous impervious clay 
strata, the water pressures the voids the lower bed may not conform the 
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ground-water table the upper bed (9). Also, the water pressure 
lower confined aquifer reduced, the pressure and level the water the 
surface aquifer will not affected (see Fig. 4). However, the lower aquifer 
will undergo compression throughout its thickness. 


Pressure Level for Ground Surface 
Water Lower Aquifer 


Drawdown Pressure Level 
Lowered Pressure Level 


Ground Water Level 
Pervious 


Ground Water 
Level (Unchanged) 


Total Pressure Line 
Pervious Lower Aquifer 


Effective Pressure Line 


Impervious 


Pressure, Tons per Pressure, Tons per 


(a) Before Pumping 
(6) During Pumping from Lower Aquifer 


From the foregoing, seen that there are two distinct effects that may 
result from the pumping liquid from the ground: 


Loss buoyancy the soil particles. This occurs the ground-water 
surface lowered, causing compression all underlying strata. 

Increased intergranular pressures pervious zone confined between im- 
pervious layers. This occurs the liquid pressure the confined 
zone reduced, and results compression the confined zone and 
adjacent parts the (impervious layers) top and 
bottom. 


The same general principles apply the effect pumping any liquid— 
whether water oil. 

Water Pumping for Industrial and Domestic water the Long 
Beach-Wilmington coastal plain obtained from two water-bearing zones. 
The upper Gaspar zone occurs the surface deposits recent age and con- 
sists from 100 200 alluvial deposits sand and gravel interbedded 
with silt and clay. includes both surface waters and confined waters. The 
lower Silverado zone confined aquifer Pleistocene deposits sand and 
gravel with intervening layers silt and clay. The zone varies from 150 
500 thickness. varies depth, being from 200 900 below the 
land surface the Wilmington section. 

The Newport-Inglewood fault forms partial barrier the flow under- 
ground water. The area the northeast the fault was originally highly 
artesian, deriving pressure from water entering the pervious strata the foot- 
hills the mountains which surround the basin. The area the southwest 
the fault was probably only slightly artesian originally. Water-bearing zones 
this area are replenished mainly underflow across the Newport-Inglewood 
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fault, through imperfections the fault barrier. result, the effect 


lowering the ground-water levels pumping much more noticeable the 
coastal area south the fault than the areas the north. 

Development ground water the Los Angeles coastal plain began about 
1870 (10). Records show that wells drawing from the Gaspar and Silverado 
zones had static water levels from above sea level 1903, gradually 
dropping about sea level the Gaspar zone and much from 
below sea level the Silverado zone (11). 

Changes elevation the ground surface caused pumping water from 
the lower strata are almost certainly part the reason for early survey diffi- 
culties. The influence probably still contributing, but moderate extent, 
with the fraction total subsidence resulting from this source amounting 
inches, not feet. The center the depression would expected lie some 
distance north Terminal Island. 

Dewatering the Site Dry Dock No. 1—Pumping water from surface and 
confined pervious strata during the construction dry dock No. drew down the 
pressure levels these beds from sea level about El. —75, based mean 
lower low water datum. The pumping covered period from June, 1941, 
April, 1942. Maximum drawdown water levels did not occur until some 
weeks after pumping was begun. After pumping was stopped, readings 
near-by wells indicated that the pressure levels rapidly rose their initial ele- 
vations, close mean sea level. 

Water levels wells the upper aquifers were affected varying amounts 
far north miles. subsidence was observed the affected area, 
varying with the drawdown pressure level, roughly inverse proportion 
the distance from the dock site. The Southern California Edison Company, 
whose Long Beach steam station about 2,000 northeast the dock, in- 
stalled fifteen observation wells, some deep and some 150 deep, meas- 
ure the reduction water pressure level the two principal pervious strata 
the local part the aquifer (7). were also made the motion 
bench marks around the plant. 

Within radius 2,000 the dock site the ground settled much 0.45 
around the casings many oil wells. most instances, the construction 
the derrick foundations permitted the ground and concrete base slabs settle 
around the casings without appreciable damage either casings concrete. 
There were few occurrences cracked slabs, and damage was reported 
least one well casing. Differential movement between the ground surface 
and the casings was observed only during the period pumping the dock site. 
Prior and since that time, the oil well casings have apparently behaved 
they were part the ground surface. 

estimate how much the total subsidence was caused dewatering the 
dock site, the theoretical consolidation resulting from the effect the observed 
drawdown liquid pressure levels was computed. 

From study the stratification and drawdown records, appears that the 
substrata affected the pumping consist sands and gravels the Gaspar 
aquifer, extending depth from about El. —70 about El. —270. The silty 
clay materials bounding, and within, the aquifer were assumed have 
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cient compressibility 0.144. Calculated contours equal subsidence 
resulting from this pumping appear consistent with the measurements motion 
around well casings. When pumping was stopped, the spring 1942, water 
levels rose their previous values, and the consolidation ceased. fact, 
slight rebound, increase volume the consolidated strata, occurred. 


‘Therefore, this settlement considered represent merely incident the 


general subsidence phenomenon the area, and its effect was felt solely the 
period from June, 1941, April, 1942. 

Pumping from the Wilmington Oil mechanics liquid flow 
during production from oil well are generally the same those for pumping 
water from confined aquifer. The oil sands forming the equivalent the 
aquifer are bounded top and bottom impervious layers silt 
and clay compacted the weight the overburden. The oil the sand voids 
initially subjected pressures which depend mainly the depth the 
strata below sea level. oil withdrawn from the well, there pressure 
reduction, which produces gradient toward the well from all sides. This 
gradient depends the rate pumping, the permeability the sand, and the 
extent and uniformity the stratum affects the resistance liquid re- 
placement the oil withdrawn. oil wells the flow complicated somewhat 
the presence large amounts gas, either free dissolved the oil, and 
the presence least two liquids, oil and water, with different physical proper- 
ties (12). 

with pumping water from confined aquifer, reduction the liquid 
pressure the oil sands produces compensating increase the intergranular 
effective pressures the soil structure within the stratum affected. Com- 
pression the oil zones results. The amount compression influenced 
the amount decrease liquid pressure and compressibility the mate- 
rials within and bordering the affected zone, top and 

Compression oil zones has frequently resulted subsidence the ground 
surface above ((5)(13)(14)(15) and unpublished reports). significant that, 
the Los Angeles basin, noticeable subsidences have been recorded first 
order leveling the areas number oil fields. Among these are the Wil- 
mington, Signal Hill, Santa Springs, Playa del Rey, Inglewood, and Hunt- 
ington Beach fields. However, first order levels have not been run across many 
fields sufficiently often detect possible The relation between the 
amount the unchecked compression the oil zones and the corresponding 
surface subsidence depends the structural strength the sediments be- 


tween. quite common for extensive oil fields experience 


surface subsidence. the other hand, subsidence large areas from any 
cause seldom noticed inland areas. Along the seashore, the relatively 
constant sea level provides automatic reminder whenever subsidence ap- 
preciable amount occurs. 

The Wilmington oil field was discovered 1932, but little attention was 
paid until the last half 1936, when carefully cored and tested well re- 
vealed the presence four major oil zones. The drilling boom started the 
Los Angeles harbor area thereafter, spreading the Long Beach harbor area 
1937 (6)(16)(17). Maximum production 100,000 bbl oil per day was 
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reached February, 1938, dropping gradually about 82,000 bbl per day 
1940, and then rising again around 100,000 bbl per day 1944. The field’s 
oil production rate 1945-1946 the highest California wide margin. 
This time sequence agrees generally with the occurrence most sub- 
sidence, apparent from comparison between Fig. and Fig. The 
anticlinal structure the oil field shows remarkable correlation with the con- 
tours equal settlement. 

ascertain the importance oil pumpage from the Wilmington oil field 
producing settlement, quantitative studies were made 1945 based actual 
measurements the several factors governing consolidation sediments the 
oil zones. These factors consist primarily the stratification the sands and 
shales, the compressibilities the sands and shales, and the reduction liquid 
pressure the oil and water the voids after different periods time. 

Stratification was determined from oil well electric logs and coring records. 
Fairly detailed information available and has been plotted contour maps 
and profiles the oil companies and various other agencies interested the 
field. The general structural characteristics are shown Fig. There are 
five major zones: Tar, Ranger, Upper Terminal, Lower Terminal, and Ford, 
that order, varying depth from about 2,200 the top the Tar zone 
about 4,700 the top the Ford zone. The basement complex igneous 
rock, schist, encountered about 6,500 below mean Each zone 
varies thickness some extent; but except for faulting the differences 


Mean Sea Level Gaspar Aquifer 


Upper Terminal 


Elevation, Feet 


Zone 


(Puente) 


Middle Miocene (Monterey Conglomerate) 


THROUGH WILMINGTON FIELD 


stratification between one point and another are remarkably small. The fault- 
ing, however, represents the most serious complication. series faults ex- 
tending roughly north-south direction divide (see Fig. the field into five 
major fault blocks, which are themselves broken innumerable smaller faults, 
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that most cross sections the field not attempt trace all details the 
complex pattern. The locations faults are determined mostly inductive 
reasoning, through study the oil well electric logs (18). 

Fluid pressures oil wells constitute one the most important subjects for 
study oil producing companies. Measurements flowing and shutoff pres- 
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sures are made determine the productivity the well, which enables esti- 
mates made future production, and also data obtained for use 
operation insure maximum total production from the field. The usual cus- 
tom shut off the well for hours longer before taking the measurements, 
allow time for the pressure build close the maximum. pressure 
recording instrument lowered into the well, and readings are taken intervals 
that pressure-depth curve can plotted. 

Initial pressures the field, determined from large number initial 
measurements early wells, were approximately the same the under- 
ground liquids were composed salt water with the surface mean sea level. 
Continued pumping has generally caused steady reduction the pressures. 
pressure any depth reduced terms pressure level, the level 
which column salt water would stand static equilibrium tube whose 
lower open end inserted one the confined sands oil zone, then may 
stated that the pumping has caused gradual reduction the pressure level 
below the initial mean sea level elevation. Typical curves showing pressure 
levels various times for different wells are.shown Fig. Pressure levels 
the various zones are not the same, because the impervious shales separat- 
ing them. has also been observed that the pressure levels opposite sides 
fault are likely differ, the displacement along the fault plane apparently 
having caused breakdown the normal sand structure and the formation 
relatively impervious materials, which act barrier the flow liquids. 

Pressure relations the oil zones are complicated some extent the 
presence gas. Initially, gas the Wilmington field believed have ex- 
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isted near saturated solution, without any appreciable free gas cap. How- 
ever, pumping continued, pressure reduction has caused the formation some 
free gas, which may have accumulated gas caps the crest the structure 
and along faults some zones. The near saturated condition the original 
oil indicated the fact that wells originally gushers, the gases 
rapidly expanding out the solution the casings, carrying the oil upward 
the typical gas lift action. Gas production from the field still fairly high. 
The free gas affects intergranular pressures generally causing smaller pres- 
sure increase than would have occurred had the gas been absent, since the 
liquid pressure reduced. Development gas cap increases the unit weight 
the sand within the depth the cap, because the change from submerged 
dry weight; but, since the gas confined beneath impervious dome 
shale, the buoyancy the liquid volume displaced gas very small unit 
weight produces lower effective pressures the strata below the shale. Unless 
the gas cap quite thick, its presence has relatively small effect inter- 
granular pressures. Since the scanty information available for the Wilmington 
field indicated large volume free gas, its effect consolidation was dis- 
regarded these subsidence studies. 

great many pressure measurements were made the first five years 

the Wilmington oil field, particularly the large property 
owners. There are many smaller operators, particularly the northwest part 
the field, known the Area,” who have made practically pres- 
sure measurements that section, although there has been considerable pump- 
ing since early the field’s history. Also, pressure measurements the large 
operators were less frequent after the beginning World War II. Since in- 
sufficient measurements were available for 1945, data were also included for 
measurements made 1943 and 1944, which were given the same weight the 
1945 measurements, preparing pressure contour maps for areas and zones 
where the later data were not available. Thus, discrepancies were introduced 
between calculated and observed subsidence for 1945, will noted hereafter. 

Compressibility sedimentary materials the pressures existing oil 
zones several thousand feet below the surface matter about which little in- 
formation available (19)(20)(21). obtain data more applicable this 
problem, was decided perform tests actual samples materials taken 
from oil wells the Wilmington field. For several reasons, shale cores only 
were selected for testing. the first place, shales, the product silts and 
clays, composed fine-grained particles, were expected much more com- 
pressible than sands. Previous tests had demonstrated fairly conclusively that 
sands, and particularly sandstones, show extremely small volume change when 
subjected very high pressures. Secondly, relatively undisturbed shale cores, 
suitable for testing, could obtained with much less difficulty than sands, be- 
cause their cohesive quality. Most the sands from the Wilmington field 
are quite friable, and break easily. 

was also felt that, rather than adopt triaxial tests, the procedure would 
more nearly duplicate the consolidation process nature the applied load 
were confined the vertical direction, the same time preventing lateral 
volume change, essentially the usual practice for standard consolidation 
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tests silts and clays. most consolidation phenomena, where the load 


applied relatively wide area proportion the thickness the compress- 
ing stratum, the volume change predominantly vertical. This condition can 
approximated the laboratory subjecting sample increased vertical 
pressure, the same time restraining against lateral expansion. 

Available samples consisted cores drilled from various oil wells. Most 
these cores were several years old, having been obtained during the period when 
the oil companies were exploring the field. few instances, fresh cores were 
obtained from wells where drilling was currently progress. The cores varied 
size from in. in. diameter. 

Equipment used performing the tests was designed accommodate two 
sets samples simultaneously. The apparatus, which pictured both as- 
sembly and detail Fig. employs manually-operated screw-driven jack 
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for applying the load, and beam gage for measuring the load. Porous 
stones, placed both ends each 1-in.-thick sample, enabled application to, 
and drainage water from, the samples. Two sets three equally spaced 
1/10,000-in. micrometer dials were used measure the compression the 
samples. determining the reduction thickness each sample from the 
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average the dial readings, was necessary measure the deflections occurring 
the porous stones and make appropriate adjustments. 

Samples tested were first placed saturated atmosphere reestab- 
lish them, nearly practicable, saturated condition. They were then 
accurately trimmed lathe—all weathered damaged parts being removed 
the process—to fit the internal diameters the sample holders, and the ends 
were faced provide samples the desired thickness. After machining, the 
samples were kept saturated atmosphere until ready placed the 
testing apparatus. 

Immediately after application small initial load, screws holding the 
sample holders position during the setting-up process were released, and 
water was applied the top and bottom the samples through the porous 
stones, which were kept constantly wet throughout the remainder the test. 
Load was applied series increments, with each increment being twice its 
predecessor, except that two more additional increments were superimposed 
the basic series increments the vicinity the preconsolidation load. 
removal the samples from the testing apparatus, the moisture content, dry 
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density, and the specific gravity the soil grains were determined for each 
sample. typical curve pressure versus voids ratio data obtained from the 
shale tests shown Fig. for 2.416-in. diameter sample. 
Typical compression-time curves are shown Fig. point Fig. 
represents the average set three dial readings, expressed terms the 
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percentage the total deflection observed the end the first 1,440 min dur- 
ing which particular load was applied. 

For the highest pressure which two the samples were subjected, was 
observed that some the material was squeezed plastic flow into the space 


N w 


Time, in Minutes 


between the porous stones and the sample ring. these instances, arriving 
the straight-line approximation, correction was made the voids ratio for 
the last point the curve proportion the amount material displaced. 
This plastic flow probably would have been avoided the diameter the 
porous stones could have been made closer that the cylindrical openings 
the testing rings. other evidences plastic flow were observed. 

For comparative purposes, the compressibility each sample can defined 
terms coefficient determined the slope the assumed initial com- 
pression curve, the change voids ratio for tenfold increase pressure. 
The initial compression curves were assumed tangent the steepest sec- 
tions the test curves, and represent approximation the probable com- 
pression curves the shales their original condition the ground. these 
compression coefficients are plotted against depth, definite variation seen 
exist Theinference that the shales encountered near the contact be- 
tween the Lower Pliocene and the Upper Miocene deposits, represented the 
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Tar, Ranger, and Upper Terminal zones, are more compressible than are mate- 
rials greater lesser depths, under the same conditions loading. 

compute theoretical consolidation accurately, would necessary 
know minute detail the variation compressibility all strata encountered 
vertical section through the oil 
zones. Suchinformation cannot feasibly 
obtained. certain that the com- 
pressibility varies rapidly stratum 
stratum, and even foot foot, ranging 
from the high values the more uni- 
formly fine-grained shales, the rela- 
tive incompressibility the cemented 
sands, which probably have properties 
close those concrete. However, 
even with the most conservative assump- 
tions, clear that compression the 
oil zone sediments can easily account 
all the observed subsidence. the 

total thickness each zone assumed 
have the same compressibility the 
tested shale samples from the zone, then 
rough computations show reduction 
thickness the oil zones amounting 
more than ten times the observed settle- 
ment, based the pressure data 1945, 
for point near the center the subsiding 
area. Data similarly derived for other 
points indicate that the distribution the 
consolidation corresponds closely the 
distribution the observed settlement. 
the subsidence was made, based the 
following procedure: Ultimate drawdown pressure levels each oil zone was 
estimated, and correction factor applied the total thickness the oil 
zones was computed. This factor was determined deducting from the ob- 
served subsidence point near the center the depression the computed 
settlement from water and dock site pumping, 1945. The difference was 
assumed represent the effect oil pumping that time. constant per- 
centage the thickness each zone was chosen make the calculated 
consolidation that particular point equal this difference. Using this 
factor, the same for all zones over the entire area the field, together with the 
average compressibility coefficients, different for each zone, and the total thick- 
ness the zones, which varied over the area the field, consolidation was com- 
puted some two hundred points. These computations were made using the 
pressure data 1945 and also using estimated pressures that may occur the 
future. Theoretical subsidence from oil pumping, 1945, shown the con- 
tours Fig. 11. 
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The many factors that influence the theoretical subsidence computed this 
manner perhaps warrant some discussion. the first place, the consolidation 
computed for the oil zones was assumed produce equal surface subsidence, 
except for smoothing out local irregularities. This assumption neglects the 
effect possible lateral stress changes, which probably may result from the 
structural strength the deep beds. The domed structure these beds itself 
may result certain amount arching action, which tends transfer laterally 
the vertical loads applied near the center the field, relieving the centrally 
located shales some the increased pressure they would otherwise ex- 
perience. This action may reduce not only the local irregularities, caused 
local variation drawdown the liquid level, but also the total surface sub- 
sidence the center the depression, spreading over wider area. How- 
ever, the same errors from this source exist the calculated 1945 consolidation 
the calculated ultimate consolidation. Since the predictions for the future 
are derived from the difference, the net error should relatively small. This 
factor may partly explain the lack agreement detail between the contours 
for calculated oil zone consolidation (Fig. 11) and the contours observed sub- 


Los Angeles River 


sidence (Fig. 1), neglecting the relatively small and uniform consolidation from 
water pumping. 

The estimates future drawdown each the oil zones probably con- 


stitute the greatest single source error the estimate future subsidence. 
They were obtained purely arbitrary manner, but were based some ex- 
tent past experience and certain limiting factors, 
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Leading petroleum engineers and geologists familiar with the Wilmington 
field agree that eventually, result pumping, oil zone pressures will 
greatly reduced. some oil fields, oil withdrawn, edge waters encroach- 
ing from the sides the structure replace freely enough prevent extreme 
pressure reductions the oil. such cases, the oil reservoirs may com- 
pletely depleted with only moderate reduction the oil pressure level—the rate 
oil removal governing large degree the amount pressure reduction. 
the Wilmington field, however, appears that reasonable pumping rates 
cannot maintained without gradual reduction the liquid pressures, since 
experience shows that edge-water replacement not sufficiently rapid. This 
may caused low permeabilities the oil bearing strata, fault barriers, 
confinement these strata within relatively impervious layers. 

With these factors mind, was concluded that ultimate liquid pressures 
may very likely reduced about 200 per in.—equivalent about 400 
water head the top each zone. The distribution this drawdown 
over the area the field uncertain matter, but appeared that the report 
prepared 1941 the Wilmington Field Engineering Committee contained 
sufficient data for reasonable approximations. Maps prepared this com- 
mittee showed outline the proved limits each zone the field, and 
assumed general that these limits enclose the area where the most severe 
pressure reduction will occur. 

The manner tapering off the assumed drawdown the edges the field 
was determined judgment alone. The nature the gradient from the large 
drawdown within the field sea level some distance from not known. 
Pressure measurements within the field show that fairly steep gradients can 
exist, particularly the less pervious sediments the Ford zone. Across 
faults there also quite commonly very steep gradient. the edges the 
field, there are very few measurements pressure, yet seems certain that the 
pressure gradients must extend for considerable distance into the surrounding 
edge waters, otherwise inward flow edge water would occur, 

The prediction, which involves the difference between estimated ultimate 
and measured 1945 pressure levels, also suffers from certain inaccuracies the 
values, which include some data actually obtained 1943 and 1944. 

should mentioned that the observed 1945 subsidence based the 
measurement relatively few points, scattered rather irregular fashion over 
area. There are number regions where observations all were 
made. For this reason, quite possible that the actual subsidence may not 
have been quite regular outline the smooth contours Fig. would 
imply. 

Compressibility the oil zone sediments should not change appreciably 
with time, the relatively brief geologic interval this subsidence. There- 
fore, the true average compressibility were known for all zones given 
point, and the increase drawdown after 1945 were the same for all zones, fu- 
ture consolidation could predicted quite accurately. However, pressure re- 
ductions expected the future vary the different zones. For example, 
fault block III, near the center the depression, much larger proportion 
the ultimate pressure reduction has already taken place the Upper Ter- 
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minal and Lower Terminal zones than has occurred the Tar, Ranger, and 
Ford zones. Thus, the proportion the total consolidation caused each 
the zones, for given increment pressure reduction, becomes important. 
This proportion depends the total thickness compressing materials each 
zone and the average compressibilities those materials. 

Thicknesses each the zones were determined principally from the con- 
tour maps prepared the Wilmington Field Engineering Committee. These 
maps showed the variation elevation the top each zone. The thick- 
nesses the Tar, Ranger, and Upper Terminal zones were taken the differ- 
ences between the elevation the top the zone and the elevation the top 
the next zone below. The Lower Terminal zone has not generally been pro- 
ductive deep the top the Ford zone, and has not been pumped for this 
full thickness. Lower Terminal zone was therefore taken have av- 
erage thickness 500 ft, based cross sections prepared the committee and 
others and various electric logs. Similarly, the thickness the Ford zone 
was taken 750 ft. The proportion each zone composed shale layers 
varies over considerable range, but averages about 50%. 

Calculated consolidation, both for 1945 and for the future, was obtained for 
each zone multiplying the full thickness the zone correction factor 
adjust the calculated change thickness, for point the center the 
depression, with the observed subsidence 1945. low ratio 
between observed and calculated unadjusted settlement may explained 
the facts that: (a) Total thickness oil zones was used the calculation; 
and each zone was assumed have the compressibility the average the 
tested shale samples. Thus, the presence sands low compressibility dis- 
regarded are the relatively incompressible sediments composed shales 
mixed with sands, which affect the average greatly. The presence many 
water-bearing strata that may not greatly affected pumping also neg- 
lected, since well casings are not perforated water zones. Also, previously 
mentioned, the tested samples were chosen for their suitability for machining, 
were very fine-grained, and therefore were more compressible than the major- 
ity shales examined the oil companies’ storage racks. 

This method leaves much desired, when compared with the degree 
accuracy usually sought consolidation studies. However, believed 
sufficiently accurate for the purpose required, because most errors inherent 
the theoretical approach are canceled the adjustment already sizable 
consolidation, representing major part the possible total, which has been 
actually measured. 

all the calculations oil zone consolidation, was assumed that the time 
lag between the imposition load, through reduction the pressure level, and 
the corresponding decrease voids ratio the compressible shales has been 
relatively small. Inspection the time-compression curves obtained the 
laboratory tests (Fig. shows that majority the curves not have the 
same shape typical time curves explained the Terzaghi theory consol- 
idation (22). 

general, the shape the time curves resembles the results certain soil 
tests which, because the character and thickness the sample, the so- 
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called time effects have practically concealed the ef- 
fects, those resulting from squeezing out water. such cases, instead 
the typical reversed S-shape the semilogarithmic plot, the time curve may 
straight, even slightly concave upward. 

The time required for primary consolidation assumed vary the square 
the thickness the compressing layer. The time required for secondary 
consolidation more less unknown factor, but generally held 
independent the thickness the layer. other words, secondary consoli- 
dation may occur nature, layer many feet thick, the same rate 
occurs the thin laboratory sample. 

Those test curves that display some similarity the theoretical time 
curve indicate that for pertinent load increments 90% consolidation takes place 
less than 400 min. this basis, the average thickness individual shale 
laminations, between sand streaks, assumed ft, then 90% consolidation 
the oil field would require about 290 days. 

For the compression oil shales considered this study, many unknown 
factors influence the time rate. Even the average thickness the compress- 
ible layers evaluate, since there are many layers varying 
thickness each the oil zones. Consequently, attempt was made 
calculate how many the time effects result from primary compression, inas- 
much the laboratory curves show this influence minor compared 
that the secondary plastic characteristics. The time required for these lat- 
ter effects matter days, instead months years, and was neglected 
computing theoretical consolidation any particular time. Since the loads 
themselves are applied gradually, over period years, the error introduced 
this procedure probably not very large. 

Influence Tectonic well known that, the coastal area 
California, motions the earth’s crust are constantly progress. attempt 
will made describe the many faults the subsiding area and around it, 
describe the seismic history the region (5)(13). obvious that 
area interlaced with active faults, where minor quakes are 
continuously the near-by Pasadena Seismological Laboratory, movements 
along fault zones are more less likely occur the result slow relative 
motion adjacent land masses. 

far known, there reliable method estimating timing the 
probable occurrence movements resulting from tectonic forces. However, 
far evidence derived from observations oil wells concerned, appears 
likely that movement along the many faults the oil field area ceased Plio- 
cene times. During these movements, the oil field was generally lifted upward, 
with downward “‘cave-ins” some the central blocks the field. general 
upward movement this section the coast evident also from the terraced 
contours the San Pedro hills around Point Fermin. this grad- 
ual uplift may still progress. 

The area that might affected tectonic forces likewise indeterminate. 
The Newport-Inglewood fault, which the major known active fault zone pass- 
ing through the Los Angeles coastal basin, lies some miles northeast the 
center the subsiding area. Although might appear that the subsiding 
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area, which coincides with the Wilmington oil field, very highly faulted one, 
compared with other areas, should remembered that knowledge the sub- 
surface details this area highly developed from observations made many 
hundreds oil wells. There assurance that many other similarly faulted 
areas the Pliocene and Miocene sediments not lie undetected merely for 
lack exploration. 

weighing the probability tectonic forces causing major part the 
subsidence, several factors must considered along with the known possibility 
such movements this region. One such factor the lack any marked 
increase seismic activity recent years the Long Beach area. Since the 
rate subsidence there has increased degree unprecedented California 
history, tectonic forces and movements underlying rock masses were largely 
responsible, the area might properly expected develop nest shocks, 
the earth stresses produced during such severe strains readjusted themselves. 
However, although innumerable shocks have occurred Southern California, 
such nest shocks has been observed Long Beach. 

Another such factor the absence damage oil well casings the Wil- 
mington field. there had been any appreciable motion along the faults trav- 
ersed the wells, pinching and shearing off the casings puncturing 
cement plugs should have taken place; and seems highly improbable that all 
the relative motion along these inclined faults could have occurred depths 
below the oil wells, without producing considerable displacements the upper 
sediments. this connection, experiences Dominguez Hills 1941 indicate 
that motions along faults little fraction inch were sufficient 
collapse oil well casings passing through the faults (23). many cases, mo- 
tions more than foot were observed, shearing off casings completely—yet 
this condition occurred oil field where surface subsidence, took place 
all, was small unnoticed. considered, therefore, that, al- 
though slow tectonic movements have occurred the area question, and 
probably are occurring present, these movements are very gradual, possibly 
the order 0.04 per less, vertically. They may either upward 
downward; and available facts are insufficient determine which, estab- 
lish the rate, concealed beneath the very much more rapid subsidence 
phenomenon centered over the oil field. 


Estimated ultimate subsidence, after 1933, from oil and water pumping 
shown Fig. 12, which some the irregularities contours derived 
mathematical basis have been smoothed out, the assumption that arching 
and other effects will nature reduce abrupt changes and transitions the 
over-all pattern. This contour map, means considered infallible, can 
serve guide for future construction this area. The most vital conclusion 
brought out the contours that all probability there will large-scale 
inundation the easterly end Terminal Island. Wherever present free- 
boards above extreme high tide level exceed ft, the chances are that there will 
overtopping the ground surface. 

Nevertheless the possibility differential movement should considered 
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the design any large structures built the subsiding area. gen- 
eral, since the depression the ground surface covers very wide area, with 
flat side slopes, the amount differential motion anticipated the length 
most buildings will very small. the buildings are sufficiently flexible 


withstand such motion without excessive stresses, particular damage will 
result. 

Continued subsidence will necessitate restudy the elevations and grades 
certain facilities. The design gravity sewers should re-examined the 
light probable future differential changes grade. Lack proper drainage 
may result some places the subsidence continues. The same problem ap- 
plies surface grades. Such sources future difficulties can thus antici- 
pated some extent, and provisions can made insure continued satis- 
factory operation. 

conclusion, this study indicates that this vicinity should continue serve 
the site for industrial development. Present structures will continue func- 
tion satisfactorily suitable preventive measures are adopted avoid known 
possibilities trouble; future structures can designed allow for probable 
future subsidence, which will continue for many years slowly reducing rate. 
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TANEOUS LINEAR EQUATIONS 
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Simultaneous equations are great importance modern structural engi- 
neering, and the purpose going through solution is, obviously, obtain 
numerical results for specific set conditions. often equal impor- 
tance, however, know the effect change the geometric and elastic 
properties the structure. Without such knowledge impossible study 
the design with view toward obtaining the most economy and certain 
other desirable qualities. Such studies, except for the simplest structures, 
can made only the expense re-solving the equations each time change 
made the coefficients, and the labor involved great that studies cannot 
ordinarily made except very roughly and sketchily. With the theory pre- 
sented this paper, the effect the solution produced change any 
number the coefficients any one equation, any number coefficients 
any one unknown, can computed quickly and simply one step. 

Although the variation coefficients the primary purpose the paper, 
other useful properties the theory, such the solution equations and the 
rapid correction errors, follow simply corollaries. the actual operation 
the method, however, the procedure the same whether coefficients are 
being varied set equations that have already been solved, whether 
the purpose solve the equations. is, matter fact, possible 
vary the coefficients endlessly and any particular therefore becomes 
merely stopping point endless series variations—from which 
follows that the method very flexible. 

Variation coefficients method solving simultaneous equations has 
several marked advantages over the classical methods and leads also the 
derivation several useful properties such equations. The numerical 


comments are invited for immediate publication; insure publication the last dis- 
should submitted March 1948. 
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work required for solution the method presented this paper very much 
less than that required standard methods; and, what perhaps more 
important, fewer significant digits are necessary the intermediate computa- 
tions obtain any desired degree accuracy the final solution. fact, 
with little care and using few fairly obvious arithmetical manipulations, 
results sufficiently accurate for design purposes often may obtained 
slide rule for the bulk the computations. The method exact one, how- 
ever, and computations may carried out any degree accuracy but with 
much less labor than that required the classic methods. 

addition, the method solution has another very important property. 
will shown, the solution falls into distinct steps, partial solutions, 
each and every one which complete itself. result, the computation 
may checked completely and absolutely the conclusion each step 
variation, that not necessary wait for the end the solution check 
for arithmetical error. Should the computer, however, either neglect 
check his work intermediate points, discover error the original 
setup the equations, the basic theory the solution such that offers 
quick method for the correction error without the necessity changing 
the entire computation from the point which the error was made the end 
the solution. Furthermore, the method self-checking many different 
ways, found that the intermediate computations some particular 
point the solution have not been carried out sufficient number signi- 
ficant digits, the defect may remedied using one the checks recom- 
pute the defective quantity. 


INTRODUCTION 


Structural analysis and many other engineering problems occasionally re- 
quire the solution simultaneous equations. Since classical methods solu- 
tion involve small differences resulting from the subtraction large numbers, 
computations must carried many significant digits, thus making the work 
very laborious. Study the effect varying the characteristics structure 
(which involves, among other things, changing the coefficients the unknowns) 
these same methods time consuming make such studies prohibitive 
except the simplest cases. The value direct and simple method for 
finding the effect varying coefficients simultaneous equations therefore 
obvious. 

Since error coefficients can interpreted mean that, effect, 
structure has been analyzed whose characteristics are expressed the incor- 
rect coefficients, the correction such errors becomes problem 
coefficients.” was meet just such problem that the methods outlined 
this paper were developed and later extended other applications. 
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obvious extension the concept involved leads directly new method for 
solving this type equations—for, will shown, merely 
stopping point endless chain variations. 

The theory presented, applies the solution simultaneous equations, 
differs one respect from the classical theories that developed from 
different point view. The usual method solution for this and other types 
equations consists, essentially, what might termed process 
ing is, the complete equations are converted such forms that 
the extraction root the elimination other unknowns will result 
simple equation with only one unknown root. This method, the contrary, 
consists process is, the solution begins with simple, 
reduced equations that are easily solved need solution, and are then 
built up, step step, the required equations converting the solutions 
conform the new equations each step. this process, necessary 
vary the coefficients the unknowns and find the resulting effect the 
values the unknowns. This process direct, and the theory for the varia- 
tion the coefficients makes the method solution almost self-evident. 


DERIVATION THEORY FOR VARIATION COEFFICIENTS 


Each one set simultaneous equations takes the form: 


which the are unknowns; the are coefficients dependent only 
the elastic and geometric properties the structure; and the c,’s are con- 
stants determined the position and magnitude the loads addition 
the properties the structure. 

For any given structure the coefficients are, course, constant, but 
will have many sets values there are loading conditions analyzed. 
has often been erroneously stated that the equations must solved for each 
one these sets constants—a statement that only half-truth. the 
work properly arranged and the equations are solved for one set values 
Cn, the solution for any other set requires comparatively little additional labor. 
practice, the equations are solved routine manner for all sets values 
the same time and the computation influence lines therefore 
great additional hardship. simple illustration such procedure pre- 
sented Table which gives solution for the simultaneous equations: 


+42+6y =|42 


Attention called Col. from which easily seen that set equa- 
tions solved for one set values solution for another set can ob- 
tained simply using the factors and performing the new constants only 
the operations indicated this column. However, the number sets 
values for large, often advantageous solve for sets constants 
such are shown Cols. and The numerical solution for ex- 
ample may then obtained setting and 31, from which (line 
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Any number solutions may therefore found either the 
operations indicated Col. using the factors Cols. and 


VARIATION COEFFICIENTS 


Remarks 
(4) (5) (7) 


should noted that these latter factors are actually solutions Eqs. for 
a,0 and 0,b. the coefficients and were changed, 
would give rise new set equations. the classical methods, solutions 
for the new set equations can found only actually re-solving the equa- 
tions. The factors Cols. and can used avoid this difficulty and 
therefore are termed ‘‘variation coefficients.” 

Table from The new set equations then becomes: 


which may rewritten 


Now let (42 treated new values is, addi- 
tional constants Eqs. for which solutions are found. Substitution 

— 1 1.875 23. 
which 


These values and obviously satisfy Eqs. and are therefore solutions 
these equations. Furthermore, since the variation coefficients are also 


Sub 
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and 


The complete solution Eqs. for all values obtained revision 

Eqs. can changed back 


TABLE 
Eqs. 2—or may any 


other equations—by using the new 
variation coefficients, since obvi- 
ously does not one No. 
set equations was solved re- 
vising another set was solved 
For simplicity, only one coefficient 
will now shown that any number 133 


can varied the same time. 


GENERAL FOR VARIATION COEFFICIENTS 


Let assumed that the equations— 


—have been solved some method, matter what, with results shown 
Table 


Also, suppose that the coefficients Ai, are increased and 
respectively. Following the procedure used varying Eqs. 


Substituting line Table 


and 
) 
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from which 


and 


previously shown (Table 1), solution for any set constants, such 


fore, 


Substituting Eqs. 15a and Kq. 14, 


(16) 
However, since any solution for 
similar manner, substituting Eq. line Table 


and 


Similarly, 

TABLE Three equations with three 
unknowns were chosen for sim- 
plicity. The derivation would 
VARIATION obviously the same for any 


number unknowns. Also, 
striction involved. Revised 
coefficients are being varied are 
given Eq. 17; revised values 
for the solutions all other unknowns are given Eqs. and 20. 
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Generalizing: 


and 


which and represent, the original and revised equations, respectively, 
the unknown whose coefficients are being varied; and represent any other 
unknown the original and final equations, respectively; the amount 
which any coefficient increased; and and are solutions for 

Application basic Eqs. the solution the equations— 


(22) 


—is illustrated Table 


EXPLANATION TABLE 


Table assumed that the set simultaneous equations lines 


(Eqs. 22) has been solved some method, matter what. The solu- 


tions corresponding the four sets constants are given the next three 
lines. Let supposed that desired reduce the coefficients these 
equations the amounts shown brackets the first three lines the 
table. For example, line the coefficient which totals re- 
reduced 2—that is, eliminated entirely. The order variation 
immaterial; arbitrary choice the coefficients will varied first. 

tion the original equations for corresponding set constants 

obtain this solution, necessary only substitute the set 


shown line Cols. and 14. 


Having computed the reciprocal the denominator Eq. the next 
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line are multiplied 
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line reasons that will obvious later, the values are also shown 

With the revised values the revised values the other 
the right-hand side Eq. that has not already been computed. Again 
the same sets constants dn. Therefore, line the 


Cols. 12. 


ing each value (in the original equations) line Eq. 
from corresponding values line give the revised values, line 


63’ 
the values the terms are obtained shown line Cols. 


tracted algebraically from corresponding values line give the revised 


The computation for the variation the coefficients now 
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and lines 10, 11, and are the solutions for the revised set equations— 


the next step Table Eqs. are further revised varying the 
referring again Eqs. 21, y’, and the case may be. The 
computations follow, step step, the procedure described detail for vary- 
ing the coefficients Since the new set equations are now being 
that is, Eqs. 23—their solutions lines 10, 11, and are used compute 
and lines 16, 17, and 18. 

These solutions correspond another new set equations— 


the last step Table Eqs. are revised further varying the 
The new solutions are,given lines 22, 23, and 
24, corresponding the new equations: 


Attention called the fact that each step Table entirely inde- 
pendent and complete. Once the solution set equations known, 
immaterial whether was obtained independently varying the 
cients some other set equations. 


SOLUTION SIMULTANEOUS EQUATIONS COEFFICIENTS 


Table the coefficients have been varied more less arbitrary 
manner. This numerical example illustrates clearly that, once the variation 
coefficients are known for any set simultaneous equations, the coefficients 
can varied will and the order the operations almost entirely arbi- 
trary. Also, the process obviously reversible—which leads directly 
another method for solving such equations. 

Table was assumed that Eqs. were solved some way, matter 
how—the solutions being given lines and series variation 
operations, all the coefficients brackets this original set equations 
were removed, resulting the corresponding revised solutions given lines 22, 
23, and this last result obviously can obtained directly 
inspection. Since the process reversible, the results shown these last 
three lines can varied back the original set equations—or any other 
set equations, for that matter. This reverse illustrated Table 
and thus another method for solving simultaneous linear equations obtained. 
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Furthermore, not necessary begin with lines 22, 23, and Table 
which are duplicated lines and Table solution any set 
equations known—from some previous computation any other way— 
can varied any other set equations desired. 
means that the ‘‘variation result from solving the equa- 
the same for any number unknowns. 


EXPLANATION TABLE 


steps which the reverse that Table the coefficients and parts 
coefficients brackets are considered removed, giving rise the new set 
equations: 


This new set equations can obviously solved inspection; the solu- 
tions shown lines and are merely the repetition Eqs. 26. 
the original set equations, necessary vary the coefficients this 
reduced set until all the coefficients and parts coefficients which have been 
removed are replaced. arbitrary choice and also furnish comparison 
between Table and Table the coefficients are varied the order 
and 

comparison Eqs. with lines and Col. shows that these values 
are the coefficients and parts coefficients which have been removed. Sub- 

From this point the end the computation the procedure exactly 
similar that used Table and needs further explanation. The final 
solution given lines 22, 23, and 24. 


RECIPROCAL PROPERTIES SIMULTANEOUS EQUATIONS 


For ready reference, and also bring out certain relationships, the opera- 


tomary consider and the unknowns and and the constants; 
that is, the equations would set this basis practical problem. 
However, once the equations are solved—considering lines Table 
equally correct, mathematically, consider and the unknowns 
and and the constants. 

From this point view, considering and the unknowns, lines 
and give the solutions for the equations lines and illustrate: 
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should remembered that, after variation operation completed, the 
set resulting equations with their solutions are complete themselves. 
The equations may not those that are finally desired, but they are never- 
theless complete. therefore follows that any set equations merely 
stopping point endless series variations. Lines and for example, 
are the solutions for set equations which the coefficient zero 
the first two and equals unity the third, but which otherwise the same 
may considered the unknowns and the corresponding solutions are 
substituted the results will check the new equations. 

The reciprocity the equations plus the fact that each step varia- 
the results are complete makes the operations self-checking. 
Furthermore, there are many methods checking that some the checks 
may used reduce the numerical work and still leave sufficient number 
for actual verification. 

consideration the fundamental theory equations will bring light 
other ways checking not mentioned. For example, the fact that rows and 
columns determinants may interchanged leads the conclusion that 
other sets equations exist which may obtained “rotating” the 
cients positions consistent with the corresponding determinants. 

Furthermore, practical problem structural analysis, most coefficients 
occur pairs consequence the Maxwell theorem reciprocal deflec- 
tions. The final coefficients” also occur pairs and the same 
order the coefficients the unknowns—furnishing automatic check, 
which may also have other uses besides checking. 

examination Table 6(c) will show that, any line, terms above the 
heavy line are merely repetitions the terms the line immediately above. 
will also noted that the values the solution line are not needed 
until the last step, those line are not needed until the next last step, 
and on, there are more than three equations. Furthermore, once value 
appears the computation, must revised every subsequent step. 
Hence, there were method obtaining these values the time they are 
actually required for the next step, they could omitted from all the previous 
computations, thus saving part the work revising them. This abbrevia- 
tion can achieved using one the checks obtain the omitted values. 

The revised variation operations illustrated Table show that there 
considerable saving work even though partly offset the necessity 
making substitution computation every step-in the process. 


THEORY ABBREVIATED SOLUTION 


The procedure followed Table for solving set simultaneous equa- 


tions may simplified taking advantage certain characteristics this 
type solution. 


October, 1947 SIMULTANEOUS EQUATIONS 1233 


| 
! 
| 
> 
7 
“ul 
ae 
le, 
- 


1234 SIMULTANEOUS EQUATIONS Papers 


Reference Eqs. and 216 will show that the revised value the 
unknown whose coefficients are being varied the particular step 
tion” under consideration, computed first and that also enters the 


TABLE RESULTS TABLE AND TABLE 


Constants 


(2) (3) (4) (6) 


+ 


(b) VARIATION COEFFICIENTS Eqs. TABLE 


+ 


formula for the revised values the other unknowns. The term may thus 
considered the critical one for its step. will shown that, sense, 
solution reduces the computation critical values only—that is, values 
the unknowns that have not been critical previous steps, and that are also 
not critical the step under consideration the one immediately following, 
need not computed. 

Table for example, the solutions lines 16, 17, and correspond 
Eqs. 23, which the values and are obviously determined the first 
two equations and are independent the third. Mathematically, this fact 
becomes evident the computation for and lines and respectively, 


Line L 
No. 
z 0 0 = +3 = +734 0 23. 
76. 
0 y 0 = = + 0 24. 
25. 
coef 
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multiplied zero; obviously the coefficient the third Eqs. 23. 
The corresponding result occurs the computation for and lines 


TABLE 6.—(Continued) 


No. Constants 
(1) (2) (3) (4) (5) (6) (7) (8) 
315 
0 0 13 13 a + 13 b c 
mei ctea 0 [y] 0 = +2 = + 7 a + 9 b ad 63 c 


Unknown whose coefficients are being varied step under consideration. Unknown whose 
coefficients have been varied previous steps. 


Since the values and are independent the third Eqs. 23, all 
the coefficients this equation may neglected the computation 
this point. However, once and are known, the values can computed 
from the third equation direct substitution. therefore follows that the 
values need not necessarily computed the previous steps. The 
values must computed this point, however, since the coefficients 
are varied the next step Table and the critical value could not 
otherwise computed. 

Similarly, the solutions lines 10, and 12, Table correspond 
Eqs. 24, which determined the first equation and independent 
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the other two. Mathematically, this fact again evident the computa- 
tion—the appropriate d,-values are multiplied zero the computation 
for all the coefficients the second and third equations may 
neglected the computation this point—that is, the computation 
for 

The values both and may computed direct substitution, once 
known. necessary compute only the values the 
unknown whose coefficients are varied the next step—in this case 
the next unknown become critical. The revised values can com- 
puted later, already explained. 

The abbreviated procedure may now summarized. Table lines 
and all the coefficients below the lower heavy line each column may 
neglected each step The values u’, which critical, are 
computed Eq. 21a; the values corresponding the unknowns whose 
coefficients have already been varied previous steps, are computed 
and the values v’, corresponding the unknown whose coefficients 
are varied the next step, are computed direct substitution. The 
values the remaining unknowns are neglected because they are not essential 
the final solution, and are not computed until they are about become 
critical the step immediately following. 


EXPLANATION TABLE 


Table the order variation the coefficients the same Table 
and 

line Col. previously explained. With this exception, nothing new 
presented line which gives the first revised value The primes are 
omitted for simplicity—the line for instance, really Since the 
solution this step will correspond Eqs. 24, the values which critical 
the next step, may computed from the second these equations— 

Values line Cols. and subtracted from the corresponding 
explained, not essential compute the revised values this point. 

the second step, variation the only necessary 
use and +-3. With this exception, the computations 
lines and 10, and lines and present nothing new. Since the solution 
this step will correspond Eqs. 23, the values can computed from the 
theses are computed lines and and are summed line 
traction the values line from the corresponding values line Cols. 
and gives the revised values line 16. 

The third and last step presents nothing since and are used 
and there substitution made; the last variation Table exactly 
the same that Table 5.. 
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MULTIPLE VARIATION 


the number unknowns large, the solution still requires considerable 
work—there will many steps variation operations theré are un- 
knowns. further simplification possible such the coefficients 
two more unknowns are varied the same time, auxiliary set simul- 
taneous equations must solved. the coefficients all the unknowns were 
varied simultaneously, for instance, the auxiliary equations would simply 
repetition the original equations. However, since fairly easy solve 
set two even three equations, does save considerable work vary 
two three unknowns The result such procedure is, effect, 
that one large set equations broken into several smaller sets. 
well known that, when similar situation occurs analysis the result 
symmetry, the resulting equations are much easier to. solve. 

Referring Eqs. 11, let supposed that the coefficients and 
however, the coefficients B,, and also are increased amounts 


and respectively. Following the same procedure used for varying the 
coefficients only one unknown: 


Substituting these values line Table gives 


Therefore, 


Similarly, substitution the values for and line Table gives 


and substitution line gives 


obtain and y’, Eqs. and must solved simultaneously; the values 
are then computed directly from 31. 

there were more than three unknowns, the equations for the remaining 
ones would take the same form Eq. 31; for example: 


Generalizing, 
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which both and represent the unknowns whose coefficients are being 
varied, and the amount which each the coefficients being 
varied. The notation otherwise the same that Eqs. and 216, and 
the primes have the same meaning. 

Revised values the unknowns whose coefficients are being varied simul- 
taneously are computed solving auxiliary set simultaneous equations 
such Eqs. 33. The remaining unknowns are then computed 34. 
33a, 33b, and can obviously extended include the simultaneous 
variation any number unknowns and the solution large set simultane- 
ous equations can thus broken into the solution several smaller sets. 
The optimum number simultaneous variations more less matter 
individual preference. The procedure illustrated Table for set four 
equations split into pairs two each. 


EXPLANATION TABLE 


The procedure for the variation the coefficients more than one un- 
known one time similar, general, that illustrated the previous 
tables for variation the coefficients only one unknown. 

Table the set four equations— 


114 


—is broken into pair sets with two unknowns each, indicated the 
the procedure and abbreviations Table suitably modified, are also used 
Table 

Remembering that Eqs. 33a, 33b, and are used instead Eqs. and 
21b, reference Table will clarify the work Table Coefficients below 
the lower heavy line, for example, are omitted from the variation computations 
both cases for reasons previously explained. However, since two sets 
coefficients are being varied the same step, the numerical equivalents 
Eqs. are shown lines two auxiliary equations are solved 
simultaneously inspection—as will explained later—and the result 
given lines and 12. Since the coefficients both and are varied 
the next step, their revised values are obtained this point direct 
substitution shown lines 16, and lines and 20. happens that 
there are other unknowns, but there were, their revised values would not 
have computed this point; they can computed later direct 
substitution when needed. The complete solution for the first step shown 

The next step similar: The auxiliary equations are shown lines and 
22, their solution (by inspection) lines and 25, and part the numerical 
equivalents Eq. lines and final solution for this step, which 
the last this computation, shown lines this the last 
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step and there are other unknowns, there direct substitution 
made. 

over-all check the entire computation may made adding the 
four equations lines Table and substituting, the resulting equa- 


Table The result this substitution (line Col. Table should equal 


Line No. Substitution factors 


286 


The over-all check lines 10, Table similarly made consider- 
lines 35, Table are considered the simultaneous equations, with 

previously stated, the solutions the auxiliary equations lines 
and 12, and lines and 25, Table may found practically inspection. 
Table and its solutions lines 35, can shown that the determinant 
Either determinant can obtained from the other. element one 
determinant equal the cofactor the element the other determinant 
divided the value the latter determinant. The cofactors second 
order determinant are simply elements—in each case, simply the element not 
crossed out. also simple matter find the value the determinant 
itself. The solution the auxiliary equations can therefore found prac- 
tically inspection. 

illustrative example, reference made lines and 10, Table 


where the determinant formed the coefficients the auxiliary equations is: 


(2) (3) (4) (6) (7) (8) 
+18 
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The cofactor the element (which and the element the 


cients lines and 12, Cols. and 

The values shown parentheses lines and 25, Cols. and 10, are 
similarly obtained from the determinant formed from the coefficients lines 
and 22, Cols. and 

both cases, the other values are obtained using the variation coeffi- 
cients obtain other solutions. For example, the set constants from lines 
and 22, Col. (that is— 226 and are substituted for line 24, 


the value line 24, Col. The primes are not indicated Table 
for simplicity. 
Similarly, substitution the constants 217 and 245 for line 25, gives 


The elements this determinant are the variation coeffi- 


line 25, Col. The other values are computed similar manner. 

will noticed that the coefficients lines Table occur 
pairs, they would practical design problem result the Maxwell 
theorem reciprocal deflections. The variation coefficients, lines 35, 
also occur pairs and the same 


practical design, the question the number significant digits which 
the computations must carried difficult one determine. is, how- 
ever, very important consideration from the point view time. The 
method variation coefficients makes the determination the required 
accuracy simple matter. 
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PAPERS 


APPLICATION GEOLOGY TUNNELING 
PROBLEMS 


Geological examination tunnel sites prior and during tunneling opera- 
tions has become established procedure. The geologist should study and ap- 
praise carefully all surface and subsurface features having any bearing the 
engineering problems associated with the drawing the tunnel. Particular 
attention should given the following: (1) Internal structures, physical 
characteristics, and distribution the rocks and overburden; (2) secondary 
structural features such folds, faults, and joints; (3) chemical alteration 
the rocks; and (4) the distribution and flow ground water. 

Suggestions for overcoming difficulties presented caving 
ground and other undesirable conditions tunnels are given, and the causes and 
effects caving selected examples are discussed detail. 


INTRODUCTION 


Geological investigations play increasingly important part the plan- 
ning and completion certain types engineering projects. Geological ex- 
amination dam sites, bridge and building foundations, and tunnel sites prior 
and during construction has become established procedure that has 
yielded substantial even spectacular returns the form safety and reduced 
costs construction and maintenance. 

The writer has spent considerable time consultant studying the geolog- 
ical features several proposed water diversion tunnels, and connection 
with this work has examined many existing tunnels and adits effort dis- 
cover and evaluate the various geological factors that are important econom- 
ical tunneling operations. Mine tunnels yield the most valuable data because, 
general, these tunnels follow cut across ore bodies formed where faulting, 
fracturing, and alteration are maximum. much pertinent in- 
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formation has been obtained from railroad tunnels, drainage tunnels, and ex- 
isting water diversion tunnels. 


THE FUNCTIONS THE GEOLOGIST 


The functions the geologist tunneling operations are manifold. Gen- 
erally, called upon examine the proposed site for tunnel determine 
located the best advantage. Such investigation entails geologic 
study the portal areas and the surface over and adjacent the line the 
tunnel. The results the study may lead one several conclusions: 


(1) The tunnel favorably situated and little difficulty should 
experienced driving it. 

(2) There are geologic indications the presence sections bad ground 
the course the tunnel. (It the geologist’s responsibility indicate the 
probable location and nature the bad ground these sections.) 

(3) The tunnel site unfavorably situated, and the cost driving the 
tunnel will very high not prohibitive. 


the location the tunnel cannot changed conditions demand that the 
tunnel driven spite the unfavorable prognostication, then the geologist 
faced with the responsibility predicting the best his ability the extent 
and character the various undesirable conditions that will encountered. 

Many tunnel sites are areas partly wholly covered concealing allu- 
vium overburden. covered areas the geologist must decide what ex- 
tent and what methods the subsurface should explored. Frequently, 
trenching drilling needed, which should supervised the geologist with 
attention minute detail. Exploratory work this type expensive and the 
geologist generally must satisfied with bare minimum information even 
though the facts obtains not begin solve the host unanswered prob- 
lems that are associated with regions thick widespread blanketing de- 
posits. Under certain conditions geophysical investigations yield pertinent 
information, and the geologist’s duty determine which geophysical method 
methods should employed and what extent. 

Once the tunnel started, the geologist should make periodic examinations 
the completed sections and the face. should maintain detailed written 
account and map showing the geological features the tunnel well the 
day-by-day progress. The study the and structures encountered the 
tunnel may lead increase the accuracy with which the geologist may 
predict the changes the rocks structures beyond the tunnel face. The 
geologist may give advice the type and the amount support needed 
bad sections and the most effective manner drilling and blasting the 
face—giving due consideration the relationship between the structure the 
rock and the ease with which mined. the event that unexpected 
culty encountered because extremely bad ground, may offer counsel 
the best method overcoming the 

The importance detailed geologic map the tunnel cannot over- 
emphasized. This map not only should show the geologic structures and rocks 
encountered, but also should serve lasting record indicate the location 
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bad sections, which may require attention many years after the tunnel com- 
pleted. 

The geologist should act advisory capacity the engineer charge 
the tunneling operation. Too often the geologist’s advice does not carry the 
weight should because either the engineer does not understand what the geolo- 
gist attempting feels that the geologist’s predictions are too vague 
and indefinite. However, this situation gradually becoming less serious and 
the engineer and the geologist are reaching amicable understanding based 
mutual respect and profitable cooperation. geologist should not attempt 
assert authority any problems facing the engineer except those which bear 
directly upon, are the outgrowth of, the geological study. the other 
hand, the engineer generally not qualified training experience make 
critical appraisal the findings trained geologist and must accept the 
advice the geologist the best that can offered. 


CHARACTERISTICS CONSOLIDATED Rocks 


Consolidated rocks are those which the constituent particles fragments 
are bound together strongly enough produce aggregate with moderate 
considerable strength. The rock mineral components may held together 
with cementing material certain clastic sediments, the contacts the 
mineral grains may interlocked intermeshed igneous rocks and 
most metamorphic rocks. 

Certain rocks are massive isotropic. The physical characteristics 
these rocks, from the points view the the tunnel engineer, 
are the same all directions within the rock. There preferred orienta- 
tion the mineral particles, and the ease tunneling such rock the same 
irrespective the direction the tunnel. this category belong igneous 
rocks such certain types granite, sedimentary rocks such certain lime- 
stones, and metamorphic rocks such certain marbles and quartzites. Com- 
pletely isotropic rocks are relatively rare. 

Most rocks possess directional structures result the parallel sub- 
parallel orientation many most the mineral particles. Such rocks are 
anisotropic and show distinct variations physical properties different direc- 
tions. Consideration the nature the anisotropy rocks important 
the tunnel engineer because has direct bearing the cost driving the 
tunnel. 

igneous rocks layered linear structures are commonly present. the 
rock contains minerals such tabular feldspar mica, there 
tendency for these minerals rotate into parallel subparallel positions 
during the intrusion and crystallization the rock. The resulting formation 
“flow layering” Rocks which show well-devel- 
oped platy structure have gneissoid aspect and commonly contain flat, lentic- 
ular, dragged-out inclusions foreign rocks called Generally, 
such rocks show strong tendency split break along planes paralleling the 
flow layering, and they part break with difficulty directions right angles 
the layering. 

Rocks containing minerals prismatic needle-like habit tend develop 
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This formation results from the parallel subparallel alinement 
the prismatic minerals during the intrusion and crystallization rock. 
not uncommon find both platy and linear structures the same rock. 


SURFACE PARALLEL 
FLOW-LAYERING 


STRUCTURE 
FLOW LAYER 


both are present, the needles prisms lie the plane the platy structure. 
Linear structure may emphasized the presence elongated, spindle- 
shaped inclusions foreign rocks. Igneous rocks which linear structure 
well-developed generally break easily parallel the lineation and break with 
difficulty right angles the linear structure. 

Fig. block diagram igneous rock showing platy and linear struc- 
tures. When the geologist describes these features, generally gives the 
strike and dip the and the direction the trace the linear 
structure horizontal plane. These data permit the easy calculation the 
pitch angle the flow lineation. Where flow lineation and flow layering are 
steeply inclined, lineation slight practical importance the tunnel opera- 
tion. However, the direction the tunnel with respect the flow layering 
merits careful consideration. tunnel being driven across the flow layering 
probably will have the following characteristics: The face will advance re- 
latively rapidly because the blasting operation the rocks will break cleanly 
and easily away from the planes weakness paralleling the flow layering. The 
walls, floor, and roof will tend rough and blocky. the other hand, 
the tunnel parallel the strike steeply dipping flow layering, the forward 
progress the tunnel will slowed because the tunnel parallel direction 
maximum strength and toughness the rock. The walls the tunnel 
will probably relatively smooth; the floor and roof, more less irregular. 
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The irregularity the walls increases with decrease the angle dip 
the layering. 

metamorphic rocks the layering commonly called and 
well developed gneisses and schists. typical gneiss consists alternating 
layers granular rocks differing mineral composition and physical properties. 
schist, although foliated, relatively homogeneous its mineral composition. 
Slate fine-grained rock with highly developed tendency split into thin 
plates. Foliated metamorphic rocks generally tend split most easily along 
planes parallel the foliation. 

The effect foliated metamorphic rocks the progress the tunneling 
operation even more pronounced than the effect layering igneous rocks— 
for the same reasons and with the same results. tunnel right angles the 
strike the foliation general constructed much more rapidly than one 
following the foliation, although the tunnel likely much more blocky 
and irregular. tunnel following the strike steeply dipping foliated rocks 
tends cave from the roof. The smoothness the walls lessens and the 
tendency the walls cave increases the dip the foliation decreases, 
the tunnel following the strike the foliation. Certain gneisses and schists 
are complexly folded and crumpled; and, where this true, the difficulty pre- 
diction the relation between the ease tunneling and the rock structure 
increases. 

one tunnel Boulder County, Colorado, the cost driving the tunnel 
was much 25% greater the sections paralleling the foliation than the 
sections cutting across the foliation. the overbreak the latter was 
about 20% greater than the first instance because the irregularity break- 
age along the walls and roof. 

Typical conditions are illustrated Fig. which shows one tunnel cutting 
across steeply dipping, foliated 
rocks, and another following the 
foliation 
layers. The latter condition 
should avoided possible. 

sedimentary rocks bedding 
stratification may con- 
spicuous, but the strength 
individual layer parallel the 
bedding not much different 


from the strength across the 
bedding. However, the rock 
thinly stratified consists 


(a) Cutting Across Strike 
Parallel Strike 


alternating thin layers different composition (for example, thin-bedded alter- 
nating sandstones and shales), acts much the same manner foliated 
metamorphic rock. sedimentary sections where shale abundant, there 
pronounced tendency toward caving and slumping, and considerable sup- 
porting timber may required the tunnel driven. 

tunnel driven through inclined layered rocks right angles their 
strike, additional factor should considered—that is, the direction which 
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the tunnel being driven relative the dip the layers. Preferably, the 
tunnel should driven that the layers will dip away from the heading. 
this done, drilling and blasting operations will expedited because the clea- 
vage the rocks tends produce overhanging face rather than one which 
sloping back toward the already completed section the tunnel. This sug- 
gestion applies equally well tunnels being driven faulted strongly jointed 
rocks. general, the tunnel should headed the down dip direction 
the planes weakness the rock. 

These conclusions purposely are generalized. attention paid the 
modifying effects fracturing, alteration, and other features because these are 
considered separately. not claimed that the generalizations noted are 
without exception because various conditions have modifying effects. 


AND 


Folding layered rocks requires consideration the geologist because, 
folded structures are present, correct interpretation permits prediction the 
attitudes the layers different sections the tunnel. foliated metamor- 
phic rocks, particularly, the prediction variation the attitudes the layers 
the tunnel level may call for detailed surface studies and all the ingenuity, 
skill, and experience the command the geologist. 

Actual incipient fractures are present all rocks. Fractures along which 
there has been little movement are called “joints.” are frac- 
tures along which there has been appreciable considerable movement. 
Joints result from various causes, most important which are expansion 
contraction the rock due thermal, chemical, mechanical causes, frac- 
turing the rock response deformation stresses accompanying move- 
ments the outer part the earth’s crust. Faults show relative displacement 
the rocks either side the fault plane. Fault displacements range from 
few inches several miles. fault may very narrow and sharply defined 
may involve rocks zone hundreds feet wide. 

Joints and faults are major cause bad, caving ground tunnels and 
deserve careful consideration. Joint surfaces commonly are rough and irreg- 
ular coarsely granular rocks, particularly the joints developed result 
shrinkage decompression. Blocks bounded such joints rarely cause 
trouble because friction along the rough surfaces prevents disengagement 
blocks from the walls roof the tunnel. However, the irregularities the 
joint surfaces have been smoothed slight movements between adjacent 
blocks the joints are smooth for other reasons, the tendency cave in- 
creased. Moreover, water has permeated closely jointed rocks and has re- 
sulted alteration and formation clay minerals along the joints, partic- 
ularly bad, raveling type ground produced. 

careful study rock mass may show that several joint systems are pre- 
sent and, moreover, that certain the joint systems are better developed than 
others. possible, the tunnel should located that its direction right 
angles the most prominent, most steeply inclined joint system. This course 
should chosen because drilling and blasting the ease breaking the rock 
the heading away from joints “slips,” advantageous. Moreover, 
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tunnel follows major joint system, the danger overhead caving con- 
stantly present. igneous rocks showing flow layering has been observed 
that very prominent joint system commonly parallels the flow layering. 

Engineers and geologists have developed real fear faults tunnels, and 
justly so. course, tunnel should not driven through along active 
fault; but, fortunately, most faults are not active. The fear encountering 
fault tunnel should justified tempered detailed study the 
faults the region where proposed drive the tunnel. Some faults offer 
difficulty, distressing problems. 

Among the faults which cause little trouble are those which show 
sharp, clean offset the rocks and are not accompanied decaying alteration 
the wall rocks. Certain faults fault zones have served localize the 
circulation mineralizing solutions and examples are known where the sub- 
stances deposited from the solutions have produced tight, well-cemented rock 
stronger than the original rock. However, mineralizing solutions may pro- 
duce the opposite effect, and broad zone within and either side the fault 
may converted into rotten, swelling, squeezing ground. Alteration 
feldspars clay minerals micaceous minerals often produces weak, slippery 
rock which will cause considerable trouble. 

Fault movements commonly develop selvages weak, pulverized rock flour 
wide fault zones where the gouge present weaving, criss- 
braided fractures, extensive caving likely result. some 
faults the rock are brecciated into loose, porous aggregates. Additional move- 
ments may cause the rounding the angular fragments. Generally, such 
faults are filled with water and the combination loose rock and water under 
hydrostatic pressure has caused more than one tunnel rerouted even 
abandoned. Such fault zones were encountered the Elton tunnel under the 


(a) 


(a) Rough and Irregular, Minimum Caving 
(b) Smooth Surfaced, Maximum Caving 
(c) Smooth and Rough, Considerable Caving 


Oquirrh Range Utah and the Leadville drainage tunnel Colorado. Fig. 
shows typical joints and Fig. shows several common types faults. 

summary, may stated that study the joint systems and. faults 
area where proposed drive tunnel should included the geo- 
logical investigations over the line the tunnel and the area adjacent it. 
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ALTERED AND SWELLING GROUND 


the preceding section brief reference was made altered rocks the 
vicinity joints and faults. was suggested that alteration might produce 
effects serious nature tunneling operations. 


(a) Clean Sharp Break, Alteration Wall Rocks 
(b) Sharply Defined, Subsidiary Fractures, Pervasive Alteration 
(c) Sharply Defined, with Subsidiary Fractures, Localized Alteration 
(d) Braided, with Subsidiary Fractures, General Alteration 
(e) Open, Porous, Filled with Breccia 


One the minerals most susceptible alteration rocks feldspar. 
Attack feldspar reactive solutions surface origin deep-seated 
derivation may cause feldspar converted from fresh, physically strong 
substance weak, pulverulent aggregate clay minerals micaceous min- 
erals, and fine-grained silica. This change, essentially chemical one, gen- 
erally most pronounced and near zones intensive jointing faulting, 
where altering solutions move with maximum freedom. 

Grain size controls the amount alteration certain extent. Fine- 
grained feldspathic rocks general alter more readily than coarse-grained 
rocks, probably part because fine-grained rocks may more closely jointed 
than coarse-grained rocks and part because the small size each grain pro- 
motes more rapid, pervasive alteration. Fault movement may pulverize the 
rocks along the fault and cause the development tiny fractures 
the rocks the vicinity the fault, thus producing ideal conditions for rock 
alteration. Moreover, faulting causes the minerals the rock along the fault 
strained and, general, strained minerals decompose more readily than 
unstrained ones. 
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the altering solutions are deep-seated origin, alteration may not change 
the strength the rock. For example, hydrothermal solutions commonly de- 
posit condiderable silica the wall rocks fissures, and the silicified rock will 
stand the tunnel well even better than the original rock. Sulfides 
deposited fractures the wall rocks along fractures tend oxidize 
exposure that mineralized areas containing abundant sulfides have tend- 
ency swell slowly and slough after exposure air. 

rock the tunnel cross section reduced the plastic flow very weak rocks 
response the pressure the overlying rocks. Rocks that act this man- 
ner generally are highly altered, are rich clay minerals, and contain consider- 
able water. The plasticity increases with the water content. 

swelling ground there actual increase the volume the rock re- 
sulting from hydration chemical reaction from combination these pro- 
cesses. Certain clay minerals and anhydrite have remarkable affinity for 
water. rocks containing these minerals are encountered the tunnel, water 
from the atmosphere and other sources tends absorbed and swelling re- 
sults. Tunnel sections have been studied where the rock when first encoun- 


tered appeared reasonably fresh and unaltered; but, after being exposed 


for short time, began swell and eventually was converted into sticky, 
plastic mass clay minerals.. This same material broke down the dump 
into almost impalpable powder. Probably, the rock was partly altered 
depth hydrothermal solutions, but the confining pressure the overlying 
rocks did not permit complete hydration and expansion. When the pressure 
was released the tunnel, the rock, metastable state, rapidly reached 
state chemical equilibrium response the release the confining pressure. 
Pressures high tons per have been encountered swelling ground. 

The most spectacular example that has come the writer’s attention was 
the swelling dike fine-grained porphyry cut tunnel gold mining 
operation. Within few days after exposure, the originally fresh looking rock 
had almost completely filled the tunnel with sticky, relatively cohesive mud. 
the Moffat tunnel water bore (in Colorado), deceptively fresh looking 
granite irresistibly splintered Oregon fir timbers, Many similar 
examples may cited. 

Rocks metastable condition may swell chemical reaction well 
simple hydration. established principle that the release pressure 
rock tends produce new mineral assemblage which equilibrium with 
the new environment. Fortunately, however, chemical changes response 
pressure changes most rocks are slow that they are importance 
tunneling operations. 


UNCONSOLIDATED 


The range composition, texture, and structure unconsolidated rocks 
extreme. Texturally these rocks grade from the finest muds coarse aggrega- 
ates such are found rock slides and bouldery glacial deposits. Unconsoli- 
dated rocks are characterized the lack interlocking contacts between ad- 
jacent grains the lack strong, binding cement matrix. Some rocks, 
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as, for example, certain sandstones, have been converted into loose aggregates 
mineral grains through the removal the cement dissolution. 

tunnel operations unconsolidated rocks generally are encountered the 
portals, where the bedrock blanketed with rock waste alluvium. Loose 
sand and gravel, wet dry, present considerable obstacle does loose mate- 
rial consisting large, round, angular boulders with cementing material. 
The difficulty tunneling through unconsolidated rocks consisting partly 
largely very fine-grained clayey materials increases direct proportion 
the amount water present. Dry slightly moist clayey aggregates stand 
fairly well portal areas; and, despite slumping tendency, the walls and roof 

_generally will stand long enough permit support and lining set before 
caving starts. the water content increases, however, the plasticity and mo- 
bility the clay mineral aggregate increase, and considerable difficulty may 
experienced excavating and supporting such ground. 

One the most vexing situations arises when tunnel driven through 
water-saturated volcanic ash beds. wet ash beds mobile, gelatinous slurry 
generally present which held with only the greatest difficulty. the 
other hand, the material compact and dry, offers little trouble except 
that may swell exposure water. 

railroad tunnel near Leadville, Colo., was driven through dry, bouldery 
glacial till without mishap. The clayey matrix the till was dry enough 
have considerable strength and serve reasonably strong binding medium. 
the other hand, near-by drainage tunnel driven through water-saturated 
glacial deposits encountered many difficulties that part the tunnel had 
abandoned and the cost the tunnel greatly exceeded the most pessimistic 
cost estimates. 


WATER 


Prediction ground-water conditions tunnels one the most baffling 
problems confronting the geologist. Any prediction must based con- 
sideration the porosity and permeability the rocks and the availability 
ground water. Rapid movement water rocks generally determined 
the presence connected open joints, fissures, faults. common ex- 
perience the tunnel engineer drive through fault containing little 
water and then encounter narrow crack section jointed ground from 
which water issues great volume and with destructive velocity. This phe- 
nomenon explained in-terms relative permeability. Many faults, but not 
all them, contain abundant gouge which effectively prevents the rapid 
movement water. Joints and small fractures, however, generally are not 
clogged with impermeable material and will permit the easy flow water. 

Faults characterized abundant breccia and open spaces generally are 

permeable and ample water. Tunnels driven limestone are likely 
encounter underground caverns filled with water. these caverns 
‘commonly determined the presence fractures which have permitted 
the easy circulation water with consequent dissolution the limestone. 

knowledge the configuration the ground-water table useful; but, 
general, the accurate determination the depths the upper surface the 
saturated zone obtained only extensive drilling geophysical study. 
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compact, consolidated, relatively impermeable rocks, the configuration the 
ground-water table great significance, because impermeable rocks 
not permit the circulation water. alluvium any permeable uncon- 
solidated rock, however, the story quite different. 


Tue TuNNEL Cross 


Too little attention has been paid the relationship between the tunnel cross 
section and the ease difficulty with which rocks the wall and roof cave. 
certain tunnels, much time and expense can saved careful consideration 
given details rock composition and structure and methods pre- 
venting reducing the hazard caving ground. the purpose this 
section briefly describe the nature the caving and suggest remedy 
certain selected examples. The generalizations not apply tunnel sections 
squeezing swelling ground. 

The caving tendency present all tunnels, but the extent caving de- 
pends the strength and structure the rocks. well known that 
average caving ground the caving tendency does not extend for great distance 
above the tunnel unless the ground almost completely without strength and 
the free running type. Observations hundreds caved sections rock 
tunnels indicate that the caved area commonly extends height above the 
roof nearly equal the width the tunnel. the rock highly fractured 
altered, the caved section may involve mass rock above the tunnel roof 
with height equal much twice the width the tunnel. the caving 
proceeds, the tunnel tends assume cross section similar that Gothic 
arch. When the arch finally developed, acts support the walls and roof 
the tunnel and prevent further col- 
lapse. Except shallow tunnels 
tunnels free running ground the caving 
effects almost never extend the earth’s 
surface. Insupporting the average tun- 
nel roof, then, the maximum strength 
the supporting structure, even the worst 
kind caving ground, general need not 
exceed that necessary support arch- 
shaped block ground above the tunnel 
roof having maximum height equal 
twice the width the tunnel. For aver- 
age rock the potential dead weight the 
roof the tunnel considerably less than 
very bad ground. However, thisdead How 

weight needs support only when the rock 

has small strength. massive, un- 

altered rocks, the potential dead weight the roof rocks more than compen- 
sated for the inherent strength the rock, the rock load reduced zero, 
and little support necessary. Also, the caving not allowed begin, 
even bad ground, the ground may held relatively weak supporting 
structure because the caved arch can develop only progressive slumping 
starting near the roof the tunnel and extending upward. 
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The shape assumed the arch resulting from caving certain extent 
function the geologic structures the roof rocks. general, the height 
the caved area greater faulted and highly fractured altered rocks 

than rocks which these features are 
absent poorly defined. Fig. shows 
the effects faults the upward ex- 
tension caving ground two selected 
examples from tungsten mines 
County, Colorado. 

driving tunnels through caving 
ground special attention should paid 
the tunnel cross section. Much time 
and effort can saved arching the 

strength there the caving rock until 
place. Obviously, not generally fea- 
sible drive tunnel which the height twice three times the width; but, 
nevertheless, profitable arch the roof somewhat and, possible, make 
the arch resemble shape the Gothic-like arch that will result from caving. 
nonlayered rocks the tunnel cross section made symmetrical possible; 
but, layered foliated rocks, asymmetrical arch may serve the purpose 
better advantage. 

Fig. shows two suggested tunnel cross sections for caving ground. 
course, spite all precautions and spite the maintenance arched 
cross section, some ground will continue cave. However, this caving ten- 
dency offset even slight degree, the precautions suggested will pay off 
the more rapid forward advance the tunnel and the need for less tem- 
porary permanent supporting structure. course, the rock the tunnel 
strong, cross section preserved which will most economical considering 
the purpose the tunnel and the immediate requirements the engineers 
driving it. 

may argued that altering the cross section the tunnel response 
changes the rocks and rock structures might the long run cost more than 
attempting preserve standard section and taking chance caving 
ground. However, éven slight concession the direction tailoring the 
tunnel cross section the rocks will pay dividends. Thus, from practical 
point view, the actual cross section employed necessarily must com- 
promise between what wanted the geologist and what desired the 
engineer. 

squeezing swelling ground circular cross section with closely spaced 
circular supporting ribs appears overcome the difficulty best advantage. 
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PAPERS 


THE GENERAL PRINCIPLES HIGHWAY 
TRANSITION CURVE DESIGN 


The theoretical and experimental background the published theory 
highway transition curve design explored this paper. The writer shows 
(1) that the theory fallacious its disregard the superelevation, and 
(2) that the experimental evidence produced for questionable. 

modified theory, called lateral force theory,” outlined. Experi- 
mental evidence produced show that the length the transition curve 
not matter great importance, and the proposal made that should 
made fixed proportion the total length the curve. also shown that 
the transition need not necessarily spiral form. 

Simplified equations for the solution the compound curve are given, 
based the radius the arc parameter, together with table spiral 
functions. There also brief description the proposed method curve 
design, illustrated worked example. 


INTRODUCTION 


designing transition curves the fundamental principle the method 
described John Oran Eichler, Assoc. ASCE, and Howard is, 
the words the 


that the rational way achieve this transition [from tangent 
curve] make the rate gain the ‘centrifugal ratio,’ P/W, con- 
stant throughout the curve. This equivalent the demand that the 
driver sharpen his curve turning his steering wheel constant 


The symbol this context denotes the centrifugal force. this theory 
also postulated that the rate change centrifugal force, for which 
symbol will used, should the same for all curves. Mr. 


comments are invited for immediate publication; insure publication the last dis- 
cussion should submitted March 1948. 


Surveyor, Dorsetshire, England. 


Examination the Transition Spiral,” John Oran Eichler and Howard 
Eves, Transactions, ASCE, Vol. 111, 1946, 959. 


972. 
Ibid., 960. 


Practical Method for the Improvement Existing Railway-Curves,” Shortt, Proceed- 
ings, Inst. E., Vol. 176, 1909, 97. 
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Shortt, England, suggested that this rate for should per sec’, claim- 
ing that this the maximum rate that will pass unnoticed. The value has 
been accepted others, although the United States values per 
and per have been proposed. 

This theory gives the same length for transition curve whether super- 
elevation present not, although obvious that this factor must have 
some effect car. 

Consideration the foregoing from the viewpoint first principles 
once raises the question whether the rate change centrifugal force 
the important effect for transition curve design. find answer 
necessary consider the forces acting car traveling canted road. 


road canted angle The angle, and all forces, are considered 
positive drawn. The centrifugal force The reactions the road 
the car, parallel and normal the surface, are designated and 
respectively, and are called the lateral and normal forces. Then, resolving 
forces and eliminating the weight the car, 

and 


which the acceleration due gravity. The two quantities, and are 
called the lateral and normal ratios, respectively. They are pure numbers, 
provided the units used are mutually consistent. 

The normal force, very little greater than under road conditions. 
Its increase due curvature very small increase already existing 
force, acting direction that cannot 
affect the stability the car the com- 
fort the passengers. Therefore, may 
neglected. 

The lateral force, the key the 
situation. the sum the components 
the centrifugal force, and the weight 
the car, along the surface the road and 
the force which tends cause the car 

skid overturn, which causes the 

passengers discomfort urging them side- 

their seats. also the force that 

applied the driver through the steering wheel, the force 
must overcome turning his wheel make the car travel around curve. 

The full implications the lateral force factor road design were 
dealt with some length previous but the most important must 


Curvature and Leeming, Road Paper No. Inst. E., October, 
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change curvature and that one change should never occur without the other. 

certain speed, called the the lateral ratio equal 
zero. this speed the car steers itself around the curve, provided the 
steering properly adjusted, and the driver need not touch the wheel all. 
The centripetal force acting still, but applied the car the super- 
elevation. the curvature and cant are changing together due proportion 
(on well-designed transition curve, for example), the wheel can seen 
turn itself the proper speed held. Alternatively, the driver reaches the 
hands-off speed superelevated circular arc, may then steer the car 
zigzag course varying the speed without touching the wheel, which again 
will turn itself. 

any speed different from the hands-off, the driver must exert some force 
follow the curve, but the direction well the magnitude this force 
will decided whether the speed more less than the hands-off. 
less must steer away, more must steer toward, the center curvature. 
Thus, the presence the superelevation does make very vital difference 
the driver, and also the passengers, who are affected similarly. Therefore, 
the rate change lateral force that the important result and should 
kept constant transition curve. 

significant effect follows from this consideration. ap- 
proximates unity. Then, any constant speed, the requirement for main- 
taining constant hands-off speed that 


which equals v?/g. For this condition circular arc not essential. Any 
curve for which possible keep the curvature and superelevation pro- 
portional will do, including the conventional transition spiral. 


Similarly, the condition for maintaining constant rate change latera 
force, constant speed 


spiral, correctly superelevated, will provide this condition, but its use 
not essential and other types curves can used. Before discussing this 
matter, however, necessary decide whether the rate change force 
important itself. 

The principle this section may summarized another manner. As- 
sume the car traveling constant speed well-designed transition 
curve spiral form. Then the rate change direction plan will 
proportional the rate change centrifugal force, but the change will 
made two agencies. The first the rate change the supereleva- 
tion. This solely function road design, the driver has nothing whatso- 
ever with it, and most cases will unaware it. The second agency 
the application force the driver the steering wheel, and this pro- 
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portional the rate change lateral force. his turn, the road engineer 
has nothing with this. 

The whole essence the theory put forward this paper that this latter 
effect must studied and allowed for design, because may cause varia- 
tion between theory and practice. However well designed curve may 
according spiral principles, the driver does not follow constant pro- 
portional the value which the spiral has been designed, not 
following the design curve. surely logical that the curve should laid 
out conform the driver’s practice; and, prefers, and habitually uses, 
some method turning his wheel other than the constant rate expressed 
constant then the type transition should chosen conform with his 
habits—even this curve not spiral form. 


THe TURNING THE STEERING WHEEL 


the rate application force the driver the steering wheel de- 
fined rate turning the the foregoing principle raises two 
further questions fact, namely: 


Does the driver turn his wheel constant rate any given transition 


curve? 
Does turn the same rate all transition curves? 


These two questions are not the same, although they have not hitherto been 
separated. The first decides the shape the transition curve, but does not 
involve any condition its length. The assumption that this question can 
answered affirmatively has formed the case for the spiral curve, although 
seen that other curves can also satisfy the condition. The second part, 
answered similarly, that all transition curves must designed pro- 
vide the same rate change force, some chosen speed. The assumption 
that both questions can answered affirmatively forms the basis the pre- 
vious theory. 

Examination the literature setting forth the case for this affirmative 
answer the two questions shows that rests entirely Mr. Shortt. 
other writers have merely repeated his conclusions; but his experiments were 
made uncanted railway tracks, and some doubts may felt once 
examining the original paper. Mr. Shortt stated that had found that 
rate change centrifugal force per sec? was the maximum which 
would pass unnoticed. This statement must involve the idea comfort, and 
has been read the past most people. However, impossible for 
observer analyze his feelings delicately and accurately able 
decide beyond doubt that the rate change force affected his comfort, 
and not the force itself. Nevertheless, there must some value the force 
that will cause discomfort, and the possibility that some lower value might 
cause some lesser degree discomfort obviously present. Mr. Shortt did 
not even discuss this possibility, nor did state how arrived his conclu- 
that the reader could judge for himself. 

Royal-Dawson’ Mr. Shortt’s principle roads, without ex- 


7**Curve Design,” Royal-Dawson, Spon, Ltd., London, 1932. 


sti 
ele 
mi 
4 D: 
‘ 
SCé 
of 
tio 
lis] 
cal 
sor 
ans 
Th: 
esse 


October, 1947 TRANSITION CURVE 1261 


perimental proof the validity the transfer, and called the rate 
‘natural rate turning.” this case the difficulty the limit the value 
the force itself, previously referred to, was avoided the expedient stipu- 
lating maximum value for the centrifugal ratio, usually 0.25. The super- 
elevation was specified some fraction the centrifugal ratio, and was then 
stated that the natural rate turning should adhered to, provided that the 
centrifugal ratio did not exceed this amount some 

This argument not very convincing, quite apart from the fact that 
contains the fundamental fallacy using and the centrifugal ratio, ignoring 
the superelevation. The theory gives the same length for transition curve, 
matter what the value the superelevation. However, Mr. Shortt himself 
stated, his original that curve which cannot given super- 
elevation more need easing the ends than one which can super- 
elevated,” thus implying that superelevation does affect the problem, although, 
mathematically, and his followers ignored this principle. 

There has also been large school among British engineers, both railway 
and highway, that has not accepted the principles Messrs. Shortt and Royal- 
Dawson practical grounds. Finally, Warren® carried out some ex- 
periments which threw doubt the validity the rate change force 
the criterion for transition curve design. experiments were small 
scale, and there was some question the accuracy the work, 
was derived from the plotted path the car. 

Thus, was thought that there was need for experimental work the rate 
turning the wheel. The writer, therefore, undertook this work, collabora- 
tion with Black. preliminary report the results has been pub- 
but summary given herein, which includes later work. 


REASONS FOR THE EXPERIMENTAL PROCEDURE 


The experimental technique had based the previous theory, be- 
cause, despite its shortcomings, this held the field, and because had found 
some measure acceptance. Therefore, the experiments were designed 
answer the two questions: 


(1) How drivers actually take curves? 
(2) What the criterion of-comfort curves? 


The two questions Section are included the first these latter two. 

Since the force the driver applies the wheel proportional the best 
method investigating these questions measure directly possible. 
That part the centripetal force applied through the superelevation (that is, 
sin being independent the driver and imposed him the road de- 
signer, irrelevant the purpose. make the conclusions general, also 
essential that the tests carried out under the most natural conditions pos- 
sible, avoid distractions the driver that would.tend influence his method 


Transition Curves,’’ Warren and Hazeldine, Journal, Inst. Municipal 
County Engrs., London, March, 1939, 1021. 


Curvature and Superelevation: Experiments Comfort and Driving Practice,” 
Leeming and Black, ibid., December, 1944, 137. 
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driving, and also that the tests should include many cars and drivers 
possible. 

The method used Mr. Warren was trace out and survey the path 
the car. This procedure, however, open number objections. the 


first place inaccurate, because involves some process double differentia- 


tion the path, determine the centrifugal ratio, from which the super- 
elevation must deducted findj. laborious every way, and involves 
considerable apparatus well technique which must infallibly distract 
the driver. Finally, the method can only useful limited number 
curves because the labor involved—resulting from the fact that detailed 
survey the curve necessary, since the superelevation must known. 
Thus, the answer obtained would limited some extent, and would therefore 
less general application. 


EXPERIMENTAL PROCEDURE 


From the considerations mentioned, was decided that the ideal method 
would use recording accelerometer. instrument measures di- 
rectly, and can determined simply measuring the slope the line 
the record. Accordingly, such device was designed and built, which, although 
simple make with the materials available wartime England, could fit into 
most cars. 

The accelerometer consists 10-lb weight hung springs from the base- 
board strong wooden case. When the instrument placed athwart the 
car, under the influence the lateral force the weight shifts over against the 
springs, moving stylus the end long magnifying lever. This records 
the movement the weight strip metallic paper traveling over drum 
driven phonograph motor. fixed stylus provides datum line. 

This instrument was carried cars traveling ordinary journeys. Before 
starting, was emphasized that the driver should not vary his usual driving 
methods any way. When the car approached any curve thought suitable, 
the operator started the motor; and, for each curve, recorded its direction 
curvature, the average speed the car, and the comfort sensations himself 
and the driver. 

the subsequent analysis, thé individual curves were identified, and 


smoothed off average straight lines this could reasonably done. 


more than 90% the curves was found that the rate change force 
could represented sloping straight lines, and the curve whole two 
such lines, with without length straight line parallel the datum— 
that is, constant them. Thus, the assumption that the driver 
turns his wheel constant rate any given transition seems justified. 

soon appeared that there was considerable variability driver’s prac- 
tice. Records made the same curve with the same car and driver dif- 
ferent times showed very marked differences even the same speed, and 
varied between such wide-limits make quite obvious that there can 
natural rate turning. other words, the answer question 
Section that drivers not always turn their wheels the same rate. 
The variability also made essential use statistical methods, that 
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extensive statistical analysis was started and large number records were 

conclusive evidence could found wide differences practice be- 
tween drivers, considering the general variability and the inevitable differences 
speeds. Most the records were taken, perforce, the British type 
light car, but fair range commercial vehicles was also included. Among 
these latter was the largest and heaviest rigid wheel base vehicles made— 
eight wheeler—loaded with tons bricks. England its speed 
limited law miles per hr, but allowing for this did not seem show 
any marked difference from the smaller cars, nor did the one American car 
was possible include. was not possible take records articulated 
vehicle, nor fast, heavy, passenger vehicle. 

marked tendency was found for vary with the value reached 
any curve. This tendency can expressed the regression equation: 


which the base natural logarithms, and the mean value 
observed, feet per This expression refers entrance transitions, 
and independent speed. The scatter was moderately wide, and the equa- 
tion represents only trend. Another way expressing this principle that 
the driver turns his wheel faster entering sharp curve than entering 
flatone. Values low per sec* were comparatively observed, 
and then only curves the lower ranges 

Another fairly marked tendency was for the driver use definite pro- 
portion the curve transition, with over-all mean two thirds. The 
use the all-transitional curve was comparatively rare. There was also 
marked tendency for the exit transition slightly shorter than the entrance 
one, although the difference not great enough practical importance. 
This latter effect probably related the caster action the steering, which 
helps the driver the exit, and opposes him the entrance. The propor- 
tion transition tended vary with the value reached, and with the speed 
the car and the length the curve, according the regression equation: 


this expression, the mean proportion transition observed (that is, 
the length the two transitions combined, divided the total length the 
curve) the speed, miles per hour; and the total length the curve, 
feet. For 0.1, Eq. gives values ranging from two thirds one 
half. should, however, viewed with caution. Because wartime re- 
strictions England, the data from which Eq. was calculated not include 
very large samples curves taken speeds higher than miles per hr, 
although, since the completion the regression analysis, more records 
higher speeds have been made, and appears unlikely that the conclusions 
will materially altered. 

The curves were also analyzed class road, including right-angle junc- 
tions and rotary intersections. This procedure gave rough classification 
the deviation angle the curve, which appears have little effect. This 
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analysis cannot made with any confidence measuring the deviation angle 
the curves themselves, because not possible decide the driver’s 
angle approach the curve, and the records gave reason believe that this 
varies fairly widely. 

recording comfort, attempt was made separate the due 
the force its rate change. The observer was merely asked 
sensation curvature was felt not, according the groups: (a) Easy, when 
sensation curvature was felt all; (b) comfortable, when sensation 
curvature was felt, but not thought uncomfortable; (c) slightly uncomfortable; 
and (d) uncomfortable. The effects the force and its rate change 
were separated the statistical analysis. 

was expected, there was considerable variation the assessment 
comfort, but the application statistical methods showed that comfort de- 
pends mainly although high rates tended reduce the tolerable values 
slightly. Speed does not appear affect the assessment comfort. 
uncomfortable noticeable value could detected low values 


CONCLUSIONS FROM THE EXPERIMENTS 


The main conclusions from the experiments are: 


(1) The driver does appear turn his wheel constant rate any 
given transition. 

(2) There such thing constant rate turning all 
transitions. The driver tends turn his wheel faster when entering sharp 
curve than when entering flat one. This conclusion accordance with 
general observation, but contradiction the previous theory. 

(3) There some marked tendency for drivers use about two thirds 
one half the curve the transition. 

(4) far comfort concerned, the important effect; and 
unimportant. 


From these conclusions appears that there sound reason for using 
the rate change force the criterion transition curve design, and 
suggested that this theory should abandoned. the other hand, may 
opinion, value which safe comfortable, and that its use convenient. 


counter this argument necessary consider the the 


compound curve. 


THe 


the condition for constant can satisfied 
curves other types than the usual spiral, this curve not unsatisfactory 
provided the superelevation properly designed comply with the conditions 
Eq. simple and easy use, and familiar, that there little 
gained abandoning it. the same time, the lemniscate has its ad- 
vocates, and, preferred, there reason whatever against its use. Never- 
theless, the spiral will used herein. Its mathematics have been dealt with 
elsewhere? some length and need not repeated. However, since has 
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been shown that the rate change force not important, also follows 
that the length the curve relatively unimportant, and may eliminated 
from the equations for the main dimensions the curve, basing these the 


Dimensions SYMMETRICAL CIRCULAR CURVE WITH 


radius. These modifications were first published slightly 
different form. 


Referring Fig. 


(8a) 


which the total tangent distance from the point intersection the 
starting point transition; the external distance from the intersection 
point the nearest point the curve; the radius the point the 
transition defined the angle and are the coordinates the same point 


for Circular Curve with Kursbatt, Highways, Bridges and Aero- 
London, February, 14, 1945, 226. 
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from the.starting point transition origin: 


TRANSITION CURVE 


Papers 


and 


these expressions, which apply whether not circular central arc 


Eas. 8d, 8e, 9a, 9b, 9c, AND 


(sec) 


Table these 
spiral functions six deci- 
mal places, for angles 25° 
half degree intervals. 
Corrections for the polar 
deflection angle Eq. 
are also given, although this 
usually negligible for small 
values 

The second form Eq. 
will found rather 
simpler for calculation pur- 
poses, spite its more 
complicated appearance. 

The Practical Problem.— 
Under practical road condi- 
tions the engineer will have 
two tangents, meeting 
angle which will have 
been measured the 
ground. there unlim- 
ited room available, com- 
plications are involved; the 
curve can made flat 
desired. Frequently, how- 
ever, some obstacle other 
restricts either the exterior 
distance (Fig. 2), the 
tangent distance The 
problem make the best 
use the room available. 
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Whichever these distances limited can substituted Eq. 
find provided the value known that the values the 
appropriate functions may determined from Table When ¢/A the 
curve will all-transitional; and, ¢/A decreases, increases, with corre- 


Values 


Values 


Qa je 
(a) Fixed 
(6) Fixed 


sponding reduction but the length the transition decreases, thus in- 
creasing Finally, when ¢/A the curve circular without transi- 
tions and minimum for the conditions, but will infinite. 

the superelevation assumed proportional the centrifugal ratio, 
quite usual practice, and the value any compound curve denoted 
and that for the corresponding all-transitional curve with the corre- 
sponding values denoted and the following expressions can 
written for the ratios between these quantities: For fixed 


and, for fixed 
For either condition, 
(11) 


The symbols and represent the values and for the value 
used, and and those for A/2. These expressions depend only the 
telative values the angles, and express the ratio the value any 
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compound curve which fixed that the all-transitional curve 
for the same conditions. Fig. gives curves for versus for 60°, al- 
though the effect small that they may considered approximately 
general. effect, these curves show the effects and caused the 

various values chosen the 
TABLE engineer, although the latter has been 

SUPERELEVATION expressed proportion transition. 
From evident that the 

Radius, higher proportions transition greatly 


between the all-transitional curve and 
0.045 one with 0.9 transition negligible, 
0.033 the ratio 1.00 1.03, although the 
justification whatever for the all-transi- 


elevation for 48, 0.025 per Thus, the case against fixed rate 


rate chosen, for fixed and fixed 

speed, long transitions and eventually 
all-transitional curves are demanded the room available decreases. 
decreases further, the ridiculous demand may made that the transition 
longer than the whole curve, even low values This argument usually 
countered ruling that such cases the must reduced— 
making nonsense the whole matter, because this theory completely ignores 
the driver, who cannot expected observe the “design speed,” is, 
course, completely unaware its existence. appears, fact, that the 
driver slows down reduce and not reduce and that, design speed 
that the rate change force should abandoned method transition 
curve design. 


ALTERNATIVE METHOD TRANSITION CuRVE DESIGN 


Although this paper only purports deal with general principles, some 
new procedure must suggested, since the effect the argument that the 
previous theory should abandoned. 

For the moment, the rule” first proposed favored. 
This rule states that curves should designed have about two thirds 
their length transitional, although the writer now inclines the belief that 


one half better proportion. The exact length the transition 


however, matter extreme importance, that may taken the nearest 
half degree. 

The proposed procedure illustrated the following example, which is, 
however, only intended demonstrate the method. The values used are not 
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put forward recommendations, nor the example intended show that the 
rule sound. The required information includes: 


The maximum speed for which the road intended, say (for purpose 
example only), miles per hr; 

The maximum value considered safe that speed, say, 0.15; 

The rates superelevation™ used, say, shown Table and 

The proportion transition, say, two thirds, that A/4, 
the nearest half degree. 


substituting the foregoing known values Eq. the minimum permissible 
radius curve will found, this example, 2,300 ft. 

Next, the deviation angle, and the greatest lengths and al- 
lowed the topography are measured the site. For the example, let 


22° 38’, and let the maximum possible ft. Then 


and, from Eq. 8b, using Table obtain the values the spiral functions, 

0.001535 0.019828 0.000030 2,570.93 ft. This value permis- 
sible, resulting 0.13, and superelevation rate 0.04 per ft. not 
permissible, some reconsideration the alinement necessary. 

From Eq. 8a, 2,570.93 (1.001535 0.200122 0.095964) 762.01 
ft, and 2,570.93 0.191986 493.58 ft. The spirals may set out 
any convenient manner, the ten-chord method being especially suitable. The 
method has the merit extreme simplicity. 


CoNCLUSION 


The whole argument this paper runs counter established practice. 
submitted, nevertheless, that the rejection the current theory based 
sound reasoning supported experimental analysis, and that the new con- 
cept far more accord with common sense and with everyday observation. 
The objection will doubt raised that takes the whole problem back 
rule-of-thumb days, and that vague and unscientific. However, does 
appear that vagueness inherent the subject, and that past theoretical and 
mathematical treatment not justified facts. 

Far from the length transition curve being theoretically determinable 
one the enthusiasts for the previous theory has claimed, 
varies widely the driver’s whim. Therefore, there seems little object 
going great lengths find it. The principal aim arrive some 
reasonable compromise. This, claimed, the proportion rule accomplishes. 
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ANALYSIS STEPPED-COLUMN MILL BENTS 


The technique applied the analysis stepped- 
column mill bents, this paper. Curves for obtaining the necessary distribu- 
tion coefficients and the fixed-end moments are included. The method permits 
evaluating the relative resistance, horizontal forces, columns having 


different stiffnesses due length, size, other factors. Several examples 


two-column and three-column bents are included describe the method fully. 

Although this method will require more time than that obtaining mo- 
ments from assumed point contraflexure, the additional time would seem 
justified, because the difficulty estimating the location such point 
contraflexure accurately for stepped columns. 


INTRODUCTION 


the design mill buildings, especially where the crane clearances require 
considerable height, the horizontal forces from wind and crane thrust are 
important factors. short buildings, simple enough transmit these 
forces the ground through bracing system the braced end bents the 
building. However, with the trend toward larger and larger industrial build- 
ings, the length increases that one more expansion joints are required; and 
horizontal forces must resisted each individual bent. Often, result, 
the bending stresses the columns become greater than the direct stresses, 
and the essentially becomes more vertical beam than column 
standpoint stress analysis. 

The simple and more familiar method analysis for such 
determine the bending moments from assumed point contraflexure. With 
column constant section from top bottom, the location this point can 
assumed with reasonable accuracy under any set assumptions regarding 
the base and top fixations. However, for column supporting crane runway, 
section” considerably more economical, most cases. Immedi- 


comments are invited for immediate publication; insure publication the last 
discussion should submitted March 1948. 
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ately, the question arises the location the point contraflexure; in- 
accuracy this assumption may cause considerable error the final result. 

hoped that this paper will present more satisfactory method for the 
analysis stepped-column mill bents. The loadings given the examples are 
used illustrations the method, and not Admittedly, 
these loadings are not the nearest approach the probable applied conditions, 
but they illustrate the various types design loads used designers. 

The method successive approximations, so-called distribu- 
tion,” used without explanation, determine the true end moments, since 
this method probably familiar most engineers and can found various 
textbooks and reference books. 

particular, the writer has used set unpublished mimeographed notes, 
Maugh, Assoc. ASCE, entitled Analysis Continuous 
Framed reference preparing this material. The discussion 
the basic theory given briefly herein Example For more detailed 
discussion the theory, reference can made “Statically Indeterminate 
also Professor Maugh. 


II. GENERAL AND ASSUMPTIONS 


With the increasing necessity for the greatest economy the use steel 
and the resultant use higher working stresses, has been more than ever 
important eliminate many uncertainties possible structural design. 
Since these columns are essentially beams, they 
should treated such. The analysis beam 
variable moment inertia generally under- 
stood; not difficult, although the work involved 
does require considerable time and somewhat 


good results, several industrial buildings erected 
during World War II. save time the analysis 
beams variable moment inertia, series 
curves presented which reduces the work that 
required for continuous beam two spans with 
the translation one joint. the well-known 
method successive approximations introduced 
Hardy Cross, ASCE,* the analysis not 
lengthy problem. 

with any design problem, certain assumptions 
must made simplify the problem and save time. Some these as- 
sumptions are known erroneous. However, the error involved 
small safely ignored. this method analysis the following as- 
sumptions were made (see Fig. 1): 


Hooke’s law applies. 
The principle-of superposition applies. 


Transactions, ASCE, Voi. 96, 1932, 
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the deflection bent, point and point move the same amount—. 
that is, point always directly above point (Unless the truss 
subject considerable deflection, this safe assumption.) 

The column simply supported points and a’—that is, the truss 
knee brace connections create effective fixity these individual 
points. Since the column usually quite rigid relation the truss 
knee brace members, this safe assumption. 

Both columns deflect the same amount—that is, any change length 
roof trusses caused horizontal forces too small significant. 

The relation between the moments inertia the upper and lower parts 
the column known. (This relation can approximated readily 
and quite closely after some experience; but the preliminary design 
indicates too large error the analysis can corrected quickly with 
new ratio moments.) 


(It can argued that, effect, the column fixed the top. points and 
are close together, true that there can little deflection between these 
points; but, this distance larger, the deflection may considerable. As- 
sumption more general and therefore safer than the assumption full fixity 
point a.) 

The explanation such design method shown more readily examples, 
and this device will used this paper. greater ease illustrating the 
principles the method, the following loading conditions were used the 


examples although realized that this not represent true 
wind action: 


(1) Beam bending column, from wind, occurs entirely the windward 
column; 

(2) Wind force assumed uniformly distributed over the entire length 
column; and 

(3) Horizontal force from crane applied entirely one column bent. 


Symbols and letter symbols this paper are defined where 
they first appear, the text illustration, and are assembled for conveni- 
ence reference Appendix 


III. EXAMPLE 


Example 1.—Consider column symmetrical bent similar those 
Fig. and fixed the base. the moment inertia the main section 
the column and that the smaller (reference) section above the step, 
let I/I, 8andz general slope-deflection equations, for moments 
beam (ab, Fig. with variable moment inertia, are: 


and 


| 7 
| 
318 
3 
ith 
ons 
as- 


1274 MILL BENTS Papers 


the angular rotation, with subscripts denoting point rotation; 
the relative displacement one end the beam the other end; the 
fixed-end moment; equals moment joint member ab; and equals 

joint member ab. the beam hinged point moment 
zero. Solving for Eq. and substituting Eq. 1b: 


Similarly, the beam hinged point 
El, 


The relation between the end rotations and the end moments can determined 
obtaining and 
the angles rota- 
tion joints and 
respectively, caused 
each case unit mo- 
and By, the angles 
rotation joints and 
respéctively, caused unit moment joint (Fig. For any end 
moments and then: 


Boa 
and 
Boa 
which 
and 


and fixed-end moments are temporarily ignored. Having 
determined the coefficients, the fixed-end moments can determined for any 
applied loads computing the angular rotations and the end tangents 
the beam considered simply supported. The fixed-end moment must pre- 
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vent any rotation the end tangents; therefore, 


and 

The angles Eqs. can readily determined the conjugate beam 


determining the reactions beam loaded with the 


(since constant). Figs. 3(a) and illustrate the use that method for 
this example. end moment ten units applied end and the length 


divided into ten units. For the load the reaction 


end moment ten units and length ten units, the end rotation end 


the same way, applying end moment ten units end found 


2 


0.05233 0.23329 


carry-over factor 0.22. Substituting these values Eqs. 


j 
(a) (c) (d) (e) 
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and 

shown Fig. the loadings this column which produce bending are 
three types—a uniform load 0.5 kips per ft, concentrated load (at the 
change section) kips, and eccentric load 120 kips. The upper load 
kips will also produce some bending, particularly the upper section; 
however, this will not especially severe, since near the center line the 

lower section. this example, 

clockwise resisting moments are 

designated “‘plus and counter- 

clockwise resisting moments, 

“minus the number within 

box indicates the ratio the mo- 

ment distributed that side 

joint; the number arrow be- 

tween joints indicates the 

over” the adjacent joint the 

direction the arrow; and, the 

completion single distribution joint, single line drawn under the 

numbers and after the final distribution (unbalanced factor negligible size) 
double line drawn under the numbers and summation made. 


a 


Symmetrical about Center Line 


Fig. 3(c) shows the for unit uniform load applied beam 


length before, the end reactions give the end rotations or, this case, 
angles and Substituting the values the end reactions Eqs. 


Fig. 3(d) shows the for unit load applied 0.3 from end thus: 


For unit moment applied 0.3 from end Fig. shows the 


Again substituting Eq. 


For 120 ft-kips Eq. yields 0.166 120 19.9 ft-kips 

For the column under consideration—with 5.63; 5.95; and 
5.95 


respectively, are 0.22; and 5.63 
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for the distribution factors point Fig. are shown Table and fixed- 
end moments for the several types load are shown Table The un- 
balanced moments end are distributed shown Fig. 5(a), assuming full 


Carry-Over Ratios 


1.06 


1.06 0.22 


Distribution Factors 
Fixed-End Moments 


End Moments 


UNBALANCED MOMENTS 
POINT 


(b) CONCENTRATED 
LOAD KIPS 


fixity the base; and the shear computed 97.0 34.8 


11.56 kips. 
like manner distributing the fixed-end moments (in ft-kips) for the 


and (in kips) 2.84. For the 


only 60% fixity assumed point 

5); and A, becomes 10.0 —— 3.00 Io/4 0.750 Io 0.84 


The horizontal forces, course, cause 
translation the top the column relation the base or, conversely 
(since assumption c—points and move the same amount), point 
moves some distance 
TABLE relation points and a’. 


causes additional end mo- 
Uniformly Distributed: 


moments are not readil 
follows: 
this case, assume units. From Eqs. the end moments for this 
translation, only, are: 


1.8 
m. J i 
8c) 
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and 


distribution, before, the final moments, foot-kips, become: 


and, kips: 
30.49 9.73 


1.01 


Values 
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For the uniform horizontal load the column, and with translation 
considered, was shown that the shear point the windward column 
was 11.56 kips. Since the total force assumed acting the column was 
kips force 11.56 10.44 kips force creating deflection the 
bent. Assuming that both columns resist this deflection equally, the force 
resisted each column 5.22 kips. Then the end moment for this 5.22 
ft-kips. 

For the 8-kip concentrated load the force causing bent deflection 2.84 
5.16 kips. Since the crane thrust spread adjacent bents the roof 


bracing, logical assume that three bents resist this force 1.72 


For design purposes, course, now only necessary total the moments 
from logical loading conditions which produce the most severe effect. 

stated previously, Example given only illustrate the method 
determining the moments column from any set loads under the condi- 
tions assumed, and not intended criterion loading conditions. For 
instance, most codes require percentage the wind applied acting outward 
the lee side the building. The method shown here would apply equally 
this condition. 

The general principles involved this example—such the determination 
constants and fixed-end moments for beams variable moments inertia, 
the distribution moments, and the computation moments from translation 
—are doubt familiar most engineers and are found various text- 
books and reference books. 


IV. Curves AND MoMENTS 


general, this method analysis relatively simple; but the deter- 
mination the constants and the fixed-end moments the construction 


diagrams somewhat laborious and time consuming; thus may have 


deterred others from such analysis. Figs. inclusive, illustrate series 
curves that will give the coefficients and fixed-end moments directly, thus 
reducing the work simple distribution described Example 
These curves are drawn groups five, each being for different location 
the change section and ranging from the top, tenths, from 0.1 
toz 0.5. 

The location the crane girder not likely below the midheight 
the column, curves were determined for lower points. The curves are all 
drawn with the varying values abscissas and the coefficients fixed- 
end moments the ordinates. Intermediate values being the dis- 
tance from the top column section the change section) can readily 
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0.16 


0.14 


0.12 


Values 


(b) MOMENT POINT 


Values 


For Force APPLIED THE CHANGE SECTION 


interpolated. two the series curves (Figs. 7(a) and 9(a)) was found 
Fig. reveals that 5.65, 5.95, and 26.5. Thus, the “carry- 
5.95 


over” from point point 1.05. Similarly, for the horizontal 


0.06 
0.04 
0.02 
S- 
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For the uniform load kips, from Fig. 0.0514 0.0514 


0.220 


(6) MOMENT POINT 


Values 


ft-kips. For the applied moment 120 ft-kips from Fig. 0.167 


7 | 4 
0.160 
0.140 
‘4 | 4 
0.060 
0.040 
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ft-kips. will noted that these values are all very close those obtained 


Example 


Thus, the use the graphs need constructing the several 


and determining the constants and the fixed-end moments from 


UNSYMMETRICAL BENT 


them. The time element thus reduced the point where this method 
analysis becomes entirely practicable and lends itself logical and accurate 


TABLE 


232.5 
12.9 76.1 


308.6 


165.75 497.2 
22.1 152.5 


649.7 


106.55 319.6 
88. 


For line 


result within the limits the original 
assumption. 

Finally, presents method ap- 
proach which columns adjacent 
bents opposite columns the bent 


are varying stiffness due mezzanine 


levels framing, heavier column section, 
other factors. 
APPLICATIONS SHOWING THE USE 
THE CURVES AND 
PROBLEMS 


three-column unsymmetrical bent. The 
solution this relatively simple prob- 
lem varies but little from the symmet- 
rical bent. Since the distribution the 
moments individual column has 
already been illustrated (Example 1), 
only the effect translation the upper 
joints the bent will shown. 

6.0 


24.0 


22.0 (from Fig. 6). Carry-over ratios are: 0.23 and 0.85. 


22.0 


ae 
wo ll N | 
0.75 
1.00 
0.235 
1.000 
0.78 
0.22 
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Also, 77.5; 12.9; and the distribution factors are 


and 
Carry-Over 0.85 1.04 0.71 
Ratios 
+356.1 


(c) SECTION 


These values fixed-end moment are distributed Fig. 11(a), with the 
132.6 43.9 


result that the shear point computed 7.35 kips. 
165.75; Kea 22.1; and the distribution factors are computed 


103.8 ft-kips 


and 


381.0 ft-kips 


for line these values fixed-end moment are distributed Fig. 11(6), 
356.1 79.4 


the resulting shear point being 14.53 kips. 
Carry-over ratios are: 0.22 and 0.71. Also, 


and the distribution factors are computed Table 3(c). 


q 
| 
4 
q 
or 
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50.2 ft-kips 


These values fixed-end moment are distributed Fig. 11(c), the shear 


From the foregoing then evident that column line takes 


7.35 
27.5% the thrust from bent sidesway; line takes 


14.53 4.89 


point being 


18.3%. 


4.—Fig. 12(a) shows simple bent with shed roof lean-to one 
side. The lean-to columns are relatively small and, assuming hinged bases, 
they provide horizontal stiffness the bent. However, since joint can 


From Roof 

12.4 Kips 
1.2 Kips 


move only horizontally (the column prevents vertical translation) the stiffness 
the lean-to roof truss insures that although joints and can deflect hori- 
zontally, they deflect the same distance. Thus, the column translation (E-line) 
this case, will combination two types shown Figs. 12(b) and 
12(c), and the total stiffness these columns much increased. 

somewhat different method required analyze this bent. First, as- 
suming that assign arbitrary value and distribute the 
moments produced. Then express the shears terms since 
distribute the moments, and express the shears terms The total shear 
the column between any two joints the sum these two shear terms. Thus 
knowing the force system which produces translation, equations can set 
for the solution and However, since column also offers resistance 
this force system, the method used Example must applied. 
horizontal forces acting only column the shears column produced 


| 
and 

6.0 Kips 

= e d | 

ac. 
ed. 

de. 
db. 


q 
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translation are the same between and and between points and 
Assume any value for this shear, and solve for and Then the total de- 
flection point the sum these 
for this particular shear. Since the tops 
ISTRIBUTION 

columns and move the same 

amount (assumption e), determine the 


shear the columns for the same 

A-value and the relative stiffness the 


ical example, first consider column 


Carry-Over Ratios 


Distribution Factors 
Fixed-End Moments 


430 430 


and the distribution factors are com- 
puted asin Table For 10: 


(a) 


For line: Assume and the same for column with 7.3. 


> 
| : 
0.71 
a 
+40.8 
1] 
e 
Ao, 
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20.0. Carry-over ratios are: 1.15; and 0.305. Finally: 


Distribution factors are Table 


These values fixed-end moment are distributed Fig. 


Carry-Over Ratios 0.5 


(a) ASSUMING 


Fixed-End Moments +73.4 
0.5 


Shear 51.6+29.7 =5.084 -29.7-6.8 


ASSUMING 
Fixed-End Moments 


End Moments 2.2 


lengths and each equal kip (forces points and total kip and 
reaction point equals kip) then: 


q 
and 
0.305 
4 a 
+28.0 
Reac 
bet 
the 
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and 
5.08 


Solving Eqs. and simultaneously: 0.913; and 0.124. Conse- 

Since both columns are assumed deflect the same amount the top, 
for column would 1.037 correspond the foregoing values. The total 
shear distributed the bent follows: 


Computation Shear 
Column 

5.08 0.124 0.63 
0.41 0.913 0.37 
Total shear bent.......... 1.10 


Column then resists 91% the horizontal thrust and column 


resists 9%. 


1.10 
Carry-Over Ratios 0.5 0.5 1.15 
0.305 


Distribution Factors 
Fixed-End Moments 


End Moments +16 


Simple Beam Shear 

Correction 


When the horizontal forces act columns and however, the shear 
between points and and between points and longer equal, and the 
relative values and change. The first step, previous examples, 
the distribution the moments from forces acting between joints, with all 
joints held against translation (see Fig. 15). 


| +3.3 106 | +10.6 -—7.1 +7.1 
9a) 
4 
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first step the solution the auxiliary force system, assume that all 
resistance translation taken the E-line columns. The force system 
the E-line columns then consists forces joints a’, and equal and 
opposite the reactions determined Fig. plus forces applied directly 
the joints such wind forces the roof. Since the base column 
hinged, half the uniform load 4.8 kips, shown acting it, assumed 
base and half the top pointe. This half point plus the kips acting 
the lean-to roof are distributed between points and column Fig. 12. 
Without another set calculations (translating points and distributing 
the moments produced, and determining the shears) the exact distribution 
between these two joints not known, but such distribution not needed for 
the present problem; only the shear needed between points and and 
between pointsaandc. Theshear between points 1.17 4.02 
kips plus 12.4 kips from the roof, total 16.42 kips. The shear between 


there relative translation between joints and the shear not needed. 
determine the relative values and set formulas similar 
except that the shear values are different: 


and 


Multiplying Eq. 10b becomes: 


solve Eqs. 10a and 10c simultaneously, and determine: 14.05; 

For A-value the B-line columns 18.35 units, the shear would 
0.96 1.76 kips. The E-line columns thus would take 
94.2% the thrust and the B-line columns only 5.8%. 


VI. 


Although most the examples this paper are not complete they show 
the general procedure for the various cases. with any method design for 
indeterminate structures, the ability make relative close initial estimate 
the members required very helpful and saves repeated computations. 
This ability, course, comes largely with experience and observation. 

Although the types problems described this paper means exhausts 
the variations that are likely occur practice, felt that practical 
approach other forms bents should not found difficult. 

The analysis probably not original, although not commonly 
used. However, without the use the curves presented herein, the solution 
would quite time consuming. Although this method solution with the 


27.63 1.76 


O 
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curves will require more time than that assuming point contraflexure and 
solving statically determinate structure, the greater assurance reason- 
ably accurate analysis should justify this extra time. 
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APPENDIX NOTATION 


The following letter symbols, used the paper and its discussion, con- 
form essentially American Standard Letter Symbols for Mechanics, Struc- 
tural Engineering and Testing Materials (ASA—Z10a—1932), prepared 
Committee the American Standards Association, with Society representation, 
and approved the Association 1932: 


nondimensional coefficient end moment the top stepped 
column (see Fig. 7); 

joint caused unit torque applied joint 

stepped column (see Fig. 7); 

coefficient, with suitable subscript, used slope-deflection equa- 
tions; thus: 

the distribution coefficient for point for given values 
I/I, and Fig. 6(a); 

the coefficient (see Fig. 6(b)) which, with C3, gives 
the carry-over ratio (the carry-over ratio from point 

the distribution coefficient for point the bottom the 
column, for given values and Fig. 6(c); 

moment inertia column section, the moment inertia for the 
reference section (the upper and smaller section stepped column) 
being designated 


Frames,” Evans, Edwards Bros., Ann Arbor, Mich., 1938. 


Mill Bents with Step Columns the Method Successive Approximations,” 
Daniel Ling, thesis presented the Univ. Michigan, Ann Arbor, 1943, partial fulfilment the 
Tequirements for the professional degree Civil Engineer. 
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stiffness ratio J/L used moment distribution, the stiffness ratio 
the upper and smaller section stepped column being desig- 
nated 

bending moment; the subscript denotes other subscripts 
refer designated joint and span—that is, general, equals 
fixed-end moment joint member ab; 

point load concentrated load; 

distribution factor; 

total weight, force, pressure; 

weight, force, pressure per linear foot; 

angular rotation the end tangents deformed beam, considered 
simply supported; 

substitution constant, function the unit rotation 

linear displacement one end beam with respect the other end; 

and 


for 


on 1 


— 
| — 
; 
| 
q 
5 
by 
p 
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DISCUSSIONS 


THEORY INELASTIC BENDING WITH 
REFERENCE LIMIT DESIGN 


Discussion 


MICHAEL HORNE 


Methods Stress Analysis Beyond the Elastic Limit.—Professor Hrennikoff’s 
method for dealing with inelastic beams ingenious, and has the advantage 
that, once the necessary tables have been compiled, can applied any 
members uniform section symmetrical about the plane bending. How- 
ever, suffers from the very serious disadvantage that for complicated problems 
the application Eqs. and becomes very involved. Many problems con- 
cerned with the elastic bending beams may solved quite simply means 
the three-moment theorem, and similar method may applied when the 
bending ceases elastic. effect, possible replace Eqs. and 
single relationship between the bending moments and deflections three 
successive points beam. 

Consider beam ABC (Fig. 30) which the bending moments points 
points and and between points and the bending moment diagram 
consists straight lines. the length, depth, and web area span 
hi, and respectively, and let the corresponding values for span 
he, and map, and are the values the function points 
and for the length AB, then, from Eq. 


Similarly, mgc and the values the function points and 
for the length BC, are 


paper Alexander Hrennikoff was published March, 1947, Proceedings. 
this paper has appeared Proceedings, follows: June, 1947, George Winter, and Drucker. 


Research Asst., Dept. Eng., Univ. Cambridge, Cambridge, England. 
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Let etc., the values the functions and corresponding 
these values Considering the center line the beam, let the de- 
flection point relative point and the deflection point relative 
point and are both measured positive downward. Referring 


Fig. 30, BC’ straight line tangential the center line point and, 
from the geometry the figure, 


Considering span AB, Fig. 30, readily proved that 


Similarly, for span BC, 


Substituting Eqs. 50a and Kq. 49: 


the beam ABC uniform throughout, he-and Hence, 


(a) 
dc 
tio 
are 
fra 
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there longer any need for the double suffixes, and Eq. 51a becomes 


Eqs. 51a and are valid for negative well positive values Ma, 
and long the following convention applied: The functions and 
are given the same sign the corresponding bending moment, and the function 
always regarded positive, regardless the sign the bending mo- 
ment. That this satisfactory readily proved the derivation Eqs. 
50a and 

Two simple examples will demonstrate the application Eq. 51. 

Example 1.—Consider the portal frame ABDE (Fig. uniform sec- 


Fie. 


tion throughout, supported pins points and Let its height 
and its span 21, and let carry vertical load point the midpoint 
span BD, and horizontal load applied beam level. Let Mc, and 
and let mg, mc, the corresponding values the function Moments 
are accounted positive when they produce tension the inner fibers the 
frame. follows from statical considerations that 


and 
the depth the section, and the web area, then, from Eq. 50b, 

and 


(d) 
1a) 
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For and can expressed terms the one unknown 
mp, and, hence 


and 


will assumed that all bending that results from axial loads may neg- 
lected, that Eq. 51b applicable. The vertical deflection point will 
denoted and the sidesway the frame, the entire 
frame continuous beam (Fig. 31(d)), Eq. may applied first the 
part ABC, then the part BCD, and finally the part CDE, thus obtaining 
three equations: 


and are eliminated from Eqs. 53e 53g, further relationship ob- 
tained between mc, and mp: 


obtaining solution for any given values and necessary as- 
sume value mp. the values and from 53c and 
and then test whether Eq. satisfied. When the correct values 
mc, and have been obtained, the values and may found from 
53e and 

Example 2.—In Fig. 32, consider beam AC, fixed points and and 
point the bending moments points and are Ma, and 


lo: 
th 
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respectively, follows from simple statics that 


(54a) 
ma, and are the corresponding values the function then, 


able four unknowns, mg, mc, and and necessary 
find two more. span imagined reflected about end into 
the position Fig. 32(c), then Eq. may applied the length 
AB, giving 


Similarly, span may reflected into the position and the application 
Eq. the length gives 


For any given value necessary solve Eqs. 51b, 54b, 55a, and 
trial and error until suitable values ma, mg, and are found. 

The Stress-Strain Relationship Steel Bending.—This subject has been 
quent discussions. stated Baker, Assoc. 
ASCE, and Roderick their discussion," 
the generally accepted theory England that 
the stress distribution similar that given the 
tensile test when made autographic testing ma- 
chine, care being taken insure that the applied 
load truly axial, but ignoring the effect strain 
hardening, shown Fig. 33. Under such circum- 
stances mild steel gives stress-strain relationship 
which the upper yield stress may between 25% and 45% excess the 
lower yield stress Professors Baker and Roderick also show that the 
assumption such stress distribution accounts very well for the behavior 
rectangular mild steel beams bending. Further work these writers, 
the Engineering Laboratory Cambridge University Cambridge, England, 
collaboration with the British Welding Research Association, has confirmed 
these results. The stress-strain relationship assumed Professor Hrennikoff 
differs from that Professors Baker and Roderick two respects: 
Professor Hrennikoff ignores the upper yield stress and takes strain hardening 
into account, assuming that strain hardening occurs bending precisely 
the same way tension test. Most the work done Professors Baker 


Strain 


Transactions, ASCE, Vol. 112, 1947, 1201. 
1238. 


| 
| . 
| 
| 
) | 
| J 
) 
7 
7 
7 
7 
q 
7 
iS- 
om, 


1298 HORNE INELASTIC BENDING Discussions 


and Roderick has been beams rectangular cross section, whereas, the 
paper Professor Hrennikoff refers also rolled beams which local yield 
due rolling stresses certain occur. Under these circumstances neglect 
the upper yield stress reasonable. However, the assumptions made with 
respect strain hardening cannot justified readily. The strains the 
end the plastic range are considerable, and under such circumstances 
should not assumed too readily that the strain-hardening effect the same 
bending tension. interesting note that the work done Cam- 
bridge indicates that the behavior both simply supported and continuous 
beams the loads causing large deflections may explained satisfactorily 
without reference strain hardening. 

The Importance Strain Hardening Limit Design.—Strain hardening, 
obviously, has important bearing attempts arrive design method 
based the ultimate capacity steel structures, but its importance should 
not overestimated. designing ordinary building frames the usual 
elastic methods, emphasis always placed stresses rather than deflec- 
tions. because within the elastic range the strains are small. 
when ultimate load method design used, the deflections become 
important the stresses, because the greatly increased strains within the 
plastic range. The structure will only useful that load which the 
deflections become excessive. The limit method design justifies itself be- 
cause the ultimate load calculated according the principle equalization 
moments must coincide approximately with the load which excessive 
deflections first occur. This because both the equalization moments and 

the appearance large deflections result from the same property mild steel— 
its capacity undergo plastic deformation without increase stress. The 
author has now shown that the plastic yield that occurs before strain hardening 
begins insufficient cause complete equalization moments. Since the 
deflections are becoming large, this stage, however, the ultimate capacity 
has been reached, that becomes pointless follow the behavior the 
structure further. The limit method design achieves its purpose leads 
reasonable estimate the load which deflections begin increase 
rapidly. Professor Hrennikoff’s work particular interest since affords 
means checking the assumption that this actually the case. The results 
obtainable this way may best illustrated taking numerical example. 

Example 3.—Consider portal frame the type shown Fig. 31(a), 
which Let the cross section the member uniform throughout, the 
web and flange areas being equal vertical load first 
applied that yield stress the extreme fibers just reached the center 
the beam, horizontal load only being applied thereafter. With steel having 
the properties assumed Professor Hrennikoff, the value just sufficient 


The entire frame remains elastic until has reached value given 


35.2 kips per in. when yield stress just reached the extreme fibers 
point Fig. 31(a). necessary obtain solutions the 


are 
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use Eqs. 53c 53h. The results are presented Table 9(a). Cols. 
Table 9(a), give the bending moments points and and the corre- 
sponding extreme fiber strains, and Cols. and give, respectively, the vertical 
deflection and the horizontal deflection the center the beam. The 


<Ips per € 
(1) (2) (3) (4) (5) (6) (7) (8) (9) 


(a) HARDENING CONSIDERED 


—26.40 44.00 —26.40 —0.66 1.10 —0.66 0.513 
35.20 44.00 —44.00 1.10 0.513 0.586 
40.00 —6.27 44.13 —46.27 —0.16 1.11 —1.44 0.515 0.677 
45.00 —3.22 44.68 —48.22 —0.08 1.18 —2.41 0.531 0.796 
50.00 0.56 45.96 —49.44 0.01 1.37 —11.00 0.568 0.970 
50.43 0.95 46.14 0.02 1.41 0.574 0.990 
51.03 1.51 46.40 —49.52 0.04 1.47 —18.06 0.583 1.022 

5.20 48.10 —18.49 0.660 1.236 
57.20 7.20 49.00 0.18 3.75 —18.80 0.727 1.379 
58.60 8.38 49.48 —50.22 0.21 18.00 0.816 1.517 

9.38 49.78 —50.62 0.23 18.46 —19.80 1.010 1.742 
65.00 12.94 50.84 —52.06 0.32 20.08 —22.00 1.597 3.824 
—26.40 44.00 —26.40 —0.66 1.10 —0.66 0.513 
50.43 0.95 46.14 —49.48 0.02 1.41 0.574 0.990 
51.03 1.53 46.42 0.04 1.47 0.584 1.022 
55.00 5.50 48.40 —49.50 0.14 2.61 0.679 
57.20 7.70 49.50 —49.50 0.19 0.840 1.515 


bending moment which strain hardening begins occurs first point ata 


value 50.43 kips per in., and subsequently point C-at value 


58.60 kips per For the sake comparison, results are given Table 


LEGEND 


ignored 


(Horizontal Load) Kips per 


based the assumption that the steel capable undergoing indefinite 
Strain without change stress after the attainment yield. Where data 
are given, the results are the same Table 9(a). This equivalent ig- 


Ss 
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noring the effects overstrain entirely, and complete failure occurs when equal 
moments sufficient produce yield stress over the entire section have been 


reached points and This occurs when 49.50 kips 


per in. and 57.20 kips per in. Fig. shows the bending mo- 
ments points and plotted against the horizontal load; and Fig. shows 


2 50 4 > H 


xX 30 bv 
B 


(Deflection) 


the horizontal load plotted against the horizontal and vertical deflections 
the center the beam. all cases the full lines refer the solution taking 
account the effect strain hardening, the dotted lines indicating what 
happens when this ignored. When strain hardening taken into account 
the deflections begin increase rapidly very soon after the value 
exceeds 57.2 kips per in. This represents very nearly the collapse load 
the frame. the collapse load obtained ignoring strain hardening 


very nearly the same that obtained when considered. The percentage 


differences the various calculated quantities 57.2 kips per in. 


hAw 
are follows: 
Description Difference (%) 

Bending moment at: 
1.02 

Deflection Point (Fig. 31): 


The percentage errors calculating. bending moments resulting from the 
neglect strain-hardening effects are negligible. The errors the calculated 
deflections are large but means unreasonably so. This interesting 
result view the considerable simplification which results when overstrain 
ignored. particularly the case highly redundant framework, 
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the analysis which would involve the simultaneous solution number 
formulas the type Eq. 51b. When overstrain taken into account, these 
must all solved trial and error, and, obviously, with more than three such 
equations solution would virtually impossible. overstrain ignored, 
however, and desired only derive the deflections the point col- 
lapse, the bending moments become statically determinate everywhere, and 
the equations may solved directly. 

The general complexity the behavior rigid structures makes difficult 
draw any conclusions the importance strain hardening limit de- 
sign. Three conclusions have been established, however, and are follows: 


(1) Tests the Engineering Laboratory Cambridge University indicate 
that the behavior steel structures the point collapse may satisfac- 
torily accounted for without any reference strain hardening. 

(2) Theoretical analysis indicates that, the collapse load calculated 
according the limit design method,.the error the deduced values the 
bending moments very small. The deflections deduced are somewhat 
excess the values obtained when strain hardening considered. 

(3) The neglect strain hardening simplifies, very considerably, the stress 
analysis structures beyond the elastic limit, particularly the theoretical 
collapse load approached. 


must emphasized, however, that these conclusions are tentative. 
particular, Professor Hrennikoff indicates, the presence semirigid joints 
may much prevent the equalization moments, and, under such circum- 
stances, the behavior members the strain-hardening range may 
considerable importance. 


18 
in. 
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DISCUSSIONS 


STABILITY THIN CYLINDRICAL SHELLS 


Discussion 


BATDORF AND MANUEL STEIN 


Esq., and MANUEL STEIN,” interesting addition 
the literature buckling cylinders torsion presented Professor 
Sturm. far the writers are aware, one else has given solution the 
problem determination the critical torque cylinders finite length 
the differential equation approach, which equations other than the simplified 
ones are employed. The results for very long cylinders that 
buckle into elliptical cross section (two circumferential waves) are all 
probability more accurate than those Mr. Donnell, since the approximations 
the Donnell differential equation were based the assumption that the 
square the number circumferential waves large compared This 
discussion primarily for the purpose calling attention the situation 
the limiting case very short cylinders, and discuss certain features the 
derivation used. 

Physically, the radius the cylinder under torsion approaches infinity 
while the length and thickness remain fixed, there obtained infinitely long 
flat plate width subjected shear. this case, the critical stress 
given 


which This agrees with the form used Professor Sturm 


paper Sturm was published April, 1947, Proceedings. 

Physicist, Langley Memorial Aeronautical Laboratory, National Advisory Committee for Aero- 
nautics, Langley Field, Va. 

Aeronautical Engr., Langley Memorial Aeronautical Laboratory, National Advisory Committee for 
Aeronautics, Langley Field, 

Thin-Walled Tubes Under Torsion,” Donnell, Technical Report No. 479, 

Advisory Committee for Aeronautics, Washington, C., 1933. 

Shear Stress Infinitely Long Flat Plate with Equal Elastic Restraints Against Rotation 
Along the Parallel Elbridge Stowell, Advance Restricted Report No. 3K12 (Wartime Report 
National Advisory Committee for Aeronautics, Washington, C., 1943 

the Stability Under Shearing Forces Flat Elastic Strip,” Southwell and Silvia 
Skan, Proceedings, Royal Soc. London, Series May 1924, pp. 582-607. 
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(Eq. 


Since the critical stress curved sheet can never below that the corre- 
sponding flat sheet, Fig. unduly conservative when predicts stresses 
this value—that is, curve should pass the left this straight line. That 
each the curves Fig. should faired into this line, rather than join 
abruptly, apparent from other studies this 

The problem buckling finite cylinder torsion defied solution for 
many years because the difficulty satisfying simultaneously the differ- 
ential equation and the boundary conditions. Professor Sturm’s solution does 
indeed satisfy the differential equations, but not the condition simple sup- 
port. That satisfies the differential equation and one the simple support 


2. 
not, however, satisfy the other simple support condition, 
can seen differentiation the buckle form assumed (Eq. 31a). 
Some interesting light thrown 
this deflection function and the TABLE RESULTS 
Sturm method solution solv- DIFFERENT 
all ing the simple problem 
finitely long flat plate shear 
this method, and comparing the Wave Wave 
his results with those obtained the Computation length, angle, 
energy method using the same de- (2) (4) 
results obtained exact methods 5.7 1.224 0.707 
Sturm method............ 6.7 0.95 1.1 
ity for simply supported and edges 0.84 
edges. These results are sum- 

Nodal angles are not given for the exact solutions because these solutions 
give curved rather than straight appears that Professor Sturm’s 
35) method exact solution for type restraint intermediate between simple 

support and clamped edges, but which the moment applied the edge out 
phase with the slope quarter wave length. 
For other than very short cylinders the difference buckling load for 
clamped and simply supported edges quite small, that not anticipated 
for that any considerable error critical stresses will result from the peculiar 
Stress Thin-Walled Cylinders Torsion,” Batdorf, Manuel Stein, and Murry 
ation Technical Note No. 1344, National Advisory Committee for Aeronautics, Washington, 
eport 


361. 
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boundary conditions implied. However, from the comparison Table 
appears that there may considerable error the buckling configuration. 
ingenious suggestion made Professor Sturm how use results 
for complete cylinders obtain approximate values for the shear buckling 
stress the curved panels between adjacent stiffeners stiffened cylinders. 
may interest that calculations have been recently published for the 
critical shear stresses simply supported and clamped curved 


Shear Stress Curved Rectangular Panels,” Batdorf, Manuel Stein, and Murry 
Schildcrout, Technical Note No. 1348, National Advisory Committee for Aeronautics, Washington, 
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DISCUSSIONS 


FORECASTING PRODUCTIVITY 
IRRIGABLE LANDS 


Discussion 


LEwIs 


ASCE.—A very interesting method forecasting the 
productivity area has been presented. However, based climatic 
conditions; and, because climate only one number important factors 
crop production, obvious that any forecast which does not include other 
factors, particularly soil characteristics, must subject wide variation 
when applied single project. 

The author suggests that the crop census data secured the United 
States Bureau Reclamation would furnish basis for very interesting 
extension his study. the limited time available the writer has attempted 
this for the projects Region the Bureau, which essentially the 
Columbia River drainage area. Figs. 2(a), 2(b), and 2(c) show the relation 
between the yield alfalfa, potatoes, and sugar beets, respectively, the 
total heat units those projects the region where sufficient acreage the 
given crops grown make the results significant. each figure has been 
indicated, addition, the comparable data for two Bureau projects outside 
Region These two projects have been chosen indicate the wide varia- 
tion which found when projects scattered over wide area are used for this 
comparison. The yield data are the 5-year averages, 1941 1945. 

will noted that the sugar beet and potato yields for Region projects 
conform linear relationship. Alfalfa yields not follow this relationship 
closely. general proposition, seems doubtful any one crop can 
used indicate the productivity over the extremely wide range suggested 
the author. well known practice, alfalfa such areas the 
Imperial and Salt River valleys grows very little during the extremely hot 
months. Potatoes have such short growing season that they can use only 
part the total heat most sections the West. possible grow 


two crops potatoes, one after the other, single season areas where 
paper late Muldrow was published February, 1947, Proceedings. 
this paper has appeared Proceedings, follows: 1947, Arthur Johnson, 
Van Loo, and Will and June, 1947, Charles Kirby Fox. 
With Region Bureau Reclamation, Boise, Idaho. 
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the total heat units amount 10,000 day degrees more. This not usual 
practice and therefore impossible tell what the total production 
potatoes might the whole the long growing season were used the 
sections; nor sugar beets require all the growing season the more 


Tons per Acre 


250 


Potatoes, Bushels per 


Sugar Tons per 


Total Heat Units,in Thousand Day-degree above 32° 


southern valleys. It, therefore, appears that crop available which could 
used show the productivity over very wide range—although seems 
probable that the hotter climates sugar cane might used for this purpose. 
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FRICTION COEFFICIENTS LARGE TUNNEL 


Discussion 


TAPLEY 


their success determining values roughness coefficients under complicated 
conditions for the various types surface the structure examined; the values 
coefficients obtained confirm, for very large-sized structures, the existing 
ideas the magnitude such coefficients. The application Eq. 
determine the measure wall roughness, shown Table 


Values 


ild 


Authors’ data (Table 13) 
Extrapolated (Fig 17) 


values 


particular interest are the values bend-loss coefficient, Cz, Table 
deriving these values, certain assumptions were made the preliminary 
calculations. Perhaps the effect these assumptions mirrored some ex- 
tent the final results; for example, the authors’ application Fig. 


paper Hickox, Peterka, and Elder was published April, 1947, 
Proceedings. Discussion this paper has appeared Proceedings, follows: June, 1947, Weston 
Gavett; and September, 1947, Barrows, Hunter Rouse, Karl Kennison, Chapman, 
Julian Hinds, and William Creager and Stephen Haybrook. 


Associate Civ. Engr., Bureau Eng., Los Angeles, Calif. 
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Creager and Members, ASCE, tacitly assumes that the minimum 
value corresponds r/D and their application Fig. from the 
same reference may predetermine the manner variation with the 
central angle. far the authors’ final results bend loss may con- 
sidered independent the assumptions, they furnish basis for developing 
data cover considerable range design conditions shown hereinafter. 

Data, developed, can applied design 

practice with less risk error, perhaps, than 
can existing data. 
With reference Fig. 16, may noted 
that points representing bends and fall 
along straight line representing r/D 13.9. 
For r/D only one point (bend was 
available. Bend however, has value 

Extrapolated (Fig 16) r/D 5.63, and extrapolating short dis- 

tance Fig. 17, another point representing 

D/r furnishing two points through which the locus 

representing r/D was drawn. Here, 1/A, 

course, equals r/L, where the length curve and Aisinradians. Calcu- 
lations are shown Table 13. 

passing curves through the origin and points representing the data 
Fig. 16, Fig. was constructed with maximum points D/r 0.2 (r/D 5). 
Passing the curves through the origin Fig. equivalent the assumption 
that infinite when constant and infinite while finite; also, 
then infinite. 


Values 
ul 


A=L/r 

angle, 

(degrees) (radians) 


Verlust mit gleichbleibendem Kreisquerschnitt,” Albert Hofmann, 
Mitteilungen des Hydraulischen, Instituts der Technischen Hochschule, Heft 1929, Fig. 


the other hand, may noted Fig. that, when very great, 
although r/D finite, not infinite; fact, the intercepts the axis 
ordinates are 0.5 and 1.5, from which and 0.667, 
respectively (for r/D 13.9 and respectively). this case has maxi- 


William Creager and Joel Justin, John Wiley Sons, 
New York, Y., 1927, pp. 119-120. 


| 
r/ d 
th: 
the 
are 
Bend 
(1) (2) (3) (4) (7) (8) 
34.26 0.236 0.598 1.67 0.200 4.24 
15.72 0.207 0.274 3.65 0.0719 4.83 the 
18.32 0.239 0.320 3.13 0.0719 4.18 
12.68 0.114 0.221 4.52 0.178 8.78 (an 
4.70 0.066 0.082 12.2 0.0719 15.2 
8.72 0.119 0.152 6.58 0.0719 8.41 wht 
90.0 0.40 1.57 0.637 2.50 
90.0 0.86 1.57 0.637 1.16 
mai 
beli 
alre 
mer 
4 


October, 1947 TAPLEY TUNNEL COEFFICIENTS 1309 


mum value depending the value r/D. The implication that for very 
long curve very long radius (but finite pipe diameter) infinite; whereas, 
for very long curve finite radius and finite pipe diameter (such coil 
many loops) the coefficient has finite value depending onr/D. This seems 
suggest that the velocity distribution becomes adjusted the curvature 
the latter case, but continues disturbed the former. 

Since the authors’ data covered central angles only little less than 
35°, was desired obtain confirmation the data Fig. the region 
larger central angles. Hofmann’s the Williams, Hubbell, 
and Fenkell data referring 30-in.-diameter pipe curve with 90° central 
angle 0.637) was used the writer calculate values for 
r/D 5andr/D The calculated values plotted exactly the loci 
Fig. previously derived from the authors’ data. was found, however, 
that other experimental data such that Balch, Davis, Brightmore, and 
others, the authors’ data considerably. 

Pending further confirmation the data Fig. 16, might worth 
noting that the loci may represented 


which and are functions r/D, and represent, respectively, the inter- 
cepts the axis ordinates and the slopes the loci. 

The straight-line equation (Eq. 12) can also represented the locus 
equilateral hyperbola whose coordinates are and and which has 


0.28 


Only that section 
0.24 


the locus for which coordinates 
are positive is, course, real. 
The data are plotted hyper- 
bolic form Fig. 18. That the 
curves pass through the origin 
Fig. follows from the fact that 
the head loss due the bend 
(and, therefore, Cg) must zero 0.04 Authors’ 

when zero. 

Although the authors’ de- 
velopment data loss head 
bends was auxiliary the 
main purpose their paper, 
believed that the data have some merit basis for design material, 
form such Fig. 16, which additional lines can added the process 
already described. Needless say, further confirmation other experi- 
mental data indicated. 


Verlust mit gleichbleibendem Kreisquerschnitt,” Albert Hofmann, 
Mitteilungen des Hydraulischen, Instituts der Technischen Hochschule, Heft 1929. 
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DISCUSSIONS 


CLASSIFICATION AND IDENTIFICATION 
SOILS 


Discussion 


JAMES STRATTON, AND DONALD BELCHER 


James ASCE.—The very nature military operations 
active theaters war precludes other than “rule methods the 
design airfields and roads, and the fluidity the military situation generally 
requires that their construction and improvement undertaken stages. 
Improvements bearing capacities must generally accomplished the face 
enemy activity and increasing severity use the engineer’s own forces. 
The military engineer must know, all times, the capacities his roads and 
airfields terms both traffic volume and loading intensities. This calls for 
rapid interpretations analyses, and therefore some standardization the 
“rule design method required insure against faulty designs 
interpretations the strengths roads and airfields which would adversely 
affect the plans the military commander. 

the spring 1942, the Corps Engineers, United States Army, adopted 
the California highway method evaluating the bearing values road 
subgrades and extended the design curves developed Porter, 
ASCE the higher loadings later required for airfields (see 
attempt was made correlate the Public Roads Administration soils classifi- 
cation with the bearing-value curves order provide design method for mili- 
tary road and airfield construction the zone operations through simple soil 
classification tests. result the inability the Corps Engineers satis- 
factorily correlate the Public Roads Administration soils classification with 
design procedures, Professor Casagrande was requested study methods 
soils classification. material for the paper was taken from the resulting 
study. 

Prior adopting the Casagrande classification, canvass was made all 
district and division engineers with the interesting result that, the forty-four 


Arthur Casagrande appeared June, 1947, Proceedings. Discussion this 


paper has appeared Proceedings, follows: September, 1947, Ralph Fadum. 
4Col., Army, Superv. Engr., The Panama Canal, Diablo Heights, Canal Zone. 


Numerals parentheses, thus: (25), refer corresponding items the Bibliography (see Appendix 
the paper), and the end discussion this issue. 
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offices canvassed, thirty nine were favor the adoption and only five raised 
objections. The objections were largely concerned with the confusion that 
would result the adoption revised soils classification when construction 
was peak. the light the favoring majority opinion, ap- 


LEGEND FOR GROUP SYMBOLS 


Clay 

Fines (material less than 0.1mm) 
Gravel 

High Compressibility 

Low Medium Compressibility 

Very Fine Sand, Silt, Rock Flower 
Organic 

Poorly Graded 

Peat 


Thickness, inches 


NO 


= 
= 

= 
oO 
3 
2 


PRA 


Flexible 


California Bearing Ratio 0.1 in.-Penetration 


RECOMMENDATIONS FOR COMBINED THICKNESS PAVEMENT 


peared that change that time would actually facilitate the transition giv- 
ing wide opportunity for correlation with the Public Roads Administration 
system classification that would aid later revisions and improvements 


the Casagrande classification. Accordingly, the Casagrande classification was 
adopted standard for the Corps Engineers. 
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Objections the soil classification method were interposed some the 
most experienced soil technicians, who condemn all soils classification methods 
the ground that their use encourages casual design conclusions whereas ex- 
haustive investigations and study are required. weighing this objection 
must noted that identified with soils classification are form 
which will always popular the profession despite 
the concern those who recognize the complexity soils and the difficulties 
predicting their behavior. The Corps Engineers does not propose sub- 
stitute soil classification method for needed investigations and tests and or- 
derly design procedures except where circumstances preclude other than rough 
rudimentary methods design. 

was the conclusion the Corps Engineers that the the 
soil classification methods proposed Professor Casagrande with the bearing 
values subgrades and base courses would helpful developing 
thumb” method design that could employed where extensive investiga- 
tions, tests, and deliberate design methods were infeasible. 

The virtues ascribed districts and divisions the United States Engineer 
Department (U.S.E.D.) the Casagrande method classification may 
summarized stating that the improved classification consonance with 
general progress soils engineering. One division the U.S.E.D. commented 
favorably the method because its mnemonic characteristics. 

The Lower Mississippi Valley Division, U.S.E.D., noted that soil deriv- 
ative loess, frequently encountered the Mississippi Valley, consisted 
from 40% 60% silt and partakes the nature both A-4 and A-5 
soil, and thus not classifiable the Public Roads Administration method. 
This division found that this soil also was not classifiable the Casagrande 
method and that further study with view eataloging this soil required. 
The same division also pointed out that soils classification correlation with the 

California Bearing Ratio (CBR) values must predicated standard 
methods soil compaction and testing for the bearing values. The United 
States Waterways Experiment Station Vicksburg, Miss., investigated wide 
variety soils and reported that the bearing values very well with 
its interpretation the Casagrande classification with one exception, which 
fell outside the estimated range CBR values specified Professor 
Casagrande for the particular soil type. 

One district suggested that the percentage soil binder and its plasticity 
index indicated numerical values the classification gravels and sands; 


for example, GW-3,1 would indicate well-graded gravel with approximately 


30% soil binder and plasticity index approximately 10%; likewise, SP-0 
would indicate poorly graded sand with fines. this plan were adopted, 
the GC, GF, SC, and groups would eliminated. For the fine-graded soils, 
this district proposed correlating the plasticity index and moisture condition 
with the classification letters. For example, C5-3 would represent ‘‘fat” clay 
with plasticity index 50% and moisture content 30%. There are 
tain virtues this suggestion which may realized indirectly adopting the 
suggestion the Upper Mississippi Valley Division, U.S.E.D., that the CBR 
number added after the letter symbols; thus, for example, SF-25 would in- 
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dicate soil sand with excess fines and would give picture its 
internal structure. This suggestion value since the CBR indicative 
behavior under load. probable that this form identification will used 
when bearing tests have been made. 

The North Atlantic Division, U.S.E.D., stated that the suggested soil groups 
are satisfactory with the exception the gravels and sands the fine groups— 
that is, and SF, respectively. These could subdivided advantageously 
into additional groups reducing the wide latitude now covered them. was 
suggested this division that additional gravel and sand group, designated 
and SL, respectively, added the Casagrande classification, include 
gravels and sands with excess fines. 

The Casagrande classification not expected local names en- 
tirely for special soils, such loess, buckshot, gumbo, These names have 


special meaning the localities occurrence and expected that they will 


used supplement the Casagrande classification. 

The foregoing discussion was prepared for delivery the Annual Meeting 
the Society January, 1943. Subsequently, the discusser’s duties deprived 
him opportunity participate engineering activities where might have 
noted the results use the Casagrande soils classification system. would 
interest learn the method classification found favor and as- 


sisted military engineers the construction roads and airfields the theater 
operations. 


the present status soil classification and identification provides well-organ- 
ized framework that necessary perspective the current situation 
this field. The trial and error process, plus native ingenuity, has produced 
“system” for each large organization dealing with soil the engineering sense. 
This-review organizes soil classifications and identification and places them 
juxtaposition that they may evaluated and their relationships noted. 
too often various methods are treated with intolerance that indicates lack 
scientific maturity perspective. Much this intolerance appears when 
something is, the system wrong. The same people throw golf 
clubs and tennis rackets because their game poor. 

Perhaps there have been many instances where each system has been applied 
too broadly—used beyond the limitations originally perscribed. Professor Casa- 
grande points this out with respect the Public Roads (PR) system. true 
other systems well and cannot considered fundamental fault. 
this respect highly desirable development Fig. would introduce 
vertical distribution based grain size. The author indicates this particularly 
under the heading, “3. Textural Soil Classifications,” where noted that 
such method inadequate for fine-grained soils. Similar limitations can 
placed several the others, particularly those Fig. under the division 
dealing with ‘‘Properties Soil its Undisturbed Condition.” 

the detailed discussions the various classifications, especially those 
under the heading, Agricultural Soil Classification,” and Geologic Soil 
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Classification,” some additional attention will add the perspective that- can 
gained from the over-all treatment the subject. Primarily, erroneous 
describe the first method Soil Classification.” Engineers 
are fact unintentionally perpetuating careless description science; this 
the proper place correct this concept the minds civil engineers. 
tell engineer that this agricultural method (as has been done for long) 
send him elsewhere looking for engineering method. Unfortunately, 
describe pedological soil classification dubious improvement. is, 
more simply, the science soil formation; the classification based all the 
properties the soils. Therein lies advantage not found other methods. 

First, and for the record stands the original paper, would more 
complete state that classification for the purpose describing soil 
and the conditions under which that soil exists. that basis there funda- 
mental leaning toward any one field application, but any field interest may 
interpret the data for its own particular use. there technical 
toward agriculture the part some pedologists, attributable the fact 
that agriculture supports their programs the allotment funds. 

regretted that influential paper such this should label the 
pedological system complicated one. the contrary, its principles are 
simple and logical, particularly when the method interpreted for engineering 
use. also possesses simplicity another order: With this system soil 
classified three dimensions, length, width, and depth. With one well chosen 
auger boring soil area can classified, whereas other methods the area 
shrinks the size the auger hole and sample has little significance beyond 
the hole from which was obtained. 

has been commonplace find the pedological classification described 
the literature classification surface soils for agricultural purposes. 
though Professor Casagrande described this way, also states that the 
underlying parent materials are included. The term soils” relative 
one. When used conjunction with “agricultural purpases,” “surface 
soil” immediately interpreted meaning most in. depth below 
the surface. For the sake clarification this point should examined some 
detail and emphasis added Professor Casagrande’s inclusion the parent 
material, which may extend considerable depths. Actually, there arbi- 
trary boundary basement beyond which lies the realm the geologist, and 
there particular need for such limitation. 

one will advocate this method for material existing depths common 
deep foundations. Therefore, depths from 100 will never 
value. Depths from are still much too great for any general 
use this method, but sand dunes that depth are found areas that 
are can classified the pedological basis. Also, there are extensive areas 
deep wind-blown silt that reach these and even greater depths (from 100 
200 ft). These may called Aeolian (geologic classification) loess deposits but 
they also coincide-with pedologic units. Along the Lower Mississippi Valley 
these are classified the Memphis soil; the north, especially Iowa, Illinois, 
and Missouri, they are described Clinton soils; and the State Washington 
Palouse. 


att 
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the surface approached, the 40-ft 60-ft depth range, the application 
pedology becomes somewhat greater although only geologic deposits that form 
the parent materials are yet under consideration. this range are included 
much rock and drift the glaciated areas. those soils weathered from rock 
there are valid relationships many areas between the soil mantle and the char- 
acter the rock such depths. the drift areas the character the mass 
material below the surface sufficiently well described, that general com- 
paction and several other requirements (26) can anticipated the basis 
the series name soil unit pedological classification). 

Each 10-ft increment, the surface approached, emphasizes two trends: 
First, that the application the pedological principles increases its greatest 
value the 10-ft zone immediately below the surface; and, secondly, that the 
dollar volume highway and airport construction increases its greatest vol- 
ume within the same general zone. 

This somewhat brief illustration intended show that, national basis, 
sliding scale applicability necessary. smaller units, such within 
state district, the scale becomes more fixed. will found that the signifi- 
cance the pedological method, applies soils the depth, function 
the type deposit formation. areas where bare rock exposed, or, where 
very shallow soils exist rough terrain, the value the method vanishes. 

Climate also important. Within the United States climatic differences are 


sufficiently great that they strongly influence soil development. the humid 


regions much water percolates downward through the soil, constantly leaching 
the soil and rock, and general building accumulations clay depths 
several feet below the surface. this process soluble materials are removed. 
the west, however, the situation often reversed. There, rainfall and melting 
snows the mountains are contrast the desert climate prevailing the in- 
tervening valleys few thousand feet below. Much the water proceeds 
ground water the valley area, where largely dissipated evaporation 
from the surface. This reversed process precipitates salts near the surface 
form the alkali and caliche soils. are also areas where the clays may easily 
become saturated with exchangeable ions and thus take characteristics 
abnormally high volume change and plasticity. Such soils occur far east 
western Minnesota, where small areas alkaline soil are associated with the 
Fargo (clay) series. 

present relatively little correlation has been accomplished this way, 
but known that pavement performance, base thickness requirements, and 
many construction details are related the individual soil series. The alkali 
attack concrete slabs, head walls, culverts, and steel pipe severe action 
many these areas (27). 

this method adaptable engineering use, then who using it, where 
being used, and how? examination this field reveals the following: 

The Portland Cement Association its work soil-cement stabilization uses 


and recommends the use the pedological method basis for obtaining soil 
data (28). 


[(26)] uses the pedological system used the Bureau 
Chemistry and Soils, soil surveys for highway purposes. Mis- 
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souri takes attitude similar that Michigan, that, the three best 
known methods classifying soils, the pedological system seems best suit 
their needs. The other two methods seem limited their scope 
the classification texture, the triaxial chart takes into account 
other soil properties and the classification developed the Public Roads 
Administration more index the behavior disturbed soils and 
does not take into account field conditions. The pedological system es- 
sentially reflection drainage, topography, parent material, weathering 
agencies, well textural class. Missouri, however, does require that 
samples each horizon each soil type taken and submitted the test- 
ing laboratory for analysis. Field density tests are also made each horizon 
order that the theoretical soil shrinkage can calculated for each horizon. 


“The Michigan State Highway Department [(26)] has found that the 
pedological scheme soil classification provides excellent means study- 
ing soils for general highway purposes, and that serves foundation for all 
other soil studies the department. The soil type, addition serving 
unit mapping, also serves unit which laboratory test results, road 
behavior observations, construction problems, and maintenance experiences 
are related. other words, the soil type serves means for classifying 

and coordinating the accumulated soil experience the entire highway de- 
partment 


The Nebraska Department Roads and Irrigation (29) states: 


“The pedological system soil classification which has international 
recognition, was selected. 

“By adopting this method classifying and indentifying soils, was 
possible make use wealth available scientific information soils, 
including the published county soil maps and reports, and supplement 
such knowledge with the test data and information obtained observation 
this survey well that which may obtained through future investi- 
gations.” 


Tennessee and Nebraska have established (30)(31) correlation between 
the soil series name and the texture the material, further correla- 
tion between these and pavement performance. 

The writer has observed distinct positive correlation between the pumping 
rigid pavements and the occurrence the Superior clay series near Superior, 
Wis.; the Sharkey clay series near New Orleans, La.; and the Susquehanna clay 
series Texas (partially reported (32)). 

There fair but incomplete correlation between the California Bearing 
Ratio (CBR) value and the soil series. This indicated the results United 
States Engineer Department tests the Clermont and Brookston soil series 
Indiana; and inference the reported CBR experiments Stockton, Calif. 
(33). 

The North Carolina Department Highways reports (32) was 
not encountered soils the Norfolk Ruston series. Pumping occurred 
the Portsmouth and Plummer series *.” occurred also the Lufkin, 
Elkton, and Coxville series. 

1935, Janda (34) reported that excessive longitudinal cracks 
particular concrete pavement were confined area Colby silt loam. 
Housel, ASCE (35), and Reagel (36) have related soil-cement con- 
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struction the soil series and Hicks (37) has established the relationship 
between the principal horizons several soil and the cement required for 
proper stabilization. Public Roads Administration (PRA) studies (38) show 
results typical soils the coastal plain (Norfolk, Portsmouth, and Ruston) 
and their relationship base course performance North Carolina. 

Professor Casagrande’s description “5. Geologic Soil 
significant. would considerable value the method could evaluated 
greater detail, perhaps less subtly; for truly “it damned faint 

The weakness the method seems lodged the broad basis the 
classification and the recognition age and general origin the most im- 
portant factors. For those reasons can seldom used engineering without 
considerable modification. evaluation this situation was concisely ex- 


pressed Karl Terzaghi, ASCE (39), summing the results investi- 
gation this field: 


are definite relations [between geologic classifications and engi- 
neering properties], because for given geologic origin and petrographic char- 
acter stratum the strength the material, rock, soil depends con- 
ditions that are not, and cannot be, recognized purely geological inves- 
tigation. For almost every geologically well-defined stratum the possible 
its physical properties from the statistical average impor- 
tant, that one cannot even make acceptable guess based 


Essentially the pedological classification includes refined geological classifi- 
cation and the parent material soil directly related materials that are 
mapped and described geologists. Olmstead (40) and others have 
shown this various sections the United States. 
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DISCUSSIONS 


MISSISSIPPI RIVER BRIDGE DUBUQUE, IOWA 


Discussion 


not quite clear why the builders this bridge 
made the bottom chord the end spans straight and level with the suspended 


floor the arched center span. Additional navigation clearance could 


gained, without increase cost, making the top chord the end spans 
straight, while arching their bottom chord. This arrangement permits the 
lowering the midpiers, the suspended floor the center span also raised 
the level the top chord the end spans. perusal the literature the 
subject shows that designers this type steel bridges prefer this arrangement. 

the stress analysis the truss the present bridge, the authors use the 
customary procedure removing first the quantities not obtainable statics. 
subsequently applying the principle virtual displacements the base 
system thus formed, the necessary elastic equations are established. However, 
the choice this base system and the method obtaining the elastic equations 
are open criticism. The base system cumbersome and the derivation and 
solution the elastic equations laborious. The coefficients Eqs. 17a, 
and 17c are such that for their solution great accuracy needed obtain useful 
results. 

Substantial savings labor could have been achieved, the authors had 
used more advanced methods stress analysis. These methods, although 
long have found, with one little favor among authors 
and stress analysts the United States. 

always possible select the hyperstatic quantities, their direction, 
and point application such way obtain many independent elastic 


paper Bergendoff and Josef Sorkin was published June, 1947, Proceedings. 
Discussion this paper has appeared Proceedings, follows: September, 1947, Jonathan Jones. 

Asst. Structural Engr., Pacific Elec. Ry. Co., Los Angeles, Calif. 

Berechnung statisch unbestimmter Siegmund Zentralblatt 
der Bauverwaltung, 1907, 


graphische Statik Heinrich Miller-Breslau, Alfred Stutt- 
gart, 1907, Vol. II, Pt. 162. 


gegenseitiger Martin Grining, Der 1921, 


Application Orthogonal Moments Problems Statically Indeterminate Structures,” 


Kine and Morkovin, Quarterly Applied Mathematics, 1944, 334. 
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equations are needed. The procedure somewhat analogous the method 
used for fixed arches and two-legged frames. this case, choosing the 
elastic center the structure the point application the hyperstatic 
quantities, the elastic equations containing one hyperstatic quantity only are 
obtained. the solution equations this type, containing one unknown 
each, slide-rule accuracy ample. 

complete analysis truss similar that the bridge under discussion 
and employing the methods previously explained given The 
base system obtained inserting three hinges the structure; one the 
crown the arch, and one each the second panel the end spans counted 
from the midpiers. The influence lines this base system are determined 
statics and utilized subsequently calculate those the hyperstatic structure. 


Rust,” Assoc. ASCE.—A statically indeterminate bridge im- 
portant the new highway bridge Dubuque, Iowa, with piers resting 
pile foundations, arouses interest about the amount and kind settlement and 
movement that may occur the piers. All foundations that depend piles 
for support are subject settlement. Appreciable settlement rotation 
the piers this bridge will introduce stresses the continuous trusses which 
may become critical, depending the amount and the nature the pier 
movements. Piles that can driven into the stream bed depth 
penetrate material which could vulnerable settlement. More de- 
tailed information about the nature the material the stream bed the 
pier sites would welcome: The soil’s sustain the pile loads, the 
method driving, the results pile tests, the probability scour these 
locations, facts about the soil support vertical planes the sides and ends 
the piers, the pile plans with spacing, and the actual sizes and pile loads— 
all would welcome information the paper about the Dubuque bridge. 

Structural engineers have believed that statically indeterminate bridge, 
especially one the size and the importance this one, should supported 
piers resting rock some other substance closely approaching rock 
its capacity sustain the loads without appreciable settlement. Although 
statically determinate simple span steel cannot specified with economy, 
under accepted working stresses, for length much more than 650 ft, canti- 
lever span perhaps some other type can specified obviate much the 
redundant effects from pier settlement. 

simple span two trusses, four supports are provided, one which 
redundant, because the portal and lateral bracing. The bracing between the 
trusses can designed lessen the effect that could result from unequal 
settlement the ends supporting pier. Settlements one pier with 
respect the others not affect simple and cantilever spans, such settle- 
ments affect continuous spans. planned observe the future the 
settlements the four piers which support the continuous trusses, may 
that the necessity for such observations will disregarded after some years 
small, inconsequential settlements. 


Statik des ebenen Tragwerkes,” Martin Julius Springer, Berlin, 1925, 370. 
Structural Engr., Chicago, 
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The writer has seen deterioration 50-ft statically indeterminate concrete 
arch spans due settlement their pile foundations, and has known the 
deterioration continuous steel structures from the nonrigidity their sup- 
ports. Allowable pile loads, due different specifications and the methods 
figuring the load sustaining capacity piles, have varied greatly that 
publication the research and field knowledge available from the Dubuque 
bridge piles would welcome. Those who have had direct charge the 
driving the piles must possess much practical knowledge the driving, 
which published would lasting value the engineering profession. 
The transfer the pier loads this bridge into the stream bed composed 
the implied material constitutes unknown discontinuity load transfer, 
which results from the imperfect and incomplete knowledge soil mechanics. 

The precedent has been set for number years sustaining continuous 
bridges pile foundations, and the precedent followed this bridge. 
seems that the capacity the piles sustain their loads without critical 
settlement should supported factual evidence well geologist’s 
advice, endeavor justify these foundations. 
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RATIONAL EXPLANATION COLUMN 
BEHAVIOR 


Discussion 
ERNEST ROBINSON 


Ernest Assoc. ASCE.—As any devotee the sport 
archery knows very well, there nothing uncertain mysterious about the 
behavior elastic column beyond the Euler length. 1915, and 1916, 
the writer published two analyses the curve whose ordinate everywhere 
proportional the curvature. The first equations for deflecting 
column connecting load and deflection both with the original length and the 
deflected length, and presented graphic solutions these equations, well 
charts for the required elliptic integrals. The second paper presented 
for the hydraulic catenary connecting the middle depth with the tension, area, 
wetted perimeter, and mean hydraulic radius well fairly large diagram 
showing the actual shape the curve. 

paper Frederick Ryder was published March, 1947, Proceedings. Discussion 
this paper has appeared Proceedings, follows: May, 1947, Robinson Rowe, and George Winter; 
and June, 1947, Robert Mains, and Chi-Teh Wang. 

Structural Engr., Turbine Generator Eng. Div., Gen. Elec. 

for Strength and Very Long Engineering News, June 10, 1915, 
pp. 1108-1109. 


Charts for Open Flumes with Catenary Ernest Robinson, ibid., March 23, 
1916, pp. 566-567. 
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DISCUSSIONS 


TESTS TIMBER STRUCTURES FROM 
GOLDEN GATE INTERNATIONAL 
EXPOSITION 


Discussion 


Carr, Assoc. has been added the 
knowledge timber design and construction this interesting study, which 
has already resulted further studies being made those who had advance 
information the results. The Timber Engineering Company, one the 
sponsors the San Francisco (Calif.) study, for three years has been carrying 
series tests full-scale trusses, has completed series tests ec- 
centric joints, and has developed, therefrom, empirical formula. addi- 
tion, has sponsored research the Massachusetts Institute Technology 
Cambridge the distribution stress around timber connectors. 

There are three treatments the paper with which the writer did not agree 
the time preparation. The first these that many places the 
paper the factors relating test load design load are given when the failure 
was not directly connected the design load used. such cases, the factor 
given represents only the minimum amount test load which the section 
the truss was submitted and does not indicate the maximum load might 
sustain. factor therefore meaningless. 

The second that the comparisons regard allowable connector values 
are based data published 1939, although the original designs were not 
prepared the basis these data, but data published 1935. 
believed that would have been more correct compare either with the 1935 
data with the recommendations published 1943, which were effect, and 
current practice, the time the report was preparation, and which continue 

The third that much the paper centers around so-called factor 
“tension perpendicular grain,” occurring particularly eccentric joints. 
The writer disagreed with this theory the time the report was process 


paper Committee the San Francisco (Calif.) Section, ASCE, Timber Test 
Program was published May, 1947, Proceedings. 


Timber Eng. Co., Washington, 
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TABLE 
SON BETWEEN 
Loap STRESS AND 
1943 RECOMMENDA- 


TIONS 
Rea 
Joint 
2 
(2) 1939 1943 


(a) 


TB-1 
TB-2 


3.86 
3.28 


2.95 
2.48 


TF | | 1.7 | 4.450 


(c) 


3.25 4.5 
TM-1 3.30 4.6 
2.53 3.1 
TM-1 3.3 4.6 
TM-2 3.6 4.9 


4.79 6.38 
J10 4.54 5.73 
3.12 5.70 
J12 4.08 4.55 
J13 5.06 
2.5 3.33 
J15 2.91 3.87 
2.93 3.96 
1.83 


(f) (Eccentric) 


shear formula. Con- 
nectors did not fail. Not true 
factors, grooves sheared and test 
setup failed. Tests not carried 
failure connector joints. 


CARR TIMBER TESTS 


Discussions 


and stated that felt the tests indicated and 
closely followed another theory, advanced 
Everett Lank, that the shear theory for notched 
beams more closely fitted the data. The authors, 
while stressing the perpendicular 
grain” theory, offer design approach. Mr. 
Lank has subsequently been proved correct 
series substantiating the shear theory, 
and empirical formula for design application 
based this theory has been developed, 
follows: 


width the section; and the effective depth 
the section, the distance from the inside 
edge the connector (nearest the unloaded edge 
of-the member) the loaded edge the member. 

Application these three comments the 
data presented brings the results given 
hereafter. order simplify the presentation, 
all values developed the paper for the 
1939 data and those based the 1943 data, pub- 
lished the National Lumber Manufacturers’ 
are given Table 26. 

The report states the first failure for truss 
TB-1 joints and L10. not clear 
just what the initial failure was. However, 
comparison current data gives the result shown 

The writer cannot reconcile Fig. with Fig. 
26, the description the French Building 
truss failure. Assuming Fig. and the de- 
scription correct, and applying the shear 
formula, the 1943 factor becomes 4.45 (Table 
26(b)) against the 1939 factor reported, based 
connector values (which did not fail) 1.7, 
appreciable increase. 

the Merry-Go-Round Building trusses the 
lumber failed and not the connector joint, hence 
the R-values are not true results. However, 


Eccentric Split-Ring Joints Yield Rational Design 


Formula,” Glass, Engineering News-Record, May 29, 1947, 
pp. 69-70. 

for Stress-Grade Lumber 
and its National Lumber Manufacturers’ 
Washington, C., 1944. 
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based current recommendations the factors Table become those shown 
Table 26(c). 

the discussion theoretical stress around connector (under the 
heading, Building (M) paper presented 
the 1947 Fall Meeting the Society tests conducted the Massachusetts 
Institute Technology, under the direction Prof. Albert Dietz, Assoc. 
ASCE, discusses data developed use stress coats and electrical strain 
gages. Current recommendations include factors developed these tests; 
connectors improved design have also been developed and are now available. 

the North Square Court truss tests the failure was due eccentric joints 


‘and hence the shear formula applies. Factors secured this method show that 


for joint the 1943 factor 2.98. Although they not reflect the true 
failure, the factors the 1943 data would correct the factors Table 
those given Table 26(d). 

Also, Table the changes factors, where applicable, based the 
1943 data would shown Tables (a), (c), and (d). The writer 
disagrees entirely with conclusion (see heading, Tests 
believed that eccentric joints may properly and capably 
handled with the shear formula. 

the basis the 1943 data, the factors Table should corrected 
shown Table 26(e). Based the shear formula application, the 1943 factor 
for test joint becomes 3.3. However, connector factors must considered, 
the factors are changed shown Table 26(f). regards the tests joints 
with ring connectors, the factors given Tables and should corrected 
the 1943 values, shown Table 

After eleven years’ experience, the writer feels that the comments con- 
clusions the tests joints with steel gusset plates are broader than the 
data available. 

However, view the results this program, this writer’s opinion that 
the designer who uses the 1943 recommendations, including the shear formula, 
assured factors which are more than safe when based these tests inas- 
much the smallest corrected true factor 2.98, whereas factor the neigh- 


borhood would considered sufficient experienced timber designers and 
research men. 


Assoc. ASCE.—The descriptions the equip- 
ment used and the methods testing are quite complete, and will valuable 
guides future tests full-scale trusses and joints involving eccentricity and 
loads angle the grain. most cases, the tests any one type were 
too few number and the variables too involved provide exact design 
procedure for similar types. The writer is, however, much pleased that this 
extensive test program has been published and particularly interested note 
that the complete report including test data, curves, field notes, photographs, 
and calculations filed several strategic points throughout the United States. 

The sudden failure under test loads the cantilever truss (TC) 
and the detailed examination joints other trusses again emphasizes the 


Prof. Mechanics and Acting Head, Industrial Eng. Dept., Univ. Florida, Gainesville, Fla. 
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need for competent engineering supervision timber structures the factor 
safety expected design provided actual construction. This 
same fact was brought out during World War when timber structures were 
built with little competent engineering inspection and investigations 
several failures showed that failure was due improper construction. 

The calculated for the Brazil Building (TB) trusses did not in- 
clude the bending effects the This fact was noted the report. 
However, since the stresses due bending could rather high, they should 
have been included the computations the The fact that the 
nature the failures the joints these and other trusses emphasizes the 


_importance secondary stresses significant. believed that many 


typical designs shallow timber trusses have not been checked for secondary 
strésses the past. The authors mention the 
TABLE need for additional tests joints which 
SON FoR shearing stress, direct stress, and 
BASED bending stress involved. This one the 
1939 1943 greatest needs the field timber engineering 
ALLOWABLE and additional experimental data combined 
Loaps stresses will aid materially the development 

more efficient and economical design timber. 
Although the data the lateral strength 
nails plywood are worthwhile and may used 
1943 1939 guide for further experimental investigations, 
are definitely needed. The tests covered the 
report not include wide enough range 
nail sizes, spacing, and edge distances 
predict the behavior nailed joints other than 

those tested.. 

noted the report, the buckling tests 
plywood panels were not exhaustive enough 
provide design details. The conclusion that stud 
spacing may increased beyond the customary 
in. in. seems feasible, but additional 
tests Douglas fir structural plywood should 
conducted along the lines investigated the 
United States Forest Products Laboratory 
Madison, Wis., aircraft plywoods. 

the joint tests the factor safety each joint was calculated the 
basis allowable loads from the TECO comparison 
the based later TECO are given Table 27. 

There are several questions that have arisen from the study this report 
and believed that the value the report will enhanced the answers 
these questions are given. The questions are: 


values over 1939 values. 


Connector Manual the Timber Eng. Co., Inc., Washington, C., 


Manual for Teco Timber Connector Manual the Timber Eng. Co., Inc., 
Washington, C., 1943. 
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What approximate method was used computing the stresses the 
members eccentric joint? 

How was the lower chord splice truss TM-2 reinforced before testing? 

Why were blocks added joint under test loads? Would this give 
true picture the type failure for this joint? 

Cannot some explanation for the large TABLE 


the statement under the heading, “II. PLywoop 
Joint Tests: Comments Tests Joints with 
tests the load was always parallel 
correet? Allow 
The values based the ultimate shear (1) (2) 
loads for various thicknesses Douglas fir 
wood, obtained the test results givenin 
report and the allowable stresses for various 480 1.22 
grades Douglas fir plywood, are given Table 
28. The report does not mention the grade 


Sheathing {575 360 
plywood used and the values Table may 


prove valuable comparison the R-values Puy, In. 
for the different grades. 


1.61 
tribution the literature and design data per- 2.02 
taining timber structures, well 


indication the need for more full-scale tests 
timber trusses, joints with combined stresses, 
nail tests, and structural plywood tests all 
trusses joints will behave exactly the same 1,090 360 
under load; but, more tests this nature are denotes de- 
conducted, the results can correlated with 


and 
un 
those contained this report. 


Stern,” Assoc. ASCE.—The report the committee contains 
much valuable information for timber designers and fabricators timber 
structures. interest learn what comments designers and fabricators 
the tested timber structures may make; whether their assumptions, theories, 
and computations were verified the test data, and which theories were 
contradicted. believed that additional valuable information could 
obtained from the designers and fabricators review their designs and 
construction logs could made available the basis the test data obtained. 

The fabricators probably can give the answer this question: Were split 
and shear plates, respectively, omitted purposely joints L7, L8, and 


truss TB-1, and joint because they were considered unnecessary? 


Associate Prof. and Head, Wood Constr., Virginia Polytechnic Inst., Blacksburg, Va. 
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Several the findings seem confirm the previously made assumption 
that, during initial loading and failures minor nature, adjustments take 
place the local and secondary stresses—often anticipated but not included 
detailed computations—which result reduction and more uniform distribu- 
tion such internal stresses and, thus, make the structure safer than was 
before adjustments took place. Although such adjustments may increase the 
load-carrying capacity the structure, the deformation and total deflection 
may also increased beyond those computed. Rolling shear the core 
within split rings may often accounted for result such adjustments. 

hoped that the appropriate industry will sponsor performance some 
the additional tests order secure the needed information. 
cannot expected that engineers architects will employ timber structures 
long some the basic design data are lacking and not available them. 
Industry may reminded that several United States technical institutions 
and universities are willing assist industry obtaining such vital data 
assigning specific problems graduate students interested timber structures 
assigning their experienced staff industry-sponsored projects. 

The statement made that member truss was replaced before 


test performance. According Fig. 25, member such designation does 
not exist. 
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